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-roposed Method of Loss Measurement 


for Semiconductor Rectifier Equipments 


I. K. DORTORT 


MEMBER AIEE 


T HAS long been recognized that the 

efficiency of a large rectifier unit cannot 
e accurately determined by input-output 
ower neasurements. Errors of measure- 
nent may be as great or greater than the 
osses in high efficiency equipments. 
[hese errors are aggravated to an un- 
mown degree by the distortion of cur- 
ent and voltage waveshapes inherent in 
ectifier operation. 

Because of this and other difficulties, 
he accepted procedure is the ‘‘segregated 
oss’ method in which the losses in each 
major component of the equipment are 
determined separately by conventional 
methods; conventional, that is, for all 
parts except the rectifier itself. 

It required many years to achieve 
easonable agreement on the method of 
Jetermining the losses in a mercury-arc 
ectifier. Almost all of the power loss 


roltage drop, having only a minor de- 
Sendence on current. The average value 
of the “arc drop’ multiplied by the aver- 
age forward current is a reasonably ac- 
irate measure of the arc loss over a 
wide range of current waveshapes. The 
preatest difficulty was caused by the 
problem of current and power measure- 
ments in single-way rectifier circuits. 
The same methods cannot be applied 
0 semiconductor rectifiers, even though 
single-way circuits are not preponderant. 
The reason is quite apparent in Fig. 1, 
which shows a typical curve of instantane- 
us voltage and current values for a sili- 
on cell. The forward current loss is 
strongly dependent on the current wave- 
shape because of the high incremental 
esistance and cannot be measured by 
average voltage drop except in special 
eases. Rectifier loss measurements are 
of necessity made with the d-c terminals 
short-circuited. In consequence, the cell 
current form factor will be substantially 
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s caused by a relatively constant forward’ 


different from that obtained in normal 
operation. 

It is necessary, then, to apply a correc- 
tion to the losses measured in a short 
circuit. Some of the previously proposed 
methods rely on corrections computed 
from the forward voltage curve. 

It was quite logical to facilitate the 
mathematical work by approximating 
the curve of actual values by a straight 
line with slope Ra, intercepting the voltage 
ordinate at a ‘threshold voltage” Ey. 
Most of the proposed methods of loss 
measurement are based on this concept. 
However, the graphical approach to Rg, 
and Eis not straightforward. 

To start with, semiconductor rectifiers 
in the class under discussion contain many 
cells, often in the hundreds, It would 
not be feasible for the equipment manu- 
facturer to measure and plot the char- 
acteristic of each cell, so he uses the cell 
manufacturer’s published curves. This 
is true even where the same company 
manufactures both the cells and equip- 
ments. The published curve represents 
the average characteristic of many cells, 
measured at a specified junction tempera- 
ture. 

The spread between cells, even if 
they have been matched, is quite wide and 
their individual characteristics are not 
necessarily parallel. Let us assume, 
however, that the curve to be used is 
fairly representative of the batch of cells 
in the rectifier and that the average Ey,’s 
and R,’s of the cells are equal to the Eu, 
and R, of the average curve. 

First, the average voltage drop curve 
must be corrected to the normal operating 
junction temperature of the rectifier. 
It is no simple task to measure or com- 
pute the mean junction temperatures of 
many cells in a rectifier, especially if the 
current division is unequal and the cooling 
is not uniform. Still more difficult is the 
task of correcting the curve. The in- 
formation necessary to do this accurately 
is not generally available. 

Next comes the problem of matching a 
straight-line equivalent to the curve. 
The uncertainty introduced by the human 
element might be overcome by a rigid 
formulation of procedure. Right or 
wrong, this would eliminate uncertainty 
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and controversy. It should be noted, 
however, that the values of Ey and Re 
are not constant, being dependent on the 
load level. This is also illustrated in 
Bicwele 


There is still another problem. A 
large rectifier consists of more than just 
a conglomeration of rectifier cells. There 
may be heavy current buses, some carry- 
ing alternating currents, some carrying 
unidirectional currents. There will prob- 
ably be many bolted joints. Some of the 
buses may be tapered or nonuniformly 
loaded. There will probably be cell 
fuses, balancing reactors, and perhaps 
control reactors. Their J?R (current? 
resistance) and stray losses or, more 
rigorously, the total effective resistance 
must be added to the incremental resist- 
ance of the cells. 


By making rather dubious assumptions 
and incongruously painstaking computa- 
tions, it is possible to come up with a value 
of effective resistance somewhat better 
than an educated guess. In very heavy 
current rectifiers, the losses in these parts 
can exceed 25% of the cell loss and de- 
serve a more accurate determination. 


Direct measurement of losses or effec- 
tive resistance of these parts is generally 
not feasible. It introduces problems of 
instrumentation and requires replacement 
of at least half of all the cells in one recti- 
fier section with solid conductors having a 
small resistance compared to the incre- 
mental resistance of the cells. The choice 
between computation and direct measure- 
ment is a tossup, somewhat in favor of 
computation. 

Other methods have been suggested 
and tried for applying waveshape correc- 
tion factors designed to avoid these 
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PEAK FORWARD CURRENT 


Fig. 1. Typical instantaneous forward voltage 

characteristic of silicon cells at 100 C junction 

temperature and straight-line approximations 
at three loads 


A—Etn=0.92, Ra =0.00088 
B—Etn=0.85, Ra=0.00108 
C—Etn =0.80, Ra =0.00146 
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tedious and uncertain computations. 
These methods involve manipulation of 
loss values measured in two configurations 
of the circuit. The variation is produced 
by replacing at least one series element 
in every series connection with a solid con- 
ductor. If there are no series elements in 
the rectifier, they must be added. This 
requires additional connections not pro- 
vided in the rectifier, and the additional 
losses must be determined and in- 
cluded in the loss calculations. 

In large rectifiers involving hundreds 
of cells these methods require many 
hours of labor and introduce errors due 
to temperature changes in various parts of 
the rectifier. This is especially true in 
fuses which are often connected directly 
to the cell leads and are heated by them. 

Note that it is not permissible to use a 
single, heavy-current jumper to short- 
circuit a block of parallel cells because to 
do so would completely change the 
geometry of the copper (Cu) and the Cu 
losses. 


IEC Proposal 


During the 1959 meetings of the Semi- 
conductor Subcommittee 22-2 of the 
International Electrotechnical Com- 
mission (IEC), a Task Force was set up 
and instructed to come up with an accept- 
able means of measuring rectifier losses 
before the end of the Paris sessions, so 
that the contemplated document of IEC 
recommendations could be completed and 
issued for balloting in 1960. The pro- 
posals submitted by the several nations 
were reviewed and a recommendation 
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Fig. 2. Correction curves for deviation from 
sine wave in short circuit 
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prepared. Towards the end of these 
deliberations, the United States delega- 
tion submitted another proposal, in- 
volving two loss measurements without 
any changes in the configuration of the 
circuit and without the need for deter- 
mining Ey, and Ry. There was not 
sufficient time for careful study and 
evaluation, but Dr. J. C. Read, head of the 
Task Force, was interested. At the last 
meeting of the Subcommittee, he pre- 
sented the group’s formal recommendation 
as well as the new proposal. He recom- 
mended that the new scheme be studied 
and developed by correspondence be- 
tween us, after the close of the Paris 
meetings, for possible inclusion in the 
document at a later date, A motion was 
made and carried to put this plan into 
effect with the Secretariat in Sweden a 
cognizant party. Our work culminated 
in the simple system defined by equation 
7 of this paper, and has been included in 
the latest revision of the IEC Recom- 
mendations on Monocrystalline Semi- 
conductor Rectifiers. 


Basis of Proposed Method 


Within the accuracy of the best straight- 
line approximation, the forward voltage 
drop of the cells can be written: 


v=EntRei (1) 


The forward current power loss of a com- 
plete rectifier structure, including stray 
losses, can be expressed in the form: 


P=A Taye: Eqy+ 
2ms(bRg+cR+dRot+ cane )) 
=A-lavg En +B-TrmsR’ (2) 


where 


Ey,=threshold voltage of cells 

Ra = incremental resistance of cells (straight- 
line) 

Tayg= average current of cells 

Ry, Ro, ...=effective resistance of uni- 
directional current and alternating 
current carrying conductors and 
components 

R’=total effective resistance, including R, 

A, b,c, d, B are circuit constants 

P=power loss of rectifier under a particular 
mode of operation 


Equation 2 is quite obvious and fine 
as far as it goes, but exceedingly difficult 
to use for reasons given previously, 
Furthermore, the circuit constants as well 
as the R’s are influenced by current dis- 
tribution, temperature variations, etc. 
Logically, the entire job belongs to an 
analog computer. In the proposed sys- 
tem the rectifier itself is the analog com- 
puter and is exactly tailored to the 
requirements. 

When a rectifier is operated under nor- 
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mal conditions of d-c output voltage a: 
current, the cell currents will usual 
have a different waveshape than whent 
is delivering the same average direct erp} 
rent in short circuit. The average vali 
of each cell current must be the sany 
except for possible changes in curte 
distribution. The average value of 
the cell currents must be exactly the sami 
The rms values of the cell currents we 
generally be quite different. If Jims is tH 
value in short circuit, the value in normi 
operation can be defined as KIims. TY 
power loss in normal operation is os 
written: 


P=(A Lave Eu )+K*(B-Trms'R’) ( | 


If the input power to the rectifier - 
short circuit is measured at two level 
rated direct current Jg and KIq, respel 
tively, then 


Py =(A-ToveE,)+(B-Pmme'R’) mG 
P2=K(ATnve'Em)+K°¢B-T'mmeR!) 
By elimination of the E wynand R’ terms: 


ple ae (a 
ie } 
From equations 2 and 3 it can be see 
that K is the ratio of rms current value 
in normal and short-circuit operation. 1] 
is therefore also the ratio of the form fa 
tors of the two waveshapes in these mode 
of operation and can be so defined. 
Assuming ‘rectangular waveshape' 
without overlap, for normal operatio 
and perfectly sinusoidal waves in shov} p| 
s 


circuit, Pa 


a: 

K=1.1 for 6-phase double-way and douliiy 

Y rectifiers 

K=1.0 for single-phase full-wave rectifier y 
resistive loading 

K=0.9 for single-phase full-wave rectifier 

inductive loading q 


In general, single-phase rectifiers fe 
ing resistive or back electromot 
force (emf) loads will, in normal operatit 
have cell currents approaching sinusoi 
rather than flat-topped waves. 
these, K=1.0 and, obviously, only ont 
measurement is required. 

When measuring the losses in heg 
current rectifiers, the short-circuit 
may account for a significant percenté 
of the total. This loss will naturally 
included in the power loss derived fr 
P, and P, and must be subtracted. If 
most easily determined by measuring 
voltage drop V’ across the short-cire 
bar, and then 

Gage! 
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In most cases JgV’ can and should be 
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ig. 3. Test circuit no. 1 direct measurement 


made negligibly small. If the ammeter 
hunt is not part of the rectifier, its loss 
hould be subtracted. 


ORRECTION FACTORS 


The validity of neglecting overlap 
equires substantiation. It is not enough 
hat this procedure simplifies things, or 
ven that it falls in line with existing 
tandards for rating rectifier transformers 
nd measuring their losses. Let us 
sume a typical 6-phase rectifier with a 
ommutating impedance of 10%. The 
ngle of overlap, wu, is then 26 degrees and 
he rms value of cell current in normal 
peration takes a correction factor for 
pverlap of 0.973.! K should theoretically 
pe reduced from 1.1 to 1.07, causing a 
eduction of approximately 4% in Pp. 
he net result is a reduction in P of less 
an 2% in the usual case. Two per cent 
-a 2% loss is truly negligible and little 
nough margin for the undetermined 
mcrease in stray losses produced by the 
teeper wavefronts of normal operation. 
0 correction for overlap is recommended. 
One other error must be considered and 
valuated. K has been assigned values 
based on rectangular cell currents in 
ormal operation and pure sine-wave 
urrents in the short-circuit tests. Most 
lectifiers do draw essentially sinusoidal 
trents in short circuit, but not all. 
pecially in low-voltage high-current 
ectifiers the waveform may be quite 
istorted. The most _ straightforward 
olution is to determine the true form 
actor of the cell current in short circuit 
md adjust K accordingly. It will be 
ecalled that K has been defined as the 
atio of the form factors of the cell cur- 
ents in normal and in short-circuit 
peration. 

Standardization bodies may prefer to 
ise fixed K factors and apply corrections 
0 the result only in special cases. There- 
ore, the following procedure is suggested 


as an alternate. Fig. 2 provides correc- 
tion factors to be applied to the normally 
obtained power loss if the form factor of 
the short-circuit current deviates too 
much from 1.11 in 6-phase double-way and 
double-Y rectifiers. There seems to be 
no need for similar corrections for single- 
phase rectifiers since deviation from si- 
nusoidal currents is not expected. How- 
ever, similar curves can be provided if re- 
quired. 

Let F’ betheactualform factor, P;’ and 
P,' the actual loss measurements, and 
let P’ be computed from them, using 
the standard K. The equation for the 
correction factor is then: 


P—P’ a —K*P,'/Py'—K) [i —(.11/F")?) 
jo (R91) Py) [Py RAR =D) 
(8) 


A +5% band of maximum allowable 
error has been drawn in Fig. 2 as a basis 
for discussion. Any rectifier falling with- 
in this band would have its losses com- 
puted directly from equations 6 and 7 
without correction. It is believed that 
most large rectifiers will fall inside this 
band. 


Loss Measurements and Circuits 


RECTIFIER TEMPERATURE 


The normal rectifier loss falls some- 
where between P; and P2, which can 
differ by no less than 11% or more than 
21% of P;. During the P. measurement 
the loss in the cells will generally be no 
more than 6% above normal operating 
loss, even though the output current 
has been increased 10%. All other 
losses will be normal or slightly below. 
Therefore, it is recommended that P: be 
measured at normal rectifier temperature 
with the rectifier carrying KJg The 
current should then be dropped to Ig 
and P, measured as qtickly as possible 
before the temperature of the Cu buses 
and heat sinks can change appreciably. 
It is further suggested that the tests be 
conducted with the rectifier adjusted to 
the temperatures that would be obtained 
in normal operation with an ambient 
air of 30 C (degrees centigrade) or in- 
coming water at 25 C, depending on the 
system of cooling employed. 


CONDUCTION PERIOD GREATER THAN 180 
DEGREES 


High-current low-voltage rectifiers, 
single-way rectifiers, and rectifiers em- 
ploying control reactors in the cell circuits 
may have a large part of the total com- 
mutating reactance in their unidirectional 
branches. In such cases, the conduction 


period in short circuit may be consider- 
ably greater than 180 degrees. The suc- 
cessive lobes of the alternating current 
will not be exact replicas of the cell cur- 
rents, since those portions beyond 180 de- 
grees will be cancelled either directly or 
in transformation. Both alternating and 
cell currents should be checked oscillo- 
graphically, and form-factor corrections 
made, if necessary. 


COMMUTATING REACTANCE 


Reactance in unidirectional circuits is 
difficult to measure. If large, it becomes 
an important factor in the regulation, 
power factor, and ripple output of the 
rectifier unit. Voltmeter, ammeter, and 
wattmeter measurements in the follow- 
ing described tests may prove useful in 
determining the total commutating react- 
ance of the rectifier. When conduction 
greater than 180 degrees takes place, the 
conventional methods of determining re- 
actance from the instrument readings 
will definitely not apply. A satisfactory 
method is still to be developed. 


Test Crrcuir 1, Fic. 3 


The most desirable way to measure 
rectifier losses is unquestionably by direct 
wattmeter measurements at the input 
terminals of the short-circuited rectifier. 
For large high-voltage rectifiers, standard 
wattmeters can often be used with- 
out instrument transformers. Industrial 
rectifiers will almost invariably require 
current transformers and low-voltage 
wattmeters. Fig. 3, illustrating a typical 
circuit, is hardly required, but is in- 
cluded for the sake of completeness. 
The unusual representation of the a-c 


Fig. 4. Test circuit no. 2 primary wattmeter 
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(A) 


(B) 


Fig. 5. Test circuit no. 3 primary-secondary wattmeter 


A—Y-connected current transformer 
B—Delta-connected current transformer 


voltmeter and ammeter connections is 
used to indicate that all 3-phase currents 
and voltages are to be measured by any 
suitable means. 


Test Circuit 2, Fic. 4 


Single-way rectifiers do not lend them- 
selves to the above method, although in- 
strumentation has been developed for 
mercury-arc rectifiers and can be adapted 
for semiconductor rectifiers. Moreover, 
both single- and double-way rectifiers are 
being built in current ratings and physical 
configurations that preclude the insertion 
of current transformer without signifi- 
cantly changing the geometry of the a-c 
connections and, consequently, their 
I’?R and stray losses. For these and 
other reasons, it may be desirable to 
measure the power input to the rectifier 
and rectifier transformer combined at the 
a-c line terminals of the transformer. 

The transformer losses corresponding to 
P; and P; must be accurately known and 
subtracted from these measurements 
before they are combined to obtain P. 
It should be remembered that when the 
rectifier is delivering rated direct current, 
Ia, in short circuit, the primary current 
is I;/K, and when the rectifier delivers 
KIa, the primary current is I,, I, being 
the rated line current (neglecting overlap). 
Because of possible deviation from exact 
values all input currents should be re- 
corded during loss measurements. Trans- 
former losses should be measured at the 
currents and temperatures corresponding 
to the P; and P» tests, or at values close 
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enough to permit corrections to these con- 
ditions with minimum error. 

The primary voltage during these tests 
will be somewhat higher than the imped- 
ance voltage of the transformer and will, 
in many cases, require potential trans- 
formers. The power factor may be quite 
low, and wattmeters of the same class as 
used for transformer loss measurements 
are recommended. 


Test Circuit 3, Fics. 5(A) AND (B) 


The preceding test circuit requires 
much extra labor and time in calibrating 
the rectifier transformer. With simple 


Fig. 6. Rectifier of 2,800 kw on which first test was made. + 


double-way rectifiers, it is. possible t 
combine the advantages of test circuit 
1 and 2 and eliminate the transforme 
losses from consideration. This is 
by connecting the current coils of the watt 
meters to current transformers in th 
primary and the voltage coils directl| 
across the transformer secondary termi 
nals. This circuit is shown for a delta 
delta transformer in Fig. 5(A). Fig. ab 
shows a Y-delta transformer with the 
rent transformers connected in delta t 
transform the line currents to the 7 
phase angles and waveshapes as the alter: 
nating currents of the rectifier. 7 
In this test, the rectifier transfo me t 
operating at very low induction, wir 
draw a negligibly small exciting curr nt tt 
It becomes, in fact, a large current trans 
former. 


Test Results 


2,800-Kw RECTIFIER 


The first test we had the opportunit‘ 
to make was carried out on an 8,000) 
ampere 350-volt single-bridge water 
cooled rectifier, close-coupled to a delta 
delta oil-immersed transformer; see Fig 
6. The loss measurements were made ox 
the primary side in accordance with tes: 
circuit 2. Although the cell buses ar: 
water-cooled, the bulk of the Cu, cc 
vection-cooled and operating at 75 C, ini 
curs most of the losses. Two sets o 
measurements were made, the first “ 
the main Cu at 24 C and estimated june, 
tion temperatures of 54 C; the se 
test with main Cu at 75 C and estima 
junction temperatures of 84 C. 
test results are presented in Table I 
matter of interest. 


oe pete aia: 
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Table 1. Loss Measurements* 


Transformer Calibration 


Test 1 Test 2 
Top Oil Top Oil 
Temperature, Temperature, 
29C 55 C 


Loading, Ia 


110% 


Loading, Ia 


100% 100% 110% 


Cu loss, watts. ..21,170. .25,620. .21,870. .26,470 
[ron loss, watts. . 100 100 100 100 


Total loss..... 21,270. .25,720..21,970. .26,570 


Loss Measurements 


Test 1 Test 2 
Pi P, P; P, 

Total loss, 

BAECS os. 56 cers 67,060. .79, 130. .69, 900. .81 ,600 
Measurements, 

3 Se 45,790. .53,410. .47,930. .55,030 
Rectifier Loss, 

P, watts 51,750 52,250 


Estimated Rectifier Losses 


Test 1 Test 2 
Estimated 
Junction Tem- 
perature, 54 C; 


Estimated 
Junction Tem- 
perature, 84 C; 


Estimated Estimated 
Average Cu Average Cu 
Temperature, Temperature, 
24C 75 C 
Sell losses, watts....44,500.......... 42,800 
Ampere-trap fuses, 
se. ice dg CORN co asad aefsv'audl ae « 1,820 
BPIOSSES... 2... es. Se OR ers eke teas ee 4,370 
Estimated total 
HOSS cndrs ab vaverafoyics AD OBO iyestesia> 48,990 
Seto, watts........ ee TEs acts 8 +3,260 
BERNE VG slisic ota ele 2 SRB oer ds See +6.7 
3 8,000-ampere 350-volt water-cooled rectifier, 


widge-connected; close-coupled to delta-delta oil- 
mmersed self-cooled transformer; circuit 2. 
Wimary wattmeter. 


Regrettably, our lack of experience, 

anning, and time detracted from the 
wecuracy of calibration of the transformer. 
osses were measured with currents and 
emperatures different from the values ob- 
ained during the rectifier loss tests. 
Sorrections were computed within the 
imits of accuracy imposed by existing 
echniques. Small errors in transformer 
Isses can be cumulative and produce a 
arger error in the final result. It is one 
f the weaknesses of the primary loss 
jeasurements method. 
This may or may not be the explanation 
f the 61/,% “error.” Actually it should 
e remembered that this error is the 
ifference between a loss computed from 
leasurements by a method we are trying 
) prove and an estimated loss based on 
mplifying assumptions. Even so, 6'/2% 
fa 2% loss is not intolerable. 

Measurements were made at two esti- 


lated junction temperatures to see how 
ee: 
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Table Il. Loss Measurements* 
Estimated Rectifier Losses 

Estimated junction temperature, C......... 87 
Estimated cell losses, watts................ 4,950 
Aiip- trap fuses;wattse.ets sacs ne Annas aoe 30 
Cu losses at 50 C average, watts............ 50 
Shutitewetesins br cwinis: arace ca tetas Coens 80 

Potal loss "wattsiiices-ccicdatcuacecensn 5,110 

Circuit 1: Direct Measurements 

PisP ep eek ss cet cated ee 4,800..... 5,440 
Rectifier loss vite ivaon terre aac... 5,110 

SILO Semrnaynrs cis taetaversteelere sated Lene 0 

Circuit 2: Primary Wattmeter 
P, Pe 

PPotal watts. amiccs stata Geceesi ine any ip POW Soar 8,200 
Wranstormervioss:.77 tac eri coke 2528055 ae 2,710 
dog? a Bs eee hile eae AA oa tee Pe 4,920..... 5,490 
Riechifiertlassy Ps > sk a os. outon as siete 5,070 
FSET OR WALES a parc tesasdions © 4 arr coerouarersr dane —40 
ERROR. Voc ans ose te ee eee -—0.8 


Dates St Yor ha AEE aR ane eters 4,840..... 5,370 
Rectifier loss.bF tyes stecas ron indore mick 4,910 
ISfrory watts. dancs weet clare enioeaes — 200 
B Sho Rare EIAs ATTN GR ARQIROS ORO CICERO —3.9 


* 2,000-Ampere 100-volt forced-air-cooled rectifier, 
bridge-connected; integrally assembled with delta- 
delta class B transformer. 


well Ey, and R’ could be determined, 
working backward from the measured 
losses. By solving equations 4 and 5 
and substituting V3 Tave for KIrms, 
the following is obtained: 


KeP, =P; 
Lj 9 
ORC A Te (9) 
KP,—K?P 
et lead oe al (10) 


Si eG ere 


Substituting suitable values for Javg, A, 
and B, and from Table I for P, and Po, 
we arrived at the following results: 


Eth, Volts R’, Ohms 
54 C Junction, 24 C Cu..... Osbviaeine ea 0.0023 
84 C Junction, 75 C Cu..... ORS 4S torenate 0.0031 


Both Em and R’ have changed in 
opposite direction to that expected from 
increasing temperatures. It is obvious 
that the loss measurements were not 
sufficiently accurate to permit this reverse 
procedure, since we are dealing here with 
the small difference of two large, nearly 
equal quantities. 


200-Kw RECTIFIER 


’ We were able to test all three test cir- 
cuits on a 2,000-ampere 100-volt forced- 
air-cooled rectifier employing a delta- 
delta class B transformer; see Fig. 7. 
The rectifier was bridge-connected. The 
results are shown in Table II. 

We take no credit for the exact co- 
incidence of estimated and measured losses 
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Fig. 7. Rectifier of 200 kw which was tested 
in all three circuits 


in test circuit 1. Had it not occurred, it 
would still seem reasonable to use the 
direct measurement as the basis for com- 
parison. P:/P,;=1.13 and P is inter- 
polated between them; any error cannot 
be much greater than the error of meas- 
urement. 

It is significant that despite the low 
secondary voltage of approximately 2.2 
volts available for the wattmeters in 
test circuits 1 and 3, the maximum spread 
for all three measurements was less than 


4%. 
Conclusions 


By interpolation between two measure- 
ments which are generally 15% and never 
more than 21% apart, the proposed 
method of measurements achieves results 
with minimum reliance on human judg- 
ment or unreasonably difficult computa- 
tions. The results obtained in our tests 
are gratifying, and it is hoped that others 
have been able to achieve similar or 
greater success. 

The testing of all three circuits on one 
unit has been particularly rewarding. 
However, it is recognized that a 2,000- 
ampere unit does not present the prob- 
lems encountered in heavy current recti- 
fiers. Nor do the results obtained sup- 
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port entirely our order of preference of 
the three test circuits: (1) Circuit 1: direct 
measurements. (2) Circuits 3A and 3B: 
primary current, secondary voltage. (3) 
Circuit 2: primary measurements. 
Accurate calibration of a rectifier trans- 
former is difficult and should be avoided 
if possible. On the other hand, single- 
way rectifiers do not lend themselves 
to testing in circuits 1 and 3, and 12-phase 
rectifiers require special consideration. 
A carefully calibrated test transformer 
would be highly desirable and would save 
much time. However, it would not be 
feasible to build one that would cover a 
wide range of rectifiers, nor would it solve 
the problem of ultracompact assemblies. 
The simplest solution seems to be the 
combined measurement of transformer Cu 
losses and rectifier losses without segregat- 


ing them. Existing and tentative stand- 
ards would require changing and a con- 
siderable saving in testing time would be 
gained. 

If the standards are to be relaxed, a 
study must first be made of the changes in 
the procedures for determining the trans- 
former temperature rise as well as the 
correction of transformer losses to a 
standard temperature. It might be 
pointed out that the standard temperature 
of 75 C for Cu losses seems to be an 
anachronism for class B and class H 
transformers. 

Let us test the combined loss method 
on the measurements made on the 2,800- 
kw rectifier. Substituting the total P; 
and P: measurements of Table I directly 
in equation 6, combined rectifier and 
transformer losses of 78,900 and 77,300 


Discussion 


L. F. Borg (Allmanna Svenska Elektriska 
Aktiebolaget [ASEA] Ludvika, Sweden): 
As secretary of the IEC Committee on 
Monocrystalline Rectifiers and Rectifier 
Equipments, I would like to express great 
satisfaction with the test method which Mr. 
Dortort has described in his excellent paper. 
A fundamentel principle in the LEC recom- 
mendations, which will be sent out for 
criticism, is that most measurements on 
large rectifier installations should be made 
on the separate units and normally as no- 
- load and short-circuit measurements. It 
was early felt that a method for short-circuit 
loss measurements on the rectifier assembly 
was badly needed and several more or less 
complicated methods have been discussed. 
The method described in the paper was pro- 
posed in order to make the measurements 
independent of any previous measurements 
on the data of the rectifier cells, which was 
an inherent weakness of the other methods 
discussed. So far as can be seen now, Mr. 
Dortort’s method is simple, practical, and 
sufficiently accurate for the purpose. 

Mr. Dortort’s comments on the impor- 
tance of overlap are essential and I wouldlike 
information on one further item in that 
connection. It seems that with bridge 
connections the short-circuit conductor to 
be applied might have an inductance large 
enough to cause a circulation current com- 
ponent not appearing in the a-c circuit. 
That is particularly likely in such installa- 
tions where the plus and minus terminals 
are situated as far apart as possible to 
minimize the risk of short circuits. If that 
is the case I think it might be valuable to 
introduce a method to correct for the errors 
caused by that phenomenon. Mr. Dor- 
tort’s observations on this point would be 
appreciated. 


H. Winograd (Allis-Chalmers Manufactur- 
ing Company, Milwaukee, Wis.): Mr. 
Dortort has rendered a valuable service by 
presenting this excellent paper at a time 
when the American standards for semicon- 
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ductor rectifier equipments are in prepara- 
tion. The measurement of losses in semi- 
conductor rectifiers has been under discus- 
sion for several years, in connection with 
international and American standardization 
work. Mr. Dortort had a leading part in 
developing the proposed method. 

While the paper deals with measurements 
for rated rectifier current, the same method 
could undoubtedly be applied for other 
currents, such as 75% or 50% of rated 
current. 

In drafting the standards for loss measure- 
ment, consideration should be given to 
possible deviations from the required test 
conditions. It might be difficult or im- 
possible to obtain the required operating 
temperature, particularly for air-cooled 
assemblies. Permissible deviations from 
the required temperature, or some reason- 
able correction factors, should be specified. 

In some cases it may be difficult to obtain 
measurements at specific current values 
required for the application of equation 6 
of the paper because of limitations in the 
current-regulating equipment or variations 
of the supply voltage. It would be desir- 
able to include alternative equations for 
measurements made at two current values 
related by a ratio other than K. 

For rectifiers in which the conducting 
period of the cells during the test is con- 
siderably greater than 180 degrees, it was 
suggested in the paper that form-factor 
corrections be made from an oscillographic 
check of the a-c and cell currents. An 
oscillographic test should be avoided, if at 
all possible, because of the added complica- 
tions. Perhaps the power factor deter- 
mined from the loss-measurement data, or 
the length of the conducting period meas- 
ured with an oscilloscope connected across a 
unidirectional branch, could be used for 
applying predetermined correction factors, 
even if they could be considered only 
approximate. This subject requires further 
study. é 

No mention is made in the paper of the 
losses caused by the reverse voltage of the 
cells in rectifier operation. In the proposed 
international standards, the reverse losses 
in the cells are assumed to be negligible. 
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watts for the high and low temperaty 
tests are obtained, respectively. ] 
tracting the transformer Cu losses 
110% Ia), the respective rectifier lo; 
are then 52,330 and 51,580 watts, aln 
the identical values derived in Table 

If the combined method is adop 
circuit 1 then becomes the prefe 
method, providing the greatest ove 
accuracy and requiring the least total 
ing time for combined rectifier and tr 
former Cu loss measurements and 
runs. The difficulties of separate transi 
former Cu loss measurements and calibra 
tion are eliminated. 


Reference 


1. Mercury Arc Powpr RectiFiers (b 
O. K. Marti, H. Winograd. McGraw-Hill 
Company, Inc., New York, N. Y., 1930. 


The losses in voltage-dividing resistors ani 
surge-suppression devices are to be deter 
mined by calculation. Guides for calculat 
ing these losses should be included in thi 
standards. 


J. C. Read (Associated Electrical Industr 
Ltd., Rugby, England): I think the p 
timely and valuable; it shows how cor 
cated this problem can be, which at 
seems so simple. Nevertheless, in spi 
the fact that I had a small part in pu 
this proposal to the IEC, I now belie 
still requires further consideration. It 
probably be clearest if we divide the subjec| 
into sections, as follows: i 
1. The suggestion that the standard 
should call for the transformer Cu loss 
rectifier loss to be measured toge 
without segregation, seems to me ope 
several objections. It would not affect 
errors discussed under sections 6 and 
below; the cost and time needed to bri 
the transformer winding up to 75 C wot 
be objectionable; and, above all, it we 
be inconvenient where the transformer 
rectifier are built in different factories, 
quite unsuitable in the not infrequent ¢ 
where they are built by different maker} 
or in different countries. I think it i] 
indispensable to have a method for dete 
mining the segregated loss in the rect 
This, it would seem, can only be done 
testing the rectifier with a suitable t 
former. The transformer used can eithe 
a test plant transformer or the rectifier’s owl 
transformer, whichever is convenient. 
2. Asarough allowance for the harmonii 
stray loss in the rectifier it has been p 
posed that the rectifier losses determined 
the test shall (as for the transformer) 
those that would exist with rectangular ¢ 
rent waveforms giving the same averag 
value of d-c output current. 
For polyphase rectifiers, this is the sai 
as assuming that the real current waveform: 
in service will increase the rectifier Cu loss 
over that which occurs with sinusoidal ¢ 
rents, by about 6% of the combined rectifie 
Cu loss plus resistive component of cell lo < 
This allowance will be too large in sot 


pases and too small in others. However, I 
think that, like the corresponding allowance 
ade in determining the transformer Cu 
os, it is fully justified on the grounds that 
a) it greatly simplifies the treatment; (0) 
ol error will be small compared with the 
otal loss; and (c) above all, it does at least 
enable tests on different rectifiers to be 
compared on a uniform basis. 
However, for single-phase rectifiers the 
surrent waveform in service, which in this 
ase depends on the nature of the d-c load 
° be supplied, may not be accurately 
own at the time of the test and an arbi- 
rary assumption must therefore be made. 
To assume rectangular current waveforms, 
as is attractive on the grounds of simplicity, 
would be to assume rms currents actually 
.. than those likely in service and would 
thus be equivalent to making a large nega- 
tive allowance for harmonic stray loss, which 
would not be desirable. However, in many 
Important cases, including that of traction, 
the current in service will not be unduly far 
irom sinusoidal. For these reasons I be- 
lieve that for the purposes of standards it 
will be best to base the efficiency deter- 
mination on the current in service being 
sinusoidal. This means that in the single- 
phase case nothing would be allowed for 
harmonic stray loss. I think this is per- 
missible, because in practice there are 
scarcely any cases where the exact value of 
the efficiency is very important for single- 
phase rectifiers, but the main need is rather 
for accurate comparative figures. 

3. For the foregoing reasons the common 
oractice of basing transformer rating and 
osses on the rms values of rectangular cur- 
‘ents giving the specified direct current is 
indesirable in the single-phase case but 
these should be calculated on the basis of 
sine-wave currents. 

4. There was great anxiety to obtain a 
nethod of loss measurement that would 
ake account of the shape of the forward- 
yoltage-drop characteristic that is shown 
ypically in Fig. 1 of the paper. This led to 
he suggestion that the short-circuit power 
nput be measured at two different values, 
esulting in equation 7. I was at first much 
ittracted by this elegant solution, but after 
utther experience with it I now think better 
cculracy can be obtained in a simpler way. 
Phe difficulty with equation 7 is that the loss 
s calculated from the difference between 
wo comparatively large quantities, both of 
vhich have to be measured at low power 
actor and, therefore, with rather low 
;ccuracy. For this reason, with normal 
tandards of testing the result can easily be 


Table Ill. Values of Rectifier Loss 
Cu Loss Cu Loss 
0.1 X Cell 1.0 X Cell 
Loss Loss 
nput power P; at 
oo, GA a ae 100 100 
nput power P2 at 
MMLC aio hah crore o sede vaya 1.1X45.45-+ ..1.1X25+ 
: 1.21 54.55 b. 20K 75 
=115.9 =118.2 
oss P given by 
equation 7 of 
DIOL ooo) nyct a te, -F ase r LOMA asin be oe 115.7 
Iternatively: 
Loss P given by 
a single measure- 
ment at 1.08 74..112.7..... raeeter 114.5 


ULY 1961 


in error to the extent of 6% or more when 
the methods of equation 7 are used. 

Assume for simplicity that Uy, equals 0.5 
of the total cell drop at peak current. Now 
in our experience, with large rectifiers the 
loss in the connections, fuses, etc., in the 
rectifier cubicle is nearly always between 
extreme limits of about 10% and 100% of 
the total cell loss. Then, in the absence of 
instrumental errors, we obtain the extreme 
values of the rectifier loss shown in Table 
III. From this table it is obvious that a 
single short-circuit test at 1.08 Iz, which is 
slightly simpler to do, would give actually 
greater accuracy than equation 7. 

5. I think that for the purpose of a 
standard specification it is best to take the 
losses in the rectifier cubicle as those 
measured at the temperature normally pre- 
vailing in the maker’s works, since this gives 
an inexpensive test and a good comparison 
between alternative rectifiers. Correction to 
another temperature would be questionable 
if calculated and costly if done experi- 
mentally. 

6. In nearly all cases the input to the 
rectifier cubicle cannot be metered directly, 
but must be metered with the aid of a main 
stepdown transformer. With double-way- 
connected rectifiers, this presents no undue 
difficulty. Hither (a) the power input to the 
short-circuited rectifier and transformer is 
measured on the primary side, and then the 
loss in the transformer and interconnecting 
leads is measured by short-circuiting at the 
input terminals to the rectifier cubicle and 
measuring the input at the same primary 
current; or (b) for a 3-phase bridge connec- 
tion the input to the short-circuited rectifier 
is measured directly as in Fig. 5. Either 
way there are no substantial inherent errors. 

However, many semiconductor rectifiers 
are connected single-way (e.g., 6-phase 
double-Y or equivalent), and with the con- 
stantly increasing voltage rating of the cells 
available the use of single-way connections is 
‘increasing. Here we must meter on the 
primary side of the transformer, as in Fig. 4, 
and then measure the losses in transformer 
and leads as in (a) above (except that only 
particular terminals are short-circuited in 
this latter test, in accordance with the well- 
known rules!), However, in this case the 
presence of reactance in the secondary cir- 
cuit can produce substantial and inherent 
errors, particularly with low-voltage recti- 
fiers. 

The assumption made so far has been that 
the current waveform in the short-circuit 
test will be as shown in Fig. 8. With 
reasonably high total reactance the slight 
distorting effect of Uy, is negligible and these 
waveforms are sufficiently closely ap- 


_proached, provided that the secondary 


reactance is negligible. However, with 
single-way rectifiers there will usually be 
appreciable secondary reactance, partly be- 
tween the two opposite secondary phases 
180 degrees apart, partly in the intercon- 
necting leads, and partly in the internal 
connections in the rectifier cubicle. In the 
extreme case where the secondary reactance 
is the determining feature, the current 
waveforms (with a perfect rectifier) become 
as shown in Fig. 9. Practical cases usually 
lie between the two extremes represented by 
Figs. 8 and 9, but consideration of Fig. 9 
will show the magnitude of the errors that 
may exist. These are due to two causes. 
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SECONDARY PHASE 
CURRENT, RMS = 


aq , WHERE Ig = 


DIRECT CURRENT 
VALUE. 


PRIMARY PHASE 
CURRENT, RMS= 


TT : 
6a2 Iqg,for 13 


PRI./SEC. TURN 
RATIO 


Fig. 8. Rectifier current waveform in short- 
circuit test when secondary reactance is absent 


If we regard the primary and secondary 
currents in Fig. 8 as unity, then in Fig. 9, 
for the same current Ig, they become 2/7 
and +/6/7m respectively. The presence of 
secondary reactance may therefore make the 
resistive component of the cell loss and the 
Cu loss in the rectifier cubicle appear to be 
only about 60% of what they should be. 

The presence of secondary reactance will 
also have reduced the losses in the trans- 
former and interconnecting leads. For a 
6-phase double-Y or equivalent connection, 
these losses are allowed for by deducting the 
measured input watts at the same primary 
current with three of the a-c input terminals 
to the rectifier cubicle 120 degrees apart 
short-circuited (and also deducting the 
small loss in the interphase transformer, if 
any). In the case of Fig. 8, this is correct, 
but in Fig. 9 it is not correct because in this 
test, the secondary Cu loss is only two 
thirds of that which existed in the secondary 
when the rectifier was present in the short- 
circuit test. This will make the rectifier 
loss appear greater than it really was. 


SECONDARY PHASE 
CURRENT, RMS = 


PRIMARY PHASE 
CURRENT, RMS = 


Fig. 9. Rectifier current waveform in short- 
circuit test when secondary reactance is deter- 
mining feature 


iB!) 


PRIMARY We SECONDARIES 


Fig. 10. Arrangement of windings in test 
transformer 


Similar errors are present in testing single- 
phase rectifiers of the single-way type. 

Mr. Dortort has drawn attention to the 
first of these possible errors, and has sug- 
gested correcting for it according to the 
actual value of the form factor that existed 
in the short-circuit test on the rectifier, i.e., 
as in Fig. 2. Such a measurement of the 
form factor appears to present serious prac- 
tical difficulties. 

It is thus worth considering what can be 
done to reduce the effect of secondary react- 
ance. Neither extra primary reactance nor 
inductance in the d-c short-circuit path help, 
for in Fig. 9 the primary current is already a 
sine wave and the direct current is already 
perfectly smooth. Since it may well be im- 
practicable to reduce the secondary react- 
ance, the most effective means for reducing 
its effect, i.e., for reducing the secondary 
current spread, appears to be the introduc- 
tion of resistance in the secondary circuit, 
since this, unlike reactance, acts on the 
whole secondary current and not merely on 
its a-c component. It is necessary, how- 
ever, to introduce enough resistance to bring 
the secondary current fairly close to the 
half-sine-wave form of Fig. 8, since as long 
as the current is as in Fig. 9, the bigger the 
loss in the secondary circuit, the bigger will 
be the second error discussed above. 

It seems to me that for the purpose of 
standardization the best course in the case 
of single-way rectifiers will be to set a maxi- 
mum allowable limit to the departure from 
180 degrees of the secondary current spread. 
Further work appears to be necessary before 
we can establish what this permissible limit 
should be. This is very inconvenient, but I 
cannot at present see any alternative. 

7. In some low-voltage heavy-current 

_single-way rectifier equipments, a further 

substantial error may be present; namely,a 
change of the stray loss in the transformer 
according to how this is operating. An ex- 
ample will explain this. 

Among our early tests to measure recti- 
fier loss, we tried supplying a double-Y recti- 
fier cubicle from an old stepdown trans- 
former that was available in the works and 
had originally been built for a synchronous 
converter. This transformer was con- 
nected star/6-phase diametric which, under 
the conditions of rectifier operation, was 
equivalent to double-Y. At the voltage 
used, core loss was negligible. The following 
transformer input readings were obtained 
in the short-circuit tests: 


Transformer and rectifier (after deducting 

loss in d-c short-circuiting loop).......... 22kw 
Transformer only (three secondary termi- 

nals short-circuited), at same primary 

PUSEECHC 5 oy. sisca <3 oa cea pie em Wen Cant AeA 27 kw 


Examination showed that this transformer 
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PRIMARY 


Si 2 SECONDARIES 


Fig. 11. | Winding arrangement in some low- 
voltage single-way transformers 


had the center point of the secondary wind- 
ing brought out half way along the length 
of the coil, as shown in Fig. 10. The 
secondary reactance was very high, but 
the surprising result could not be accounted 
for by the effects discussed in section6. The 
explanation appears to be asfollows. When 
operating with the short-circuited rectifier, 
the currents (Fig. 9) in the two secondary 
part-windings in the same stack added up, 
so far as their a-c components were con- 
cerned, to a uniform complete sine wave of 
current occupying the whole length of the 
stack. But when only the transformer was 
short-circuited, across one Y, the secondary 
current occupied only half the length of the 
stack and thus half the length of the primary 
winding. In this latter condition the leak- 
age flux path was quite different, and in the 
transformer concerned accounted for a 
great increase of stray loss. 

This example cannot be altogether dis- 
missed as a freak case. Many low-voltage 
single-way transformers have the wind- 
ings arranged approximately as Fig. 11. 
This is obviously part way towards the ar- 
rangement of Fig. 10. Such transformers 
are often quite satisfactory for their purposes, 
and the working loss in them is correctly 
measured by the existing standard method. 
However, when used to measure rectifier 
loss they may give rise to serious error. 
This error is caused by the wide spread of 
the secondary current. 
same safeguard as was proposed in section 6 
will also guard against this. 


8. It is possible that still further points” 


need to be considered. However, sofar as 
we have yet gone, I now believe for the 
above reasons that for the purposes of stand- 
ardization, (a) the loss in the rectifier cubicle 
should preferably be taken as the input power 
to the cubicle measured in a single short-cir- 
cuit test, at normal factory temperature, at 


1.0874 for polyphase or at Ig for single-phase ~ 


rectifiers; (b) in the case of single-way rec- 
tifiers the secondary current spread in the 
test should not depart from 180 degrees by 
more than a definite amount (to be specified 
after further investigation); and (c) the rat- 
ing and losses of single-phase rectifier trans- 
formers should (unlike polyphase) corre- 
spond to sinusoidal currents giving the speci- 
fied direct current, not rectangular currents. 


REFERENCE 


1, RECOMMENDATIONS FOR Mercury Arc Con- 
VERTERS. Publication 84, International Electro- 
technical Commission, Geneva, Switzerland, 1957, 
clause 342, 


I. K. Dortort: It was obvious from the 
start that the proposed test methods have 
certain inherent weaknesses which require 
further study and experience. The three 


Consequently the 


_and the secondary of a single-way rectifi 


: 


gentlemen who have submitted discusalie 
have pinpointed some of these problems, 
particular they have questioned the re) 
tions under which tests are conducted, th 
marriage of rectifier and transformer, an 
validity of tests in which the condu 
period is greater than 180 degrees due 
inductance in the circuit elements. 

Mr. Borg has also brought up the inte 
ing case of extended conduction due 
inductance in the d-c circuit. In a 1 
paper,! I presented an analysis of a rectifie 
with all commutating reactance in th 
circuit but with the d-c short circuit in 
ing a large inductance. It was shown t 
in the ideal case conduction would last 2 
degrees. While the paper does not evalu 
the effect of finite d-c inductance, our ex 
ence leads us to believe that a separate t 
method to satisfy this condition will n 
necessaty except for extremely special c 

In the first place, it appears that the | 
inductance would have to be of the ordex 
of magnitude of the leakage inductance o# 
the transformer to produce conducting pex 
riods appreciably greater than 180 degre 
Second, losses in the circuit will furt 
decrease conduction beyond 180 degre 
In addition to the extended overlap, | 
d-c inductance will of course tend to m 
the cell currents flat-topped, but this e 
if appreciable, can be corrected in accordai 
with Fig. 2 of the paper. 7 

Actually, a high inductance in the sho 
circuiting bus is more likely to occur in 
double-Y rectifier rather than in a brid 
connected rectifier since the d-c termir 
are likely to be much further apart, an 
further contribution is made by the leak 
reactance of the interphase transforme 
Moreover, a double-Y connection is mora 
likely to be used for high-current low-voltage 
rectifiers in which the effect of d-c inductaneg} 
becomes more pronounced. 

More serious, and more prevalent, is 
presence of inductance in the unidirecti 
circuits, i.e., the circuit elements. T 
the condition referred toin the present pap 
Application of a correction based on fo 
factor from the curve of Fig. 2 seems to B 
the only practical solution and will give g 
results for the rectifier itself. However, 
extended conduction period reduces the 
losses in the transformer and makes pro 
calibration of the transformer more 4 
cult. Dr. Read’s suggestion of placing 
upper limit on the extended conducting fp 
riod is well taken. Just what this li 
should be is difficult to say at the mome 
Considerable study and experience will bh 
required to settle this question. Form-fa 
tor correction should be quite satisfact 
for the Cu loss in the unidirectional circui 


transformer. The Cu loss in the primar 
as well as in the secondary, of a double-wa 
transformer will be reduced due to cance 
tion of d-c components. Purely as a gue 
a conducting period of 190 degrees is 
gested as a starting point for discussion. 

Although no investigation has been made; 
it seems likely that by far the largest number 
of industrial rectifiers will have conducting} 
periods of less than 190 degrees. The form 
factor correction on the rectifier should bei 
adequate, and correction of the transform 
copper loss can be estimated fairly well from) 
the primary current readings taken during] 
the tests. | 
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Mr. Winograd questions the use of the 
scillograph for determining form factor, 
nd suggests determination of conducting 
ngle by means of a cathode-ray oscillograph. 
t appears to me that measurement of the 
onducting angle is not sufficient for deter- 
nination of the form factor when the angle 
s appreciably greater than 180 degrees. It 
yas not the intent to specify a magnetic 
scillograph for this work. A plot or pic- 
ure of the waveform displayed on a cathode- 
ay oscillograph should be sufficiently ac- 
urate. I certainly agree that further con- 
ideration is required but do not feel that 
york on the standards should be delayed. 
Phe effect of form factor has been of only 
econdary importance in the rectifiers tested 
o far. 

Regarding Mr. Winograd’s comment con- 
erning P,; and P, measurements at other 
han Igand KJ, it is obvious that loss meas- 
irements at three-fourths, one-half, and 
me-fourth load can be made easily and 
rithout too much loss of time by substitut- 
g fractional values for Jg but retaining the 
me ratio of direct current in the two 
easurements. The difficulty of maintain- 
ig exact values of direct current presents 
somewhat more complicated procedure. 
, instead of Jz and KJq, the test is made at 
a and bKIg, equation 6 of the paper 
anges to 


_a(K*—a)P,—bK"{K—b)P, 
abK(bK —a) 


a=6 so that the ratio of the two currents 
still K, the above expression is only 
ightly simplified to 
mi k*—a)P,—K\K—a)P; 
a?K(K—1) 


The specification of ambient-air or raw- 
oling-water temperature at the time of the 
Sst poses a difficult problem. In a recir- 
lating cooling system, whether air or 
ter, the temperature of the equipment 
n be controlled without too much difficulty 
rovided the raw-cooling-medium tempera- 
tes are low. If the raw-cooling-medium 
mperatures are higher than specified, 
ere would be no recourse but to wait for 
ore suitable conditions or to find some 
eans of reducing these temperatures. In 
Tect-cooling systems recirculation can 
ten be improvised. In any case, I fully 
tee that tight specifications could be 
erous and it might be wise at this time to 
ovide fairly wide latitude in the standards 
that tests can be made under prevailing 
mditions except in extreme cases. When 
tisfactory methods are worked out for 
aking temperature corrections, the per- 
issible band of ambient temperatures 
ight be broadened. 

While most standardization committees 
e agreed that the reverse current losses 
the cells are to be neglected in computing 
Jencies, it is obvious, as pointed out by 
r. Winograd, that voltage divider resistors, 
din some cases surge-suppression devices, 
t be taken into account. 
e rms value of the reverse voltage across 


At no load’ 


a circuit element in a 3-phase bridge-con- 
nected rectifier is 0.662 g. In the double-Y 
connection the rms value of reverse voltage 
under the same conditions is 1.33Eg,. From 
these values and the known value of resist- 
ance across a circuit element, the loss of a 
complete bridge rectifier is 8E£y.2/3R and for 
a double-Y rectifier it is 32Eg.2/3R. For 
the sake of simplicity it is suggested that 
nonlinear surge-suppression resistors be 
treated as linear elements since in most cases 
their nonlinearity becomes effective only at 
excess voltage. 

Dr. Read has brought up several points 
that were bound to cause concern and dis- 
cussion. While it seems to me that the 
combined transformer and rectifier loss test 
is the most advantageous, restriction to this 
method in the standards would be unwise and 
unfair. We should not impose higher testing 
costs by making it mandatory to marry the 
transformer and rectifier on the test floor. 
However, the increased cost is not chargeable 
to the time required to bring the transformer 
up to temperature. If the rectifier and 
transformer are not brought together, the 
transformer must still be subjected to a heat 
run unless it is a duplicate of a previous de- 
sign, previously tested. During the recti- 
fier test with a calibrated test transformer, 
we must still allow enough time to stabilize 
the temperature of the test transformer in 
order to make its calibration valid. 

With reference to single-phase rectifiers, 
there can be no serious objection to making 
K=1.0 for all but highly inductive loads. 
Even then, we would have no serious objec- 
tion to K=1.0 in order to put all tests on a 
comparable basis, if the standardization 
committees decide to do so. 

The suggestion to substitute a single loss 
measurement at 1.08 J, for the two measure- 
ments at Iz and KI, does not seem to be 
justified by facts. The statement is made 
that under normal standards of testing the 
proposed method of testing is subject to at 
least 3% error because the power loss is the 
difference of two large quantities, the P» and 
P, terms in equation 6. These terms have a 
ratio of approximately 2 to 1 and errors of 
measurement can therefore be assumed to 
be doubled in the result. Dr. Read has 
offered two examples to show that a single 
measurement at 1.08Jz would be more 
accurate. If we check P at the two theoret- 
ical limits of E,,=100% and E,,=0, we 
find that P=P, and 1.21P, respectively. 
Had we made a single measurement at 1.08 
Ig we would obtain an error of 8% in the 
first case and of 11% in the second. It is 
true that these represent extreme theoretical 
cases, but by now we have seen a variety of 
rectifiers in which both of these limits have 
been approached. 

After Dr. Read submitted his discussion 
of this paper, he carried out additional 
investigations and transmitted the results 
to me. Using the same copper-loss ratios 
as in the earlier investigation, he introduced 
into each case values of Z,, equal to 25 and 
75% of the total cell drop. In three of the 
four cases the computed error for a single 
measurement at K=1.08 ranged from 0.17 
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to 2.7%. When E;, represented the major 
portion of the total loss, the error was 
found to be +4.2%. 

It should be noted, however, that equa- 
tions 6 and 7 of the paper contain the 
differences of two quantities having a ratio 
of, roughly, 2to1. These quantities in turn 
contain P, and P; which differ by approxi- 
mately 15% and would normally be meas- 
ured with the same instruments and instru- 
ment transformers, all set for the same 
ranges. Therefore, their errors will not be 
accumulative. Instrument error will there- 
fore affect the end result to no greater extent 
than the instrument error of a single meas- 
urement at K=1.08. The single measure- 
ment would introduce an additional error 
above the instrument errors. 

Returning to the question of conducting 
periods greater than 180 degrees, it is not 
difficult to account for the negative rectifier 
loss reported by Dr. Read in a measurement 
made with a transformer having the winding” 
arrangement illustrated in his Fig. 10. We 
need not go into the nebulous realm of 
stray losses to account for this astonishing: 
result. Working only with J?R, a quick 
analysis shows that because of the greatly 
extended conduction period brought about 
by uncoupled secondaries the secondary cop- 
per losses obtained during the short-circuited 
rectifier test will be 39% less than when the 
transformer is tested alone in the approved 
manner. The primary losses are down 60%. 
It is unlikely that any rectifier transformers 
designed for industrial use, whether single- 
phase or 3-phase, would ever embody the 
construction shown in Fig. 10. The wind- 
ing arrangement shown in Fig. 11 is more 
reasonable and the effects of loose coupling 
would be reduced to about 1/4. Even this 
winding arrangement is contrary to usual 
practice. 

I wish to thank those who have con- 
tributed discussions on the proposed test 
method, and particularly Dr. Reed. It is 
necessary to bring to bear on this subject 
all the technical knowledge and the experi- 
ence that is being accumulated in various 
organizations. The need for standards is 
urgent and cannot wait until a perfect 
method is developed. Fortunately, the 
industrial applications in which standards 
are most needed and to which the bulk of 
power rectifier production is devoted, can be 
covered by standards which do not contain 
all the refinements discussed here. Ques- 
tions of combined tests and ambient tem- 
peratures are not peculiar to rectifiers and 
can be handled for the time being by reason- 
ably nonrestrictive test codes. Limits must 
be assigned to avoid application of these 
standards to rectifiers for which they are 
not suitable. As the art develops, the 
standards can be modified to cover rectifiers 
which must be excluded at this time. 
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Comparison of Steel and Aluminum 


Subway Cars 


L. W. BARDSLEY 


MEMBER AIEE 


ILD STEEL has been the usual 
basic structural material for sub- 
way-car strength members with aluminum 
generally utilized for decorative or low- 
stressed sections. 

The characteristics of aluminum alloys, 
i.e., high strength-to-weight ratio and 
excellent resistance to corrosion, appear 
to offer them a place in subway-car con- 
struction, where power, brake-shoe wear, 
and way-maintenance costs reflect car 
weight, and where corrosion resistance 
is a vital factor in body life. 

The Toronto Transit Commission 
(TTC), by purchasing at one time a num- 
ber of otherwise identical steel-shell and 
aluminum-shell subway cars, established 
a means whereby practical “in-service” 
‘comparisons of steel and aluminum in sub- 
way-car construction became possible. 

The TTC contracted for 104 steel sub- 
way cars for delivery in 1953-54, to 
operate its Yonge Street line, opened in 
1954. These cars were manufactured by 
Gloucester Railway Carriage and Wagon 
Company of England. During the course 
of car construction in 1952-53, the Com- 
mission was approached, through the 
manufacturer, by the Aluminum Develop- 
ment Association (ADA) of Great Britain 
to consider the construction of a number 
of these vehicles in aluminum rather than 
steel. 

As a result of discussion with the ADA 
and the car builder, the Commission re- 
vised its order to 100 cars of steel and four 
of aluminum with the general dimensions 
as noted in Fig. 1. By this means, since 
the cars were identical in all respects 
other than the aluminum application, 
in-service testing and experience would be 
able to determine if the operating gains 
through the use of aluminum would 
justify the increased initial capital costs. 

The advantages of aluminum compared 
with mild steel (SAE 1020) are: a much 
more favorable strength-to-weight ratio; 
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superior corrosion resistance; and its use 
in an unpainted condition. From these 
considerations, operating savings could be 
expected (1) in areas affected by weight, 
and (2) in maintenance costs through re- 
duced corrosion of body members. 

In any comparison of vehicles where 
performance is in question, one of the 
fundamental factors for consideration is 
the relation of schedule speed and power 
costs. Also, an examination of any 
operating line must consider passenger 
density, headway, cars per train, and 
schedule speed. Headway and train size 
are influenced by passenger density and 
associated with service facilities such as 
station lengths, station exits, entrances, 
andsoon. , 

Under a given set of conditions concern- 
ing passenger density, headway, and train 
size, a variety of scheduled speeds, num- 
ber of trains, and, consequently, cars re- 
quired for service can be ascertained. 
However, train speed variations reflect 
variations in power consumption. Foran 
economic analysis, therefore, the initial 
capital cost for cars required for a service 
should be analyzed in respect to annual 
operating costs of train crews, power con- 
sumption, brake-shoe life, and allied 
factors. 

In comparing vehicles such as steel and 
aluminum subway cars where a weight 
reduction is obtained by using aluminum, 
either of two choices may be made: 


1. Choice 1 would be to maintain the same 


10'V/ie” FREE 
93/6 COUPLED 


operating schedule with the aluminum or! 
lighter subway car because, for any givent 
weight reduction, the lighter car can main- 
tain the schedule of the heavier car at a: 
lower operating cost. / 


2. Choice 2 would be to utilize the weight! 
reduction to obtain a faster schedule speedi 
for the same power usage. The fastert 
schedule could result in a reduction in thes 
number of trains operating with a conse~. 
quent saving in train crews and initial capi-- 
tal outlay for rolling stock. | 


The TTC’s aluminum subway carss 
operate in mixed service with the steel 
cars and must, therefore, operate on thee 
same schedule. Thus, the comparison ofi 
performance must be as in choice 1, thatt 
is, determine the savings in operatingg 
costs while maintaining the same schedule:! 

With these factors in mind the Com- 
mission decided to establish comparatives 
tests on certain components and oper 
ating cost areas where definite savingss 
might be expected as follows: in fae~ 
tors influenced by weight reduction: 
(1) tractive power costs; (2) journal andd 
drive unit bearing life; (3) wheel ands 
brake-shoe life; (4) permanent way lifee 
and maintenance. In factors influencedd 
by corrosion resistance: (1) body main 
tenance costs. : 

By assigning values to the saving fae+ 
tors determined and amortizing the initial 
vehicle costs over a 30-year period, ag 
realistic picture of the financial relation 
ship of the steel and aluminum subwayy 
cars can be obtained. ¢ 

The relatively short in-service life ob 
the vehicles does not permit completed 
actual analysis. However, the calculatecc 
effects together with actual results whereq 
possible have been compiled. 


Structure and Weight 


The steel and aluminum subway cars} 
are basically similar in body appearance 


SPRING 
BUMPER 


Fig. 1. Toronto Transit Commission subway car plan and elevation 
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Fig. 2. Interior 
body shells 


A — Aluminum 
cars 
B—Steel cars 


Fig. 3. Under- 


frames 


A—Sills and 
body bolster of 
aluminum cars 
B—Sills, body 
bolsters, and 
heating ducts of 
steel cars 


no 


excluding the essential difference in 
‘material. 

_ The entire body shell of the aluminum 
car, including side frames and equip- 
ment mounts, is constructed of alumi- 
‘num alloys, with the one exception of 
‘the body bolster which is steel. This 
latter becomes necessary because of space 
testrictions and vehicle standardization 
coupled with the low elastic modulus of 
aluminum. 

Additionally, because of the small 
‘number of cars involved, the aluminum 
members were fabricated from standard 
‘shapes in many instances where, if 
“quantity permitted, a more advanta- 
-geous strength-to-weight ratio might have 
been obtained through the use of special 
extruded shapes possible with aluminum. 
- It would appear that, without limita- 
tions imposed by the steel prototype 
ood the necessity for standardization for 
“equipment mounting, a lighter aluminum 
‘ear could be produced within the dimen- 
‘sional framework imposed by the subway 
‘Structure. 

_ The aluminum alloys used in the alu- 
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(A) 


minuin car design varied with requirement. 
For example, the main structural and 
other highly stressed members were con- 
structed from an aluminum magnesium 
silicon alloy containing a small percentage 
of copper and chromium. The exterior 
panels, which were to be unpainted, used 
an aluminum 31'/.% magnesium alloy, 
cold-rolled for strength. 

The steel cars were constructed basi- 
cally from mild steel, although aluminum 
was used in some cars for roof sheets and 
doors and in all cars for trim and other 
decorative purposes. The steel cars were 
also painted in the standard Commission 
color scheme. 

Figs. 2 and 3, depicting underframe and 


(B) 


side-frame construction, illustrate varia- 
tions in design between the aluminum 
and steel cars. 

Car weights of the 100 steel cars varied 
because of an increasing use of aluminum 
in nonstrength locations as car construc- 
tion progressed. Table I lists car weights. 
for all car types used in the Toronto sub- 
way. 


Power Consumption 


Considering the probable effect of 
weight reduction on power consumption, 
the ratio of power to weight saving is 
an important factor. The fundamental 
equation of tractive effort: weight < accel- 


Table |. Toronto Transit Subway Cars 
Number in Basic Structural Weight, 

Car Nos. Group Type Material Lbs 
5000-5029........ BO cecis aeetacevanioane driving, motored|..< 2. .2 js 605 ce Stealisc# Aydt ccdoalene anaes 85,525 
5030-5099........ WOW ae Satter. cvs tte: (chaighalss sreVolKoaeobn, Son ping Oooo noS steel... ‘ 83,670 
5LOO-5LOD:. see. a0 (CYS geet te dnvingsimotored.... saen. see ole ENikorieehb oleae gsang ganas 73,250: 
PLAS epL LE, cehauetace a G6, owhalee teat Grinxime motored. a-tisiersteraici<s oases steel . 82,726 

(dynamic braked) 
5200-5227........ A: Se Pa erty EC nondriving motored...........-.. SEEEL 6; ssoyehepresicys. wadbels eu 76,720 
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Fig. 4. 


(A) 


SECONDS 


(8) 


Speed and current time curves; 10,000-lb load; down 0.56% grade; run 


length, 2,145 feet 


A—Type G2 aluminum cars; 73,250-lb tare 
B—Type G1 steel cars; 83,670-Ib tare. 


eration/g, would suggest that power con- 
sumption and weight should vary directly. 
However, it is equally apparent that other 
factors, such as losses due to rotating 
equipment within a vehicle; train resist- 
ance as a function of vehicle shape, 
weight, and speed; and friction at bear- 
ings, contribute to modifying the simple 
relationship of force and acceleration. 

Two methods have been employed by 
the Commission in attempting to obtain 
power consumption comparison of the 
steel and aluminum subway cars: (1) 
calculation and (2) direct measurement 
(in service through the use of installed 
meters). 


Calculation 


For comparative calculation of the 
effect of weight reduction on power con- 
sumption certain conditions were estab- 
lished. 


1. Length of run assumed to be 2,145 feet 
representing the average between-station 
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subway distance on the Yonge Street line. 


2. Anaverage run on tangent track of 1.2% 
up and 0.56% down. 


3. Acceleration rates approximately the 
same on both car types. 


4. Tare weight of 7steel car, 83,670 lbs 
(pounds). 


- 5. Tare weight of aluminum car, 73,250 


Ibs. 

6. Assumed passenger weight per we 
10,000 Ibs. ; 
7. Same schedule speed. . q 


8. Same motors and field shunting. 


The cars under consideration were 
identical in respect to body profiles 
mounted underbody equipment, trucks 
and motors. For comparison, scheduled 
speed or run time was held approximately ; 
constant for each car. The motors 
were identical in characteristics and an 
shunting; the only effective variable ot 
than weight which might affect energ g 
consumption was accelerating ampe 
and coasting. 

Within the framework of identical time q 
and distance, a number of variations - ii 
accelerating amperes and coasting can be 
obtained. For practical comparison hows 
ever, the most realistic assumption woulc 
be to maintain the same acceleration as 
mutch as possible in both cases. ; 

On the basis of the conditions outlin dai 
performance speed time and curreng 
time curves for the aluminum and stee 
cars were developed as noted in Figs. 4(A 
and (B). From these curves the data 
Table II are derived. 


Direct Measurement 


i Mis i tle 
— + 


With a view to subsequent power-con 1 
sumption tests, the Commission during 
car construction arranged for the inclusio# 
of eight resilient mounted watt-hov 
meters in four steel cars, nos. 5096-5095 
and four aluminum cars, nos. 5100-51 r 
for measuring traction power cone 
tion. 

The cars entered service in March 1954 
and monthly readings of the meters wert 
taken. The average power consumptio’ Of 
obtained through these meters betweed 
March 1954 and January 1956 are give 
in Table ITI. | 


5 
Table Il q 

Weight Per % Weight Power Consumption, Powe: 
Car Type Car* Reduction Kw-Hrt/Car Mile Roduet op 

; 3 

Steel cars 5096-5099...... 02670 Sy Stee ete: 2 Te eae agrees 4 
Aluminum cars i 
DIOO STORM. say as, aca ey als 832508 ce, sole seas Bi Oieveicccvaiiese. «\ Soci aon jeer CL ee re fee Portas 8.3 ; 

*Weight is tare plus 10,000-1b load. 
+ Kilowatt-hours. 
Table III 
i % Powe “ 
: % Weight Power Consumption, Consumpti 

Car Type Weight* Change Kw-Hr/Car Mile Chang o 
Steel cars 5096-5099........+45. $3 '.67.0:. Ret th o/s re at) ‘- 
Aluminum cars 5100-5103....... MS. 2O0 en kee anit cties iene ae meneame rt AQs daha nee eee 8.3 


*In-service weight not known. 


Bardsley—Comparison of Steel and Aluminum Subway Cars 


Table IV. Load Conditions Journal Bearing Life 


— 


® Operating Time 


Steel Car Weight, Lbs 


Aluminum Car Weight, Lbs 


Total Weight* 


Total Weight* 


N) Tare* Passenger Tare Passenger (F) 
51 21,000... 105 ,010 21,000 94,2 

ODS So ada e ; (OWS ade ,250 
MEU sae SaOl Oke vas see TL 2OOR cimrct res Qos 2LOM Ry aan TORZOO Gc s.0 3 E2005. 22. 84,450 
.96 Mi OOO Ska wscae 91,010 GO00 fort r. 80 , 250 
.61 A LOO ni: cusyecexers 88,210 BEZOO™ pistes 77,450 


Steel car weight is average for 100 steel cars. 


The comparison between the calculated 
id direct-measurement power-consump- 
on figures indicates identical answers. 
his accuracy is, we feel, coincidental. 
irect measurement of power includes 
atiables of passenger loading, train 
akeup, and actual coasting. To deter- 
ine these figures theoretically, partic- 
atly based on an average run calcula- 
on, the accuracy obtained is extremely 
iestionable. Coincidence in this case 
probably the answer. 

Coasting is another extremely impor- 
nt variable. It may be seen from the 
itves in Fig. 4 that for a very small 
crease to run time, an appreciable 
nount of coasting can be inserted 
hich results in a major reduction in 
ower consumption. 

For the purposes of this study, how- 
yer, the actual power saving of 0.40 kw- 
‘/car mile will be used. 


earing Life 


Although both the motor-drive bearings 
id the load or journal bearings are 
tected by weight reduction, the possible 
ing in life of drive bearings has been 
mitted as a negligible factor because of 
ie bearing cost. 

The analysis of the weight effect on 
earings is, therefore, confined to the 
uurnal bearings, of which there are 
ght pairs or 16 bearings per car. 

The general antifriction bearing for- 
ula is: 


C 
=(£).0/8 


here 


=minimum bearing life, millions of revo- 
lutions 

=constant for bearing assembly 

=equivalent dynamic load 


‘The minimum bearing life, L, for the 
oronto aluminum subway cars has been 
culated as 1.61 times that for the steel 
‘The load conditions and calculation 
lethod to determine this factor are il- 
strated in Table IV and Appendix I 


‘The minimum bearing life L, since it 
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represents revolutions, in the case of road 
wheel can be transposed into miles by 
calculating the number of revolutions per 
mile by the 30-inch wheels used on the 
subway cars. This approximates 150,000 
miles for the steel cars and 241,000 miles 
for the aluminum cars. 

The bearing industry subscribes to the 
theory that “average” bearing life, on 
which costs are based, can be assumed to 
be five times the minimum or what is com- 
monly referred to as B-10 life. On this 
basis the average expected bearing life on 
the Toronto subway cars under discussion 
should be, for steel cars: 750,000 miles, 
and for aluminum cars: 1,207,500 miles. 

If it is assumed that car miles per year = 
55,000 and the number of bearings per 
car=16, then annual bearing saving per 
year per car=0.47R, where R is the bear- 
ing cost. 

Factual data are not yet available be- 
cause of the short service life, so the 
theoretical figure only is used in subse- 
quent evaluation. 


Brake-Shoe Life 


Within the normal operating speed limit 
of 50 mph (miles per hour) on the Toronto 
subway, the maximum kinetic energy for 
an 83,670-lb steel subway car with its 
maximum loading of 31,000 Ibs (222 
passengers) is determined from the kinetic 
energy formula 


Work /Sec. in IO ft Ibs. 


Fig. 5. Relation of 
Diamond S_ brake- 


KE=~ y 
2g 
where 
W =weight, lbs 
V=velocity, ft/sec (feet per second) 
This results in a KE value of 960,000 ft- 
lbs per sec, or 60,000 ft-lbs per sec per 
shoe. 

Within this range of work per shoe, 
experiments conducted at the University 
of Illinois Experimental Station have 
indicated that the relationship of brake- 
shoe wear, in terms of lbs per 100 million 
ft-lbs of work, approximates a linear rela- 
tion to the average rate at which braking 
is done; see Fig. 5 and Appendix II. 

From this relationship it can be derived, 
as noted in Appendix II, that brake-shoe 
wear varies directly as the square of the 
car weight for a vehicle equipped with 
cast-iron brake shoes. 

It would be expected therefore, that in 
the aluminum cars cast-iron brake-shoe 
life should be increased by the square of 
the ratio of car weight or: 


Brake-shoe life, aluminum car 


2 
= MEW) ste easaae Xbrake-shoe life, 
W aluminum car 


steel car 


Assuming a 10,000-Ib passenger load 
in each case: 


Brake-shoe life, aluminum cars=1.24x 
brake-shoe life, steel car 


Actual tests during 1955 and 1956 
showed the following: average brake-shoe 
life, aluminum cars=7,550 miles per shoe; 
average brake-shoe life, steel cars=6,220 
miles per shoe; therefore actual brake- 
shoe life of aluminum cars = 1.214 brake- 
shoe life of steel cars. 

The discrepancy between theoretical 
and actual results may be due to the 
assumption that the brake-shoe wear 
varies linearly with the rate of doing work 


shoe wear to work Pounds of Shoe Wear per IO ft Ibs. work 


done in braking 
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and validity of the assumption of pas- 
senger loading. 

Within the test period, because of a 
change in service train size requirements, 
the aluminum cars were mixed indis- 
criminately with steel which resulted in 
higher work requirements for the alumi- 
num car brake shoes and hence decreased 
life. 

For analysis, however, the agreement 
appears satisfactory although the in- 
service performance record will be used to 
compare savings which come to 21.4%. 
Since the TTC average brake-shoe con- 
sumption per car per year for 104 cars is 
118, the annual saving per car per year 
for the aluminum cars may be stated as: 
0.214X118XA=25A where A repre- 
sents the cost per shoe. 


Wheel and Track Wear 


The effect of weight reduction on 
wheel and track wear is difficult to deter- 
mine. Tangential shear components in 
propulsion, dynamic loading in motion, 
and tangential effects on wheel and track 
are all effective factors in wheel and track 
wear. 

In considering the variation in life of 
wheels and track work due to weight re- 
duction, it is probably safe to assume that 
since acceleration and braking rates are 
maintained constant on the steel and 
aluminum cars the main governing factor 
would be dynamic effects. For the wheel 
these effects are generally estimated 
at 100% of static weight and therefore 
would be expected to vary directly with 
weight change. 

It has been determined from service 
operation that 1/16-inch diametral wheel 
wear represents 5,000 car miles for the 
aluminum cars and 4,400 car miles for the 
steel cars. This means the weight reduc- 
tion due to the use of aluminum has re- 
sulted in an increase of wheel life of 
13.7%. 

Present experience indicates wheel 
replacement for a steel car to be 21/4 years 


or a per-car annual wheel replacement of 
32/9 wheels per year. Assuming Z =cost 
of the wheel, the annual replacement cost 
of steel car wheels per car =32/9 Z. 

For the aluminum cars, the annual sav- 
ing per car is 13.7% of this cost or (0.137 X 
32/9 Z) =0 49 Z. 

The variation in wear rate on wheels 
may reflect similar variations in the rate 
of rail wear under similar track mounting 
and operating conditions. With a mixed 
operation over common rail, actual 
measurements are not possible. For 
this study no theoretical evaluation for 
track wear has been attempted and thus 
no value has been used in considering 
relative economies between aluminum 
and steel cars. 


Body Maintenance 


Resistance to corrosion of body com- 
ponents is a major factor in body main- 
tenance. The corrosion rate in outer 
sheets, body posts, sills, purlines, carlines, 


- and other members is the prime factor in 


establishing economical vehicle life. 

Protective coatings are applied to mild 
steel surfaces during car construction 
to minimize corrosion. Some of these 
coatings, such as exterior paint, require 
frequent renewal and examination to 
afford adequate protection. 

With aluminum a choice of corrosion 
resistive alloys such as magnesium and 
magnesium silicons with high corrosion 
resistance eliminates the necessity for 
protective coatings. 

It is possible to utilize the aluminum 
alloy in an unpainted condition. [If color 
is desired, suitable paint finishes can be 
obtained. The natural aluminum oxide 
film over which paint is applied is inert 
and this tends to maintain adherence of 
the paint to the parent metal better than 
on mild steel sheeting. 

The short in-service life of the subway 
cars has not allowed sufficient time to 
develop a complete picture of the relative 
body maintenance costs of steel and 


Table V. Car Operation Savings Accrued Through Use of Aluminum 


Savings 
Annua) Saving, 
By By Dollars per 
Item Calculation Measurement Factor Value Car 

POWEr Ss eure rarsismiseleus oco.e whe cate shane ese 0.40 kw-hr/car mile.......$44.70/kw-year.. .. 278.00 
Bearings i ae mrtotenere cise: Quad ARGS atcteys ctavaraehernetee Mem tees Cheers ete R = 59.00... .. 28.00 
Brake SHOES sy ciris wiley try ayers aloes seats Wale 2OTAG A=4.87..... ..109.00 
Wheels Sin jatcariccrcntessccts citas Cromreita st 0.49 Z. Z=102.00 -. 50.00 

pubtotal'of:savings dileito weight reduction.......— cso ee nee .. 465.00 
Body maintenance...... BBB 5 OO cok ccReu cnsvciinrate vlerstuvcvees te, Sle nactncte MRT ER EL en Fe 66.00* 
Total allsavitigs |... sesdertiouteucns £(s) Hjs"R: gaia! (a. a iio. mjibisa: age vo ca: ice 9 las @ehte, a] aieaval cartels a oe ee aaa eraiort a. ae 531.00 


* The present Commission contract for power calls for a demand payment plus an energy charge over a 


30% load factor. 
watts used at the car. 


In brief, the contract is evaluated to represent an annual charge of $44.70 for d-c kilo- 
Thus from kw-hrs per car mile and using the average schedule speed of the vehicle a 


kw demand figure is obtained. This demand figure multiplied by the annual demand charge represents 


annual power costs. 


Average for a system in rush hour including turn-around time 15.56 mph. 
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aluminum subway cars. However, ox 
the painted steel cars, corrosion at river 
heads in certain areas has necessitate« 
paint renewel at these sections. 4 
It appears that the steel subway can 
will require repainting after 6 yeateiny 
service. The repainting. cost would t 
approximately $600.00 per car representa 
ing a charge of approximately $100.00 pei 
car per year. t 
Future experience may show that th) 
figure $100.00 for the reduction in bods 
costs by the use of aluminum is either to, 
low or too high, but at present it appeari 
reasonable for use in this analysis. ; ] 


Conclusions 


The savings summary, Table V, i 
cates that the operating cost reductiog 
obtained through the use of aluminum i 
car construction is $565.00 per car p 
year. Not all possible savings have | 
included, so that the figure represen s 
minimum. 7 

From the savings picture two pertiner 
facts can be drawn: 


1. The weight saved in car construction 
that is by reducing weight of a car bi 
10,414 Ibs, represents an annual saving ¢ 
$465.00 which is equivalent to $0.05 a 
of weight saved. } 


2. If cars are amortized over 30 years ai 
money obtained at 5% interest, the annual 
savings of $565.00 represents the ability t 
pay approximately $8,700 more on original 
capital investment. 


It is of interest to note that on 
aluminum and steel cars in question, 


aluminum cars was $3,600 higher etal id 
steel. Subsequent tenders for steel anij 
aluminum cars increased this price diffe 
ence to $7,700. However, the perfornt 
ance analysis indicates the aluminum 
operating costs would justify these diffez 
ences. 
It should be realized that the additigl 
capital investment which could be justif 
fied by the annual saving is obviously } 
direct function of current interest rat 
and would, of course, have to be co: 
sidered at the time of vehicle purchase, 
It should be noted that the weight amt 
body maintenance savings attributed t] 
the use of aluminum cars are not specifi} 
cally confined to cars constructed with th 
metal. An analysis could be made Act 
stainless steels, wherein the inherent corre| 
sion resistance and extra strength wo 1h 
bear similar analysis. 
The comparisons outlined for # 
Toronto Transit Commission’s alumin' i 
and steel subway cars indicate the veri 


Whether weight saving is accomplished 
either structure redesign or the use 
alloy, steels, aluminum, or similar 
etals, the obvious fact is that every 
und of weight in a vehicle should be 
refully scrutinized. 
On the cars described, certain savings 
uld be obtained through the use of 
mamic braking, by a reduction in 
ake-shoe and wheel costs, and weight 
duction in car heating equipment. 
1ese factors would, of course, have to be 
nsidered in any car purchase. 
Present experience would indicate the 
e of aluminum in subway-car construc- 
m does result in an appreciable operat- 
g cost reduction and should be given 
rious consideration in the future. 
It is obvious that the acquisition of 
aminum subway cars has provided a 
uy to obtain factual data on the merits 
aluminum in car construction which will 
rve as a useful guide for future vehicle 
chases. 


Appendix | 


The following calculations are based on 
e Roller Bearing Engineers Committee 
port.! See Table IV. 


q \ F)N 
ae 100 


1ere 


=car load 


Ly 1961 


F'n = effective load 
N=operating time load is applied 


From this: 


Fy, steel=92,410 Ibs per car or 11,551 Ibs 
per wheel 

Fy, aluminum =81,879 lbs per car on 10,235 
lbs per wheel 


For individual bearing analysis L =(C/P) 
10/3 where C=81,500 lbs for TTC bearing 
assembly. 


P=XVF,—YF, 
where 


X =radial factor 
Y=thrust factor 

V =rotation factor 
F, =radial load 

Fy, =thrust load 


For Toronto subway cars P=approxi- 
mately 1.89 effective load per wheel, there- 
fore, for steel cars: 


P=1.89X 11,551 =21,840 
and for aluminum cars: 


P =1.89 X 10,235 = 18,940 
L steel=(C/P) 10/3 =(81,500/21,840) 
10/3 =80.48 
L aluminum =(81,500/18,940) 10/3 =129.1 


Therefore: 


Life rate of aluminum =Z aluminum/Z steel 
=129.1/80.48=1.61 times that of steel-car 
bearings. 


Appendix Il 


The following data are taken from the 


University of Illinois Station data and 
graph, Fig. 5. 


Pounds of brake-shoe wear/100 million ft- 
lbs work=K ft-lbs/sec 


If Ws;=work/wheel/stop in ft-lbs and 
X =Ilbs of brake-shoe wear/stop, then: 


x W. 
—— ah 
108 ft-lbs t 


and 
xX 
108 ft-lbs 
pounds shoe wear/stop X 108 ft-lbs 
= = 
Thus, 


tX108 ft-lb 
However, (W,) =kinetic energy =1/2( W/g)- 
V? where W=weight per wheel. 


Thus, it is noted shoe wear varies directly 
as the square of the vehicle weight. 


lbs shoe wear/stop =K 
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Predictive-Control System Application 


H. CHESTNUT 


MEMBER AIEE 


REDICTIVE control is used in this 

report to describe a form of auto- 
matic control in which the manipulated 
variable operating the controlled system 
is actuated by an estimate of the error 
which will exist at some future time. 
Repeated estimations of the future error 
are obtained by predicting ahead, on a 
fast-time base, both the reference and the 
controlled variable as well as some of their 
lower-order derivatives. Data smoothing 
and prediction of the reference may be 
required, and use of a model of the con- 
trolled system is needed. By virtue of 
being able to use fast switching at high 
energy levels and continually predicting 
a new value of future error, it is possible 
to obtain some of the desirable features 
of the bang-bang “‘optimum’”’ control as 
well as the dynamic programming con- 
cepts. 

One form of the idea of predictive con- 
trol has been described by Noton and 
Coales,! and work has been done on this 
method of control by J. C. Lozier of the 
Bell Telephone Laboratories. In a sense, 
it is similar to the application of some of 
the prediction concepts used in fire-con- 
trol and missile-guidance schemes. 

Use of predictive control appears highly 
attractive with the availability of small, 
high-speed logic devices and power 
actuating means. With predictive con- 
trol, any disturbances or reference input 
value, whether large or small, can bring 
to bear the full power capability of the 
power element. As such, with deter- 
ministic inputs, faster speeds of response 
for a given input and power source can be 
obtained than for a linear system. 

By virtue of the repetitive prediction 
calculations, including in the model the 
major nonlinearities or time variations 
of the process if necessary, higher gain 
with adequate stability can be obtained 
than would be possible with linear 
means. The equivalent speed of response 
of the system is markedly dependent on 
the speed of the repetitive predictions as 
well as by the ability of the system to 
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sense its own performance including its 
derivatives. 

The use of prediction for the reference 
signal as well as for the model of the 
process permits the optimum bang-bang 
concept to be applied to systems in which 
there is no constant reference input. As 
a result, the control can achieve a greater 
correction in a shorter time. Further, 
evidence points to the fact that an exact 
model of the process is not necessary and, 
as illustrated in this paper, a simpler 
second- or third-order model using 
equivalent time constants may be used 
with considerable effectiveness. 

Since the predictive control is a non- 
linear one, the amount of overshoot to 
step changes in reference input or dis- 
turbance is not proportional to the 
initialerror. The overshoot increases less 
rapidly with increasing input signals than 
for the linear control and, as such, more 
satisfactory performance is obtained 
where rapid synchronizing is required. 


Applications 


The characteristics of predictive con- 
trol make it highly attractive to a large 
number of computer-control applications. 
Although eminently suited for on-line 
digital computers, the method can also be 
used with analog computing means. 
Each computing means has certain fea- 
tures that make it attractive for partic- 
ular cases. Typical application areas 
where predictive control appears to be 
well suited are discussed in the following. 

Space navigation and rendezvous mis- 
sions can employ a form of predictive con- 
trol either to minimize energy expenditure 
or to minimize the time required to accom- 
plish a correction. Not only may the 
guidance be made to achieve a certain 
position and velocity but also it will seek 
out a means whereby these may be ac- 
complished within the heating or other 


REFERENCE r(t) AND REFERENCE RATE r(t) 


PREDICTED REFERENCE r(t¢) AND LOGIC 
PREDICTED REFERENCE RATE r(tr) 


SMOOTH 
AND 


FAST-TIME 
PREDICTOR 


Fig. 1. 


—=—a=d 


ACTUATING | 

4 orth CONTROLLED OUTPUT c(t)AND | 

CONTROL FS CONTROLLED | | CONTROLLED OUTPUT RATE é(t) 
SYSTEM ; 


TRIAL ACTUATING 
SIGNAL M(t¢ ) 


structural restraints, if this is possibl) 
with the power or other conditions ave 
able. 

Landing of the vehicles to presctilfes 
flight paths, within structural or aerct 
dynamic limitations, can be done eff 
tively with predictive control. Sines 
errors in the initial conditions or dici 
turbances during landing can b 
eliminated as rapidly as can be accony 
plished with the full torque or force cap 
bility of the controlled system, this con 
trol means has performance supefici 
to linear systems. An example of th 
application of predictive control to ft 
aircraft landing problem is discussed i 
this paper. | 

The ability of predictive systems t} 
operate with inexact models makes -| 
highly desirable for use in dynamic chen 
ical-control processes. The work .| 
Eckman and Lefkowitz? has shown 7 
results to be relatively insensitive” 4 
model inaccuracies. Thus the contra 
chemical processes in which the characte 
istics change with time or other parant 
eters in ways difficult to measure, 
be made to advantage with predict 
control. If necessary, it should be poss 
ble to combine predictive control witi} 
slower speed changes in the model, < 
through such self-checking procedure] 
as those of Eckman and Lefkowitz? __ 

From studies of the sort peal 


the material that follows, it appears t 

predictive control can be used to aq 
vantage in general, on-line, computer con 
trol.4 ; j 


Description of Predictive-Control 

Systems | 
St 

Predictive-control systems are simild 
to nonlinear optimum switching system 
in that they control the polarity of ef 
actuating signal such that the out 
of a plant, or other controlled systen 
synchronizes with a given reference inpp 
in minimum time. The means by wha 
this controlled switching is perforr 
is by the use of a fast-time model al 
control logic as shown in Fig. 1. | 
| 


actuating signalcan be a constant, a ram) 
a set of pulses, or whatever predeted 


| 


- 


p S2 6 PRESENT STATE CONDITIONS 


PREDICTED OUTPUT c(t¢) AND 
PREDICTED OUTPUT RATE ¢ (ts) 


FAST-TIME 
MODEL 


Block diagram of predictive control system 
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INITIAL AIRCRAFT ELEVATION 
eye RATE OF DESCENT 


a7 CONTROLLED AIRCRAFT 
POSITION 


t 


LEVATION 


REFERENCE 


HORIZONTAL DISTANCE —» 


lined characteristic is obtainable. The 
ontrol logic determines the polarity and 
uration of the actuating signal that is 
pplied to the controlled system and also 
he switching sequence by observing the 
ehavior of the predicted error and error 
ate in the event of reversing the actuating 
ignal polarity at the present time. The 
redicted error and error rate in the 
uture are obtained by smoothing and 
rediction of the reference input and by 


he fast-time simulation of the controlled | 


ystem, 

Fig. 1 shows that the actuating signal, 
Z(t), is applied to the controlled system 
ia switch S, and the trial actuating signal 
A(t) to the fast-time model via switch 
2. Switch S» also applies to the present- 
tate conditions existing in the controlled 
ystem at the present time on a variable- 
requency repetitive basis. During syn- 
hronization, switch S, operates much 
nore frequently than S,;. During close 
ollowing, S, and S: operate at approxi- 
nately the same rate. 

The control logic maintains a given 
etuating signal polarity until such time 
s the predicted error and error rate reach 
ero together. At this time, switch S; 
perates to reverse the actuating signal 
larity. 

Consider the case shown in Fig. 2 which 
lescribes the controlled landing of an 
ircraft and employs a variable-frequency 
epetitive prediction basis. The refer- 
nee shows a certain desirable path to 
ollow in order to touch down on the run- 
vay at the right point with an acceptably 
mall downward velocity. At point a the 
ircraft is too high and descending too 
lowly. The control logic immediately 
alls for a negative polarity actuating 
ignal which causes the aircraft to follow 
he trajectory as shown. The dotted 


ULY 1961 


Fig. 2. Predictive 

control applied to 

automatic landing of 
aircraft 


qs TOUCHDOWN 
YY 


RUNW 


curve proceeding directly downward 
from point j indicates the trajectory in the 
event that the polarity is never reversed. 
At point 6, where the error and error rate 
become of the opposite polarity, the fast- 
time model determines the trajectory in 
the event of polarity reversal at that time. 
The dotted curve starting at b indicates 
this predicted trajectory which is ex- 
tended until the predicted error rate is 
equal to zero. Since this trajectory does 
not intersect the reference where the pre- 
dicted error rate equals zero, this is not 
the proper time to switch. A second pre- 
diction is made at point c with the same 
conclusion. The prediction is repeated 
until finally, at point 7, the predicted 
error aad error rate reach zero (or some 
acceptably small value) together. The 
control logic switches the actuating signal 
polarity at this point. The trajectory 
therefore follows the solid curve to point 
k, Beyond k a small hunting exists 
about the reference. 


The error phase-plane portrait for the 
situation shown in Fig. 2 is shown in Fig. 
3. At starting point a, the error is 
negative because the output is greater 
than the reference and the error rate is 
negative because the negative output rate 
is less than the negative reference rate. 
The actuating signal is held negative 
through point j. At point 0 the first 
prediction is made to observe the value 
of e(ty) when e(t;) reaches zero in the 
event of a polarity reversal. Successive 
predictions are made until finally, at 
point k, the predicted error changed its 
sign from minus to plus when e(t,) 
reached zero. This information is used 
by the control logic to switch the ac- 
tuating signal polarity at point 7. 

The actual error phase-plane trajectory 
passes from point a through 8, c, d, 7, and 
k. Along this curve time takes on real- 
time dimensions. The predicted trajec- 
tories shown dotted in Fig. 3 are obtained 
from the fast-time model and along these 
trajectories time takes on computer-time 
dimensions. 


It will be noted that the predicted tra- 
jectory starting at b is much shorter than 
the one starting at d. This occurs be- 
cause the next prediction is started when 
the previous predicted error rate reached 
zero. Thisresults in a variable-frequency 
prediction rate, which means that, for 
large error and error rates, the prediction 
rate is slowest. For small error and error 
rates, the prediction rate is fastest. This 
situation is desirable because it results in 
tightest control when it is needed most. 


Beyond point & there exists a limit 
cycle about the origin. The magnitude 
and frequency of this limit cycle depend 
in large part upon the prediction repeti- 
tion rate. Because this rate is relatively 
high, the magnitude of oscillation can be 
small and the frequency high depending 


é (t) (ERROR RATE) 


— ACTUAL TRAJECTORY 
——-— PREDICTED TRAJECTORY 


e (t) (ERROR) 


Fig. 3. Phase-plane portrait 
for aircraft landing of Fig. 2 
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COMPUTE e(t) AND é(t) FROM 
r(t), F(t), c(t) AND C(t) 


@ DETERMINE SIGNS OF e(t) AND é(t) 


® COMPARE SIGNS OF e(t) AND é(t) 


DIFF. SAME 


INTRODUCE PRESENT STATE 
CONDITIONS AND REFERENCE! 
INPUTS INTO FAST-TIME 
PREDICTORS 


@ 


PREDICT FUTURE ERROR 
@] — TRAvEcToRY usING 
SGN M(tg)= SGN e(t) 


@| STP Preoicrine | 
WHEN @(t¢) = 0 


HOLD SGN M(t) 
= SGN eltfllartsy = 0 


MAKE SGN M(t) = SGN e(t) 


Fig. 4. Logic flow chart 


on the nature of the controlled system and 
the speed of the equipment to perform 
the switching functions. 


Sophistication of the control logic can — 


provide means for reducing the limit- 
cycle magnitude or for removing it 
entirely. For example, in one of the 
cases described later, use is made of the 
principle of reducing the magnitude of the 
limit cycle by reducing the amplitude of 
the actuating signal when the present 
error and error rate are less than certain 
threshold values. A number of other 
more elaborate methods might also 
be employed to reduce the limit-cycle 
amplitude. 

The preceding description employed 
the variable-frequency prediction method 
which is associated with analog computing 
equipment. In similar fashion, the fast- 
time prediction and comparison, if done 
on a digital computer, will probably be 
done on a time-sampled basis and the 
comparison made at discrete instances 
of future time. In this event, the time 
when this prediction is started is more 
rigidly prescribed and a somewhat differ- 
ent format for performing the control 


logic must be employed. rit) elt) SGN e(t) |} 
S ae ce | 
Control Logic : as. On aa ® SYSTEM |é(t) 
r(t) + e(t) SGN é(t) | WHEN SGN e(t) Z(t) gz(t) 
eds sat 03 rece #sonatt) gt | | ; | 
The heart of the predictive-control sys- =~ Giri ta =f PRESENT 
tem is the logic section. In many ways \ 2 CONDITIONS 
the control logic in a predictive system aie ane 4 adalat NO vee | 
a P SMOOTH HOLD \ = 
and the controller in a linear feedback AND \ k 
5 a cenit : FAST-TIME (3) \ - 
ystem perform similar functions except rica? \ aS c(tf) 
that each accomplishes its ultimate goals —t Fagen elh abel | 
based on different judgment rules. el | 
In its simplest form the control logic ht 
is just a means of mechanizing the switch- 
ing criteria that were introduced in the Fig. 5. Analog mechanization of control logic for predictive control 
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previous section. This procedure is 
summarized in a logic flow chart which is 
shown in Fig. 4. From this chart it can 
be seen that the first step in applying the 
logic is to compute the present values of 
e(t) and e(#) from the present values of 
r(t), r(é), c(t), and c(é). After determin- 
ing the signs of e(¢) and e(¢) they are com- 
pared for the same or different sign. 

If the signs of e(¢) and e(t) are the same, 
the controlled variable is moving away 
from synchronization; thus the sign of 
M(t) is made the same as the sign of e(t). 
This will insure that the controlled system 
is driving toward synchronization as 
quickly as possible. The choice of polar- 
ity for M(f) is apparent so that no pre- 
diction is necessary. It is of interest to 
note that a linear control or human opera- 
tor would also react in the same manner. 
An example of this condition would be 
segment a through 6 in Figs. 2 and 3. 

When the signs of error and error rate 
are different, the controlled variable is 
moving toward synchronization. Jf the 
actuating signal, M(t), is switched at 
precisely the correct time, a deadbeat re- 
sponse will result. Thus the predictive 
system logic calls for repetitive computa- 
tions to determine if future synchroniza- 
tion would occur if the polarity of the 
actuating signal were switched at the 
presentinstant of time. For this informa- 
tion a fast-time scale model is used to 
predict the future response of the con- 
trolled variable and a fast-time scale 
predictor (or extrapolator) is used to 
predict the future values of the input. 
Since both the model and extrapolator 
are synchronized and predict at the same 
fast-time scale, a subtraction process is 
used to determine when the future error 
rate is equal to zero. When this occurs, 
(at time fy), the sign of future error is 
stored and used as the basis for selecting 
the sign of the actuating signal M(#) for 
the next instant of time. If this value of 
future error falls short of synchronizing, 


‘puted, at 2 their signs are determined, at =| 


no actionis taken. If, on the other hand¢ 
the future error overshoots synchronizing ¢ 
the polarity of M(t) is switched imme4 
diately. Thus the switching action wili 
occur properly if the polarity of the acq 
tuating signal is made the same as the] 
future error signal evaluated when the 
future error rate equals zero. An ex 
ample of this type of operation is alsa 
seen in Figs. 2 and 3. Predictions fo 
trial switchings are made at points [ 
through j. When the response reacheg 
point j, the logic indicates future syny 
chronization will occur if M(t) is switchee| 
immediately. Finally, it may be seer} 
that synchronization does occur at pointy 
k. | 
When a digital computer is used td 
perform the logical processes just ded 
scribed, the computer flow diagram is very 
similar to the logic flow chart. 
such finer points as establishing predic 
and sampling rates need there be any 
difference. On the other hand, * 
circuitry for the analog computer logiei 
not as obvious when given only the logic 
flow chart. Fig. 5 shows a typical fo 
of analog logic. | 
By comparing Figs. 4 and 5 it is ap? 
parent that the mechanization of Fig §} 
performs the logic of Fig. 4. Thus al 
1, present error and error rate are com 


the signs of error and error rate are com; 
pared. If the signs are the same, switch} 
S; is energized and at 7b the sign of thal 
error is used for the sign of the actuating 1 
signal applied to the controlled cyst 
If the signs are different, at 4 the presen — 
state condition switches open and at £ 
the fast-time model predicts the future 
error trajectory. When the predictect 
error rate becomes zerO, at 6, the sig | 


system. Thus the analog logic mecha) 
nizes the logic flow chart. 
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ig. 6. Typical aircraft landing responses to initial errors in elevation 


and rate of descent 


It was noted earlier that a variable pre- 
ction rate is an excellent method of 
ducing the size of the limit cycle about 
e point of synchronization. This mode 
operation in an analog simulation re- 
lires a suitable time delay to insure a 
If-triggering action when prediction is 
scessary. Fig. 5 was described for this 
‘pe of operation and the delay unit fol- 
wing step 6 make possible a reiterative 
op between steps 6 and 4. Thus the 
1alog simulation can easily work in the 
ariable prediction mode with only the 
, switch pickup time (6 milliseconds) pro- 
ding sufficient delay. 

Thus it can be seen that the mechaniza- 
on of the control logic can be accom- 
ished by either digital or analog means. 
1 theory, both methods obtain the same 
sults since they use the same set of 
gical criteria. In practice, the digital 
ymputer is more accurate and the ana- 
g computer is considerably more eco- 
mmical. 


imulation Results 


As a typical example for which the pre- 
etive-control method can be used to 
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Fig. 7. A\ircraft landing response to disturbance 


advantage, a simplified aircraft landing 
problem was chosen. Fig. 3 depicts the 
case where the reference is an exponenti- 
ally decaying curve and where the aircraft 
was initially too high and descending too 
slowly. This particular problem was 
chosen because the reference is time-vary- 
ing and the results can be readily inter- 
preted for the various effects of initial 
condition errors, model inaccuracies, load 
disturbances, vehicle acceleration con- 
straints, and magnitude variations in the 
actuating signal. Such effects as ground 
cushion, nonlinear lift and drag, etc., 
were not included in the description of 
the aircraft kinematics. 

Before describing the results obtained 
with the use of a predictive control for 
this aircraft landing problem flying a 
prescribed exponential flight path and 
having a +2 g lateral acceleration limit, 
it is interesting to note some other con- 
trol approaches that might be employed 
to solve this problem. 

These approaches might include use of a 
linear system, or a final-value control sys- 
tem in which the touch-down point is the 
final value, or some sort of a nonlinear 
adaptive system. 


ELEVATION ERROR - FEET 
fS) 


20 «Fig. 8. Aircraft 
landing error re- 
sponse for in- ay 
accurate model 
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If one were to use a linear-control sys- 
tem, one finds that the actual flight path 
lies under the desired exponential path. 
This, of course, means that the plane 
touches down before the proper landing 
position and does so at a higher than 
desired touch-down rate. The presence 
of a rather low loop gain and a limited 
lateral acceleration aggravates these con- 
ditions. Although one might bias the 
desired path to compensate for the con- 
trol-system lags, the amount of bias will 
change with aircraft speed and other 
control parameters. 

The use of a final-value control-system 
approach means that the controller must 
consider at each instant of time, the 
plane’s path from the present to the time 
of touch-down. This requires that the 
amount of information handled by the 
controller be greater than for the pre- 
dictive control which only need consider 
the plane’s flight path from the present 
until the time in the future when the error 
rate is zero with maximum control effort 
used throughout. Using the final-value 
approach, there is a possibility that wind 
gusts may later throw one off-course such 
that the acceleration required to land 
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properly is greater than that which is 
available. 

Use of a nonlinear or adaptive-control 
system has been shown to improve the 
flight-path performance over a similar 
linear system but the performance still fell 
short of the results described later for a pre- 
dictive system. Naturally, some non- 
linear systems are better than others and 
extensive effort was not made to obtain a 
“best” system. Frequently, considerable 
effortis required to obtain a suitable adap- 
tive control whereas the predictive-control 
design means tend to be more direct and 
require less time. 


SECOND-ORDER CONTROLLED SYSTEM 


‘In its simplest form, the controlled sys- 
tem of Fig. 1 for the aircraft landing 
problem was taken to be of the form: 


Ga 28 
Mis) s(Ts-+1) 


Fig. 6 shows the type of response obtained 
for the case where the aircraft initially 
has initial condition errors of +20% 
in elevation and +20% in rate of descent 
and for a model that has the proper 
time constant and gain. For purposes 
of simplicity, the reference was taken as 
an exponential in time. This assumes 
that the horizontal velocity remains con- 
stant. Examination of this figure shows 
the fast synchronizing and close follow-up 
characteristics of predictive control. 
Fig. 7 shows the rapid recovery from 
load disturbances. It is characteristic 
of a predictive control to provide fast re- 
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Error responses when neglecting various third-order system terms in model 


sponse characteristics in removing initial 
errors and disturbance effects and to fol- 
low closely time-varying input commands 
because it is capable of using full power 
capabilities for even minor corrections. 
Since the logic network is not affected by 
past history, except through the present- 
state conditions that are repetitively ap- 
plied to the model, predictive controltends 
to be an optimum type of control in the 
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ELEVATION ERROR- FEET 


sense that it tries to reduce an ertc 
and its derivative to zero in minimur 
time. Since the parameters of the con 
trolled system will, in general, not t 
known exactly, consideration was given 
to the effect of inaccuracies in the tim 
constant and/or gain of the second-orde 
controlled system. Fig. 8(A) compare 
the error response of the controlled sy, 
tem for the model having its time constan 
50% fast, 50% slow, 100% slow, as we 
as the same value as the actual controle 
system with the initial elevation 209) 
low. | 

Fig. 8(B) illustrates the error respons] 
when the model has an inaccuracy ¢| 
+20% in gain and time constant fd} 
comparable initial conditions of Fig. 8(AN 
It can be seen that these cause rele] 
tively little difficulty. The basic effec 
of model inaccuracies can be summarizes 
as follows: 2 


= | 


1. If the system is more sluggish (lowd 
gain and/or slower time constants) thay 
the model assumed, an overshoot whos 
magnitude is related to this degree of in 
accuracy will occur. 


2. If the system is faster (higher ail 
and/or faster time constants) than th 
model assumed, a deadbeat system having 
the basic response of the model will result 


If the exact system characteristics 
either unknown or time variable, it is bes 
when designing for minimum settling 
time and minimum overshoot, to hay 
the system somewhat faster than tk 
model. 
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Fig. 10. Error responses for fourth-order system and various second-order models 
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GHER ORDER CONTROLLED SYSTEM 


[he representation of a controlled sys- 
n as a second-order system is generally 
rather simplified assumption. In the 
tilts that follow, the controlled systems 
asidered have transfer functions of the 
m: . 
com K 

— M — s(Tes+1)(T3s+1)( Tus +1) 


(2) 


ere 7,;=0 represents a third-order sys- 
n. 

It is possible to control in a near 
timum manner such a third- or fourth- 
ler system with only a second-order 
del. It should be pointed out that for 
third-order system a phase-space form 
model and logic is necessary for true 
timum performance; however it was 
cided to investigate what penalties 
are would’ be by keeping the much 
apler second-order logic-control system. 
It can be seen in Fig. 9 that the effect 
adding an additional significant time 
ustant to the system with no model 
mpensation, is to increase initial 
ershoot and settling time and, if 
ried to an extreme, even causes sus- 
ned oscillation of the system in the 
ady-state condition. Thesmall steady- 
ite errors shown in Figs. 9, 10, and 11 
re caused by drift in the unstabilized 
alog computer used to evaluate these 
os. \ 

By increasing the dominant time con- 
int ‘of the model it is possible to com- 
nsate for the effect of an additional 
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system time constant as shown in Fig. 
10. Thus, as the model time constant 
gets slower, the responses are less oscilla- 
tory although somewhat more sluggish. 
In Fig. 11 the model time constant was 
made twice as large as the dominant time 
constant of several system configurations. 
From the transient responses it is seen 
that all systems have the same initial 
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transient response even though the trans- 
fer function of the controlled system is 
quite different. However, it is also seen 
that even though the fourth-order system 
has an almost deadbeat initial response, it 
eventually remains in a sustained limit 
cycle. 

Finally, in Fig. 12, the effects of com- 
pensating a third-order system with a 
third-order model using second-order 
logiccan beseen. The stabilization effect 
on the response by including an additional 
time constant in the model with present- 
state conditions applied to it is very 
noticeable in this figure. 


One of the decided disadvantages of a 
bang-bang type of control system is the 
excessive wear and power drain caused 
by the power element banging back and 
forth, even when the error is small. 
For this reason the magnitude of the 
actuating signal was reduced by 75% 
when the error and error rate became less 
than 8 feet and 8 feet/second as shown in 
Fig. 13. The results obtained were what 
would be expected in the sense that full 
power capability would be applied during 
large errors while smaller power signals 
would be applied during times of small 
errors. If at any time a load disturb- 
ance or change in input would increase 
the error beyond this narrow band, full 
power would be immediately applied for 
correction. 

It was seen from this relatively simple 
example of an aircraft landing problem 
that the predictive control will: (1) follow 
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Fig. 12. Error responses for fourth-order system with second- and third-order models 
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variable actuating signal 


closely a time-varying input; (2) react 
in a near optimum manner to correct any 
errors that might be caused by either 
initial conditions or load disturbances; 
(3) work successfully for a wide varia- 
tion in model inaccuracies; and (4) 
operate at reduced power capabilities 
when the errors were less than a specified 
minimum. 


Conclusions 


The experience gained over the last 
2 years in applying predictive control to 


various control applications indicates that 
this form of control has definitely desira- 
ble characteristics for many high perform- 
ance control systems. Because a maxt- 
mum actuation signal is used, the maxi- 
mum force capabilities of the controlled 
system are used to remove initial condi- 
tion errors and disturbance effects. This 
allows rapid synchronization with and 
close correspondence to a wide range of 
command signals. The conventional sta- 
bility considerations do not seem to be 
a major consideration. Rather, the basic 
limitation lies in the speed of computation 
and in the accuracy of the measurement 
of the present-state conditions. 

A major disadvantage of this form of 
control lies in the need for a fast-time 
simulator and a fairly elaborate control 
logic. A trade-off must be made between 
rapid, accurate response and equipment 
complexity. 

The simulator may be of either the 
analog or digital form. The analog ap- 
proach offers economy, small size, reli- 
ability, and fast response. The digital 
approach offers accuracy and is most prac- 
tical when the computer is used for a 
multiplicity of purposes. 

The logic used in this paper represents 
an improvement over past techniques in 
that it is somewhat simpler and uses a 
variable-frequency repetition rate for pre- 
diction which is highest when the con- 
trolled variable is close to the command 


Discussion 


J. F. Coales (University of Cambridge, 
Cambridge, England): This is an excellent 
and concise description of the single-variable 
predictive-control system. The simplifica- 
tion of the logical rules and the provision 
for resetting as soon as the error rate is re- 
duced to zero are valuable refinements. 
In our early work we had established that 
the accuracy of the plant simulation could 
be quite crude, mainly because, as the in- 
stant of changeover is approached, the time 
over which the prediction is made is re- 
duced so that at the only time that is opera- 
tionally important, the actual instant of 
changeover, the prediction error is small 
even when the prediction is crude. It is 
most valuable to have quantitative results 
on this important matter. Because of the 
short prediction time at the actual instant 
of changeover, the prediction of future 
input and input rate can also be quite crude, 
so that simple linear extrapolation is 
usually adequate when the input is not pre- 
determined as in the controlled landing case 
when the input (or desired output) is a 
random process, it is possible to design 
quite simple predictors to give the most 
probable future values of input and input 
rate, provided that the autocorrelation func- 
tion of the input is known. It is probable 
that, for this purpose, the autocorrelation 


© 


Sd 


functions of most inputs met in practice 
can be adequately represented by exponen- 
tial functions, in which case quite simple 
Laguerre function predictors can be used to 
give the most probable future values of the 
input and its derivative. 

We have recently done some work with 
predicted changeover when the load is sub- 
jected to random disturbances such as 
gusts of wind on a radar aerial. Such dis- 
turbances can result in great changes in the 
trajectories in the e/é phase plane. This 
can introduce large errors and unsymmetrical 
limit cycles which, when smoothed out, re- 
sult in a steady-state error. The use of 
Laguerre function predictors can become 
important under these conditions because, 
with a large unknown torque opposing mo- 
tion, the time from changeover to alignment 
may be increased several times and the ac- 
curacy of prediction must therefore be im- 
proved. 

In Cambridge we are now turning our at- 
tention to multivariable systems and it is 
hoped that it will be found possible to use 
simple logical rules to drive these to the 
optimum as fast as the constraints allow. 
If the systems can be made noninteracting 
by the techniques of Boksenbom and 
Hood,'s? there is really no problem, since 
all controls can be driven to their desired 
values independently and simultaneously. 
This, however, is not necessarily the best 
procedure because the networks required to 
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signal. This effect improves the speq 
of response and accuracy of the syste ma mi 
reduces the magnitude of any limit eye 
that may-exist. si 

It was found that the model need n 
be an exact replica of the plant to provic ic} 
good control. Time constants and gain 
can be in error by 2 to 1 and more witi 
out excessive loss of performance. ; 
general, to avoid overshoot, the fast-ting 
model should be made to simulate a § sy 
tem slower than the actual controll id 


system, | | 
| 


When approximating higher 
systems with a lower order model, 
model may be made more effective w 
its most dominant time constant is prij 
portional to the sum of the time constanij 


=» 


in the higher order controlled system 
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make the parameters noninteracting ma 
make the system inherently more slugg! 
The requirement is therefore to find me 
of applying predictive control to interactii 
multivariable systems by the applicati4 
of simple logical rules, and it is this that 
are investigating. 
It does not seem to be generally appp 
ciated that the limit cycle, which is ine “i 
in 


t 
} 
4 


i” 


uloxe 


ble in “bang-bang’’ control of this 1 
need only appear at the output of the re 
There is no objection to putting a Filth 
between the relay and prime mover 
smooth out the oscillation, but this will, 7 
course, introduce a delay. Howeverg ti 
authors have shown that adding an add 
tional delay has little effect, even if mi 
simulated, provided it is only about 7 y 
that of the dominant time delay. ; 
to reduce chatter to about the same : 


as that in a linear system, it is nece 
only to insure that the prediction time a 
delays around the relay circuit are 
ciently short to keep the limit-cycle pet 
short enough that it can be smoothed % 
with an integrator of a time constant 1/2 
that of the dominant time delay of the plan 
In general, the filter can be at a low pow 
level and, if its characteristics are inclue 
in the high-speed simulator, it will have 
effect on performance unless the limit « bi 
is too great, in which case it may increase. ti 
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in C. Lozier (Bell Telephone Laboratories, 
ippany, N. J.): The authors are to be 
gratulated on their success in mechaniz- 
a practical form of optimum switching 
h a minimum of analog computing 
ipment. There have been many papers 
the theory of optimum switching, all 
olving about the idea that an uth order 
tem with a limited driving function can 
duce a deadbeat step response (the error 
l all of its n-1 derivatives simultaneously 
ng to zero) in minimum time with n-1 
ersals of the maximum driving function. 
wever, as far as I know, no one has pre- 
ted a practical realization of the general 
imum switching case for higher order 
tems. 
The simplifying aspect of approach used 
the authors is the concentration on null- 
just the error and its first derivative, 
ich is done here with a single reversal of 
driving function. It would appear that 
a, practical system where noise is present 
Ss simplified optimum should realize most 
the advantage gained from the optimum 
itching approach. I would be interested 
having the authors’ comments on this 
nt. 
Vly associates and I have been doing 
lar work using a digital computer 
chanization and have found it advan- 
feous to be able to predict the system re- 
mse with both polarities of drive. This is 
pful, for example, in cases where the first 
que reversal was premature. (A pessi- 
stic model will cause the system to 
itch early.) Under these conditions the 
tem will undershoot with a correspond- 
ly poor transient response unless it can 
ik with the same polarity of drive that it 
using, see that it will cross the error axis 
the phase plane short of the origin, and 
itch the torque back again. Predicting 
th the opposite polarity of torque will not 
lp in this situation for the simple reason 
ut the predicted trajectories will not cross 
s error axis, and hence the switching cri- 
ion has no information to go on until it 
too late to prevent the undershoot. The 
rves for’ T=7.5 and T=10 in Fig. 8(A) 
monstrate this effect. Perhaps this 
lity to look both ways could be added 
thout too much complication. 
One of the most interesting aspects of an 
timum-type system is its small-signal 
rfformance after the transient phase is 
er. This criterion is one that will work 
th any amplitude of input signal. Fur- 
srmore it is one that seeks to drive the 
stem so that the input and output meet 
igentially; one would thus expect the 
all-signal performance to be very good. 
fact, if the input were constant, such a 
stem could conceivably oscillate about the 
sired value at a frequency determined 
marily by the response characteristics of 
» high-speed model. It is clear that such a 
forcement of the stability considerations 
m the properties of the actual system is 
ated to the model’s ability to predict. 
Unfortunately, with this simple criterion 
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the prediction is only to the crossing of the 
error axis in the phase plane, and thus be- 
comes negligibly short when the error is very 
small. Consequently one would expect 
that some intermediate oscillating condition 
would obtain between that dictated by the 
dynamics of the system and that dictated 
by the dynamics of the model. Further- 
more one would expect that this effect would 
be greater in higher order systems where 
the need for a lead in the switching times is 
greater. A full-scale analysis of this prob- 
lem would be very difficult but it would be 
interesting to have the authors’ experience 
with this small-signal mode of operation. 
In particular, it would be interesting to 
know what oscillating frequency was ob- 
tained for the fourth-order system and how 
it compared with the oscillating frequency 
of the system without prediction. 


L. F. Kazda (University of Michigan, Ann 
Arbor, Mich.): The authors are to be 
complimented on their excellent paper on 
the application of predictive control to 
practical systems. They have taken known 
concepts and have investigated the many 
ramifications of these concepts which are 
necessary before this approach can be ap- 
plied to a practical system. Review of the 
paper reveals that the authors have not 
defined the type of system which they 
considered amenable to their method of 
analysis. For example, are pure stochastic 
inputs prohibited? What happens to the 
performance of this type of system if the 
aircraft suddenly receives an extremely 
large sharp-edged-vertical gust, which may 
result in a major change in the performance 
characteristics of the aircraft? Would the 
performance deteriorate or would the system 
become unstable? In short, has stability in 
the large been given any detailed study in 
these applications? Given a particular ap- 
plication, it would have been helpful if the 
authors had outlined a format which could 
be followed by others in applying the ideas 
developed in this paper to a new problem. 


Frederick A. Russell (Newark College of 


Engineering, Newark, N. J.): The authors 
are to be congratulated on an interesting 
and well-presented description of an ap- 
plication of a predictive-control system. 
This is a form of automatic control which 
holds great promise for further development, 
especially in situations where a computer is 
required for other purposes. 


While the authors discuss the attractive 
characteristics of this system in resisting 
load disturbances, they give only passing 
consideration to the effects of noise in the 
reference signal, and perhaps some comment 
on this phase of the problem would be 
useful. I have worked with a very similar 
control system in which reference signal 
noise could be the determining factor in the 
size of the limit cycle. When e(t) and é¢) 
have been reduced to small values, it can be 
seen from Fig. 3 of the paper that an im- 
pulse of noise can readily cause an incorrect 
switching decision on the part of the logic 
block of Fig. 1. 

If random noise is present in the reference 
signal, r(t), it will be amplified in 7(¢) by 
differentiation. If we imagine the control 
system stationary at zero error and zero 
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error rate, it is easy to see that the control 
logic will be actuated by the spurious noise 
into making random decisions as to the 
direction in which the output should be 
ordered to move. An average of a large 
number of these excursions can be found 
for a given amplitude of noise, and this mo- 
tion can be called a limit cycle. 

The amplitude of this limit cycle due to 
noise depends also on the torque exerted by 
the bang-bang controller, and one great 
advantage of limiting the driving torque 
to a small value as the error and error rate 
become small is that the system will havea 
smaller limit cycle due to noise, although 
it becomes more sensitivie to load dis- 
turbances. 

Filtering, especially of 7(t), can be used to 
reduce the amount of noise. However, the 
filtered reference signal is a delayed version 
of the reference signal which the control sys- 
tem is required to follow; hence filtering also 
tends to cause incorrect decisions on the 
part of the control logic, because decisions 
are being made on “‘stale”’ data. 

We conclude then, that in the presence 
of a noisy input, a compromise must be 
made among three parameters: the amount 
of filtering, the amount of torque reduction, 
and the sensitivity to load disturbances. 
Optimization of these factors is an area 
which should be very attractive for further 
research. 


W. M. Gaines (General Electric Company; 
Phoenix, Ariz.): The use of a mathematical 
model in process-control computers has 
proved to be a very effective means of ac- 
complishing process optimization and con- 
trol. This paper is another excellent 
example of how a mathematical model can 
be used effectively to improve the perform- 
ance of a dynamic system. 

It is useful to separate the types of 
optimization into two classes: (1) steady- 
state optimization; and (2) transient opti- 
mization. 

In “steady-state” optimization the objec- 
tive is to determine the optimum settings 
of the process-control variables to obtain a 
long-term “optimum” performance. Once 
the model has been developed there are a 
number of methods of determining the best 
performance. Linear programming can be 
used when performance can be approxi- 
mated by linear relationships. Gradient 
methods (methods of steepest ascent or hill 
climbing) or other iterative techniques can 
be used in nonlinear cases. (Such methods 
have been adapted to utilize the process 
itself as the model and by repeated per- 
turbation and measurement to achieve 
optimum performance.) A review of the 


_ literature will indicate that the majority of 


applications of process-control computers 
has been in the area of steady-state opti- 
mization. 

In transient optimization we are in- 
terested in minimizing the cost during a 
transient. This may be a transient intro- 
duced by external disturbances or it may be 
the result of a long-term management deci- 
sion to change product characteristics. 
Mathematically the result is the desire to 
minimize an integral of the form shown in 
equation 3: 


I(m, h)= Sar (Qn, Ins Mn, Mn), r)dr (3) 
where Q, and M, are the desired output 
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INPUTS OUTPUTS, 
Daca PROCESS dn 


Fig. 14. Nomenclature for equation 3 


and input. The nomenclature is defined in 
Fig. 14. 

To date there has been very little general 
analysis of this type of problem which has 
led to practical on-line solutions. Dynamic 
programming offers a method of solution but 
in most practical situations results in a very 
formidable requirement on computer speed 
and memory. The work of Dr. Merriam! 
provides a possible solution in those in- 
stances where a more restrictive class of 
transient optimization can be met. The 
greatest progress, however, has been made 
in the restrictive class of problem wherein 
the objective is to minimize the elapsed 
time of the transient. The paper is a signifi- 
cant contribution as is the work which has 
been done at the Case Institute.2 Hope- 
fully, we will see this type of effort con- 
tinued and further progress in the analysis 
of the more general problem of transient 
optimization. 

The paper has several subtle implications 
which can be of considerable use in indus- 
trial process control; and I would like the 
authors to comment upon the following 
items: 

1. In many industrial processes the input 
variables are not changed in the stepwise 
manner indicated in the paper. They are 
limited to rather slow rates of change either 
because of physical constraints upon the 
variables themselves or because of process 
limitations which make it mandatory to in- 
crease or decrease the variables more 
slowly. It would seem in instances where 
this is the case that the predictive technique 
could still be utilized effectively. I would like 
the authors’ comments upon any experience 
they might have had with cases where input- 
manipulated variables are not changed in 
a stepwise manner. 

2. In most process-control situations we 
are concerned with a multivariable problem 
and there is usually considerable coupling 
between the variables. It would seem that 
in this instance the dimensionality of the 
process model would similarly have to be 
increased. While it is still possible in this 
instance to represent the dynamics by 
second-degree equations, in all likelihood 
the model would now constitute a set of 
simultaneous second-order differential equa- 
tions. I would like to ask the authors how 
this particular problem can be handled and 
what complexity it may introduce. 

3. Difficulty can be encountered in on- 
line applications of process-control com- 
puters in obtaining sufficiently fast sampling 
and real-time computational rates to obtain 
the degree of stability desired. (Such 
sampling rates can be determined in linear 
operation by the use of z-transforms or other 
linear methods of stability analysis.) It 
appears quite possible that by using pre- 
dictive-control systems the sampling rate 
required during the course of a transient can 
be drastically reduced. The amount of re- 
duction appears to be a function of the cor- 
rectness (or exactness) of the process 
model; i.e., if the model were perfect it 
would be necessary to sample the output 
only once. I would like the authors’ com- 


136 Chestnut, Sollecito, Troutman—Predictive-Control A pplication 


ments on this particular potential of the 
predictive technique. 


REFERENCES 


1. Use oF A MATHEMATICAL ERROR CRITERION 
IN THE DeEsIGN oF ADAPTIVE CONTROL SYSTEMS, 
C. W. Merriam, III. AZEE Transactions, pt. IL 
(Applications and Industry), vol. 78, 1959 (Jan. 1960 
section), pp. 506-12. 


2. OptimizING CONTROL OF A CHEMICAL PROCESS. 
D. P. Eckman, I. Lefkowitz. Control Engineering, 
New York, N. Y., Sept. 1957. 


A. M. Hopkin (University of California, 
Berkeley, Calif.): The development of the 
general concept of ‘“optimum’’ relay control 
systems has resulted in much more paper 
analysis than in actual hardware. This 
situation probably exists because, in spite 
of the attractive features of the optimum 
relay type of system, its practical realiza- 
tion usually seems to involve too much 
complication and cost to compete econom-~ 
ically with alternative systems. With 
this fact in mind, the authors are to be 
congratulated on seeking to find systems 
which combine sufficient simplicity with 
sufficient performance to be practical. 

The authors point out that they are fol- 
lowing up the work of Noton and Coales 
(reference 1 of the paper) and of Lozier. 
It might be of interest to note that about 8 
years ago Otto J. M. Smith suggested the 
use of a somewhat similar procedure. This 
suggestion was not widely circulated, but 
one of his students did write a thesis con- 
cerning an analog computer study of the 
idea,1 and one of this discusser’s students 
later wrote a thesis on a similar topic.? 

The authors are quite careful to indicate 
that they are approximating the actual 
system behavior by a simple model, and 
that the criterion which determines the 
switching-control logic applies strictly only 
to the control of systems described by 
second-order equations. I would like to 
consider the example given with the purpose 
of indicating situations where the second- 
order criterion, based on error and error rate, 
might well prove inadequate. 

In the example in the paper, the aircraft 
transfer functions, equations 1 and 2, were 
assumed to have only real poles and to 
possess no zeros. I feel that possibly a 
more realistic transfer function for the 
elevation control of a tail-steered aircraft 
would be of the approximate form: 


_C(s)_K (s+Z1) (s+Zz) 
CO RO) et ee as 


where it is assumed that the elevator is 
positioned by a hydraulic cylinder. For 
equation 4 it is possible that one of the zeros 
may be in the right-hand s-plane, and & 
may be 0.2 or less, resulting in lightly 
damped complex poles. 

For this case both the presence of lightly 
damped complex poles and the presence of 
zeros, particularly those in the right-hand 
half-plane, render the higher order system 
less likely to give good behavior with a 
switching logic based on a second-order cri- 
terion. These two factors also require that 
care be used in selecting a simple fast-time 
model to predict system behavior. The 
procedure followed by the authors of the 
paper of carefully checking system behavior 
say on an analog computer, would most 
certainly be called for. 
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Irmgard Fliigge-Lotz (Stanford Universi | 
Stanford University, Calif.): After readij 
this very interesting paper, I have two qu 
tions which concern an extension of ¢ 
method. . | 
1. I wonder how predictive contr 
works for a system which has complex ; 
real and complex poles. For the secon 
order case the optimum switching curve hij 
a rather simple character for real charat| 
teristic roots. For complex roots, howeva 
this curve is more complicated (see Bl 


Tsien!). I would like to know how pred 
tive control works in such a case, and pa} 
ticularly, whether it will be necessary | 
limit the magnitude of the initial distur 
ance in such a way that only the first an) 
of the optimum switching curves near ti 
origin of the phase plane are used. ; 

2. What experience do the authors hag 
with systems with zeros in the transfs} 
function? 7 
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J. J. Jonsson (Polytechnic Institute | 
Brooklyn, Brooklyn, N. Y.): The con 
bination of prediction and optimum Huy 
time using a bang-bang controller shouil 
make this system extremely attractive - ; 
future applications. It therefore represe: ny 
a valuable contribution. Two comment | 
seem pertinent. 

il, Fig. 38 shows the actual trajeciee 
from 7 to k. The Paper seems to indicati} 
and it would be true in practice, that t 
model trajectories cannot be accomplishal 
in zero time. Therefore, if the system 
state 7 when the model trajectory is i 
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Fig. 15. Model and system trajectories wil 
null error of E 


1—Model trajectory just short of origin 

2—Model trajectory commanding switch of t 
controller 

3—Actual trajectory 
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. 16. Model and system trajectories with 
null error of E/2 


See subcaption of Fig. 15 
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17. Discrete time model trajectory 


ced, the actual system will continue up- 
rd and to the right until the model 
ches point k and the system logic com- 
nds a change of controller sign. In fact, 
he model trajectory passed just to the 
of the origin by a very small amount, 
entire new model pass must be made 
ore switching can eccur. This will result 
aA maximum null error of E as shown in 
. 15. For a given model speed, E may be 
uced by a factor of 2 with a slight 
dification as indicated in Fig. 16. If 
system is at state 7, the actual trajectory 
sredicted ahead on a fast-time basis to 
nt j’, mid-way between 7 and m. A 
tehing command will result from model 
jectory 2, giving trajectory 3 with a 
ximum null error of H/2. Extending the 
> j’ prediction beyond half the j to m 
france does not seem advisable since the 
sibility of a null error of less than — E/2 
n exists. 

. As so many others, this paper is pri- 
rily concerned with large signal perform- 
e and, therefore, the challenge of the 
0 Joint Automatic Control Conference 
sents itself again, namely: How may 
describe or anticipate small-signal be- 
dor near the origin of the phase plane? 
isider the model trajectory of Fig. 17. 
> shortest possible fractional part of this 
lectory, j to p, represents an elapsed 
e T, the entire trajectory being made 
of » parts. Such a limiting value of T 
| exist in a practical sense in the case of 
ital simulation, and probably for analog 
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switching 
curve 


Fig. 18. System switching curve and possi- 
ble trajectories 


Computer decision synchronized with required 
system switching 

+ Proper switch point 

@ Other possible switch points 


simulation also. This implies a maximum 
limit cycle frequency for the actual system 
of 1/2T; that is, at its highest possible 
frequency of steady-state oscillation the 
model must command two symmetrical 
controller switchings per system cycle. The 
actual system trajectory for this condition 
is shown in Fig. 18, curve 1, where 


fa= ; =the model frequency 


f=actual limit cycle frequency 


The system under consideration is third 
order of form 


A 
s(s+a)(s-+b) 


It is conceivable that the system may 
also lock in at ‘ 


pois ie dae etc. 
AGS 


to give the symmetrical trajectories shown 
in Fig. 18, switching decisions being made 
at the end of every second, third, fourth, 
etc., model cycle T, respectively. Discrete 
orbital modes, not necessarily circles, are 
illustrated together with the proper steady- 
state switching curve. 

The effect of sampling the actual system 
state variables must also be taken into ac- 
count, i.e., KG(s) is replaced by a sampling 
transfer function KG*(s) where 


+00 


KG*(s)= - e) KG(s+j nwm) 


n=—o 
with 
T =sampling time 
®m=model frequency 


The effect will be that of requiring a shift, 
primarily in phase, of the switching curve of 
Fig. 18 to anticipate the sampling influence. 

In general computing, sampling, and logi- 
cal decision to switch is done on an arbitrary 
time base, with discrete switchings being 
possible at points 1, 2, 3, etc., with respect 
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Fig. 19. Discrete computer decisions asynch- 
ronous to required system switching 


Fig. 20. Continuous switching curve as 
modified by sampling and relay holding 


A—Switching curve for continuous system 
B—Switching curve with effect of sampling 
and controller hold included 


to the proper switch points X of Fig. 19. 
Therefore, there is no assurance that the 
control logic is synchronized to the points X 
of the phase plane. Switching should not 
be done at point 1 because this is too early 
and yet point 2 is too late. On the average, 
the time delay is one half of a computing 
cycle. If the delay is referred to the actual 
system then 


Average controller holding delay 
=99° 
fm/2 


That is, if f=fm/2 the delay is 90 degrees; 
for f=fm/4 the delay is 45 degrees, etc. 
Such a delay must be anticipated by shift- 
ing the switching curve of Fig. 18. 

The complete system transfer function 
must include the original transfer function, 
the effect of sampling, and also the effect 
of controller holding. While the gain is af- 
fected very little, a significant change 
occurs in the phase characteristic. In fact 
at f=fm/2 the transfer function for the 
third-order system cited above lies at 270 
degrees. Curve B of Fig. 20 indicates the 
new average switching contour. It is im- 
portant to note that this is an average curve 
since the actual controller holding delay 
may vary anywhere from zero to the full 
time T, or zero to 180 degrees. 

It is therefore apparent that switching 
criteria near the origin of the phase plane 
may be quite different from those used for 


137 


large signals. What this should be is not 
altogether clear because of the arbitrary 
time relationship between the discrete 
model behavior and the fixed switching 
points of the actual system. The (7—1) 
switching usually associated with an m-order 
system is replaced by a single switch in the 
é—e phase plane in this paper when driving 
to zero from large errors. For third-order, 
and higher, systems one would hope that 
most of the stored energy in the system 
would have been removed with this tech- 
nique which indeed represents a very com- 
mandable practical approach to the n- 
dimensional switching problem, but the 
same reasoning does not apply at the origin. 
Indeed (unless zero forcing function is ad- 
missible as a controller possibility) it would 
be desirable to have the system oscillate 
with minimum error, namely, trajectory 1 
of Fig. 18. But what will be the actual path 
in the e—é phase plane? It is suspected 
that the mode of operation will skip from 


one trajectory of Fig. 18 to another in some 


statistically random fashion. Additional re- 
search is necessary to estimate the limits of 
error and error rate. There will be some 
tendency to ignore this problem by revert- 
ing to a linear mode of operation near the 
origin or else simply increasing the model 
simulation frequency. The first solution 
still leaves the above questions to be 
answered. The second solution has economic 
implications and physical limitations since 
systems are now in existence where com- 
munication between computer and the 
actual control system is done with discrete 
pulses at a repetition rate which cannot be 
economically increased. 


H. Chestnut, W. E. Sollecito, and P. H. 
Troutman: We are pleased with the large 
number of interesting and worthwhile dis- 
cussions on our paper. In a number of cases 
the discussers have had a long and scholarly 
association with this problem, and we are 
delighted they are willing to share their ex- 
periences and comments with us. 

We agree with Mr. Coales that the bang- 
bang aspect of predictive control is not 
essential when the errors are small. Use of 
a filtering scheme such as Mr. Coales sug- 
gests, or of a decreased bang-bang ampli- 
tude signal, or perhaps going to a very small 
linear control range in the vicinity of zero 
error are all possibilities that might be em- 
ployed for the small-signal case. 

We are interested to hear of the work that 
Mr. Coales is doing on the subject of dis- 
turbances which may result in large changes 
in the trajectory. We feel that these 
random disturbances can best be considered 
as specific inputs (pulses, steps, or ramps) 
rather than that one should design for them 
on any over-all statistical basis. 

We also feel that the work on the multi- 
variable systems is extremely important and 
are glad to learn that Mr. Coales and others 
are doing work on this problem. We, too, 
are concerned with the possibility of using 
predictive controls for this application. 

Mr. Lozier has pointed out that our work 
as described deals entirely with the noise- 
free reference signal condition. We feel 
this to be valid because in many industrial 
processes the signal-to-noise ratio can be 
brought more fully under the designer’s con- 
trol and thus it is worthwhile to study the 
noise-free case first. However, as he 
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points out, the presence of noise would tend 
to complicate the logic involving higher 
order prediction terms and as such the 
simple first-order prediction could be of 
practical importance for signals containing 
noise. 

We are interested to learn that Mr. Lozier 
and his associates have been doing a digital 
computer mechanization of the predictive- 
control problems and found it advantageous 
to use prediction with both polarities. We 
feel that for some conditions it might be nice 
to try to use both polarities of drive at each 
prediction time but for reasons of sim- 
plicity we did not feel it worth the added 
complexity. We have, however, incor- 
porated in our signal polarity for the under- 
shoot conditions described by Mr Lozier, 
the ability to see that we will cross the error 
axis in the phase plane short of the origin 
and we do switch the torque back. For the 
case described in Fig. 8(A), the real problem 
is that the curves for 7=7.5 and T=10 de- 
scribe a situation in which the model is off 
by 50 and 100% respectively. 

The problem of the small-signal limit 
cycle was not significant because the oscilla- 
tions in the controlled variable were within 
prescribed acceptable limits. Our . feeling 
has been that to reduce the oscillations we 
would go to either a lower elevator force on 
the airplane or to a small range of linear 
operation about the zero point. 

In the case of analog prediction, the fre- 
quency becomes extremely high for the 
small-signal condition as was described by 
Mr. Coales. Hence, the small limit cycle is 
attributable to a combination of the predic- 
tion time and the time constants in the 
process. 

Although the significance of the closed- 
loop oscillating frequency asked for by Mr. 
Lozier was not studied in detail, it is our 
feeling that the oscillating frequency of the 
system without prediction of necessity 
would be much lower than that which will 
be realized with the prediction employed. 

In regard to Professor Kazda’s inquiry re- 
garding stochastic signals, it is our purpose 
in the predictive control to operate with sig- 
nals in which the deterministic portion is 
much larger than the stochastic one. Our 
emphasis has been upon handling individual 
disturbances as they arrive and reducing 
them as quickly as possible. Fig. 7 of the 
paper indicates the nature of the treatment 
of disturbances. Our investigation of the 
effect of model inaccuracy was directed at 
determining the effect of changes in the 
processes. Our consideration of stability in 
the large has been directed at determining 
the effect of large initial errors, disturb- 
ances, and inaccuracies in a model of the 
magnitude of those studied in the paper. 
The various values for these initial errors 
for model inaccuracies were studied for the 
purpose of determining the maximum devia- 
tions which might be expected and evaluat- 
ing the performance on the basis of these 
magnitudes. 

Regarding Professor Kazda’s request for a 
general format to be applied by others, it 
was our intent that the procedure outlined 
in Figs. 1, 4, and 5 would serve this purpose. 
Our objective was not to consider the general 
mth order system, and therefore the logic 
and general approach may appear somewhat 
limited. 

It is our feeling that in the control of 
many systems commonly encountered in 
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practice probably one or two integratici 
due to the process might reasonably be ¢ 
pected. The remaining time constan 
delays would probably be caused by 
time of response due to the controls ox 
which the designer has presumably — 
influence. Because the designer has th 
fluence, it is reasonable to expect that | 
would choose them in such a way that th 
would not become dominate in influencj 
the over-all response of the control and 
process. For this reason our attention 1 
been concentrated on relatively low-oré 
process systems and a general conside “4 
of the mth order system has been given ; 
attention. i 

It is of interest to note that the obser 1 
tions of Mr. Russell are so similar to thi 
of Mr. Coales when he discusses the | 
fects of random input signals and load c 
turbances for the small error region. 
also agree with Mr. Russell's conclusi« 
that a compromise must be made ama 
input filtering, torque reduction for sm 
errors, and sensitivity to load disturbana 

Mr. Gaines has very nicely summariz 
the state of the art concerning proc 
optimization and control. We will answ 
his questions in the order in which 7 
sented them. 

It is true that the input-manipu 
variable need not be changed in the st 
wise fashion described in the paper andi 
the predictive technique can be effectivy 
utilized where it is mandatory that 
variables increase or decrease less abrup? 
In such instances it is necessary to kra 
the time relationships of the manipulad 
variables in order to so generate them in 
predictive portion of the control. 

We are in whole-hearted agreement t4 
the multivariable problem is importa} 
It also appears that the model dim 
sionality would have to be increased, 
that another set of logic would probaé 
have to be formulated. Mr. Coales | 
indicated that his group is currently 
vestigating this area of interest. 

It appears that where the process cami} 
adequately described by linear my 
it is possible that, by the judicious w 
predictive control, the sampling rate 
quired during the course of a transient 
be reduced. This possibility exists becad 
in the calculation of the predicted respon} | 
the exponential form of the time reat 
can be used and the response at requli 
times in the future can be calculated ¢ q 
tapidly without recourse to time-consuny 
integration procedures. 

The discussions by Mrs. Pligge-Lot z 
Mr. Hopkin both raise important questil} 
in application areas in which we made | 
limited investigations. Since their excel 
discussions are closely related, both | 
be answered simultaneously. Ty] 

Using predictive controls for second-o 
systems having lightly damped complex] 
pairs should result in good swite i 
characteristics since the predictor will ii 
immediately upon entering a new quad 
whether switching should occur immedi 
(ie. the trajectory overshoots the ot i 
or be delayed until the predicted traje ct 
passes through the origin. Thus i i 


} 


q 


Third- and higher order systems 
definitely moré difficult to control. Ree 
studies were made of a fourth-order sy 34 


ares of the optimum eniiching curves. i 


ntaining two poles at the origin in 
dition to a lightly damped complex pole 
ir. It was found that simple phase- 
ane switching logic was incapable of 
oducing a stable system since the higher 
tivatives were uncontrolled; neverthe- 
3s, it was of interest to note ‘that during 
ch revolution of the diverging limit cycle 
e trajectory would pass through the 
ase-plane origin. Thus it appears ap- 
ent that additional studies should be 
ade to include higher order systems. 
A transfer function having a zero in the 
s-plane was investigated and 
little difficulty. However, 
rticular care must be taken to insert the 


proper present-state conditions into the 
model prior to each prediction. Zeros in 
the right-hand s-plane were not investigated. 

Mr. Jonnson’s discussion indicates that 
he has done considerable work in the area of 
predictive controls. The idea of shifting the 
prediction time from j to 7’ is an effective 
way of counteracting time delays in the 
system. 

Mr. Jonsson’s analysis of the small- 
signal operation has served to clarify the 
nature of the problem when digital com- 
puters are used for predicting system 
responses. For analog prediction, this 
prediction time becomes very short and 
additional considerations enter the picture 


concerning the small-signal limit cycles. 
We agree that small-signal performance is 
somewhat different in nature from large- 
signal performance. 

We also feel that introducing a linear 
range of operation for small errors is 
attractive when physically practical. For 
the case where a linear range is impractical, 
it is necessary to find other ways to reduce 
the small-error limit cycles when extreme 
accuracy is a requirement. 

In closing we should like to re-emphasize 
our appreciation for all the constructive 
comments offered by the discussers. We 
feel that their contributions have added 
considerable value to the paper. 


ANY LARGE materials handling 

machines and processing machines 
quire smooth, accurate stopping. In 
e past, accurate stopping has usually 
en accomplished by the use of hydraulic 
akes. Recently, an electromagnetic 
ake has been developed to perform the 
me functions and to overcome some dis- 
vantages of hydraulic brakes, such as 
e need for hydraulic lines and slow re- 
se of parking torque by a foot-operated 
draulic pump. The electromagnetic 
ake can be readily controlled from small 
ot devices and can be fed through 
lector bars so that the operator need 
be on the same structure as the drive 
brake. Thus, the new brake offers 
eased flexibility in application. 


inciples of Operation 

For most applications in industry a 
ast also provide a holding torque when 
drive is stopped and must apply 


ergency braking torque if failure of 
otrol power occurs. Thus, the main 
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Controllable Torque 
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requirements for a controllable torque 
brake are: 


1. With no power applied, the brake must 
exert a holding or parking torque to prevent 
the equipment on which it is mounted from 
being moved by external forces such as wind. 


2. Means must be provided for releasing 
the holding torque to permit the drive to 
run and to reapply the holding torque when 
desired by the operator. 


3. If the power supply to the drive is inter- 
rupted, the brake should apply emergency 
braking torque and holding torque. 


4, With the holding torque released, 
means must be provided to apply retarding 
torque in the amount and at the time de- 
sired by the operator. 


The first three of these requirements are 
the same as those for the spring-set 
magnetic-release brakes used widely on 
crane hoists and other industrial drives.1~* 
The fourth requirement, controllable 
torque, may be accomplished by coupling 
an electromagnet to the brake mechanism 
to apply torque in proportion to the 


Fig. 1. Electrically energized 
controllable torque brake 
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force of the magnet. Thus, the marriage 
of a spring-set magnetically released 
brake and a magnetically applied brake is 
indicated. The brake of Fig. 1 is an 
embodiment of this principle, and the 
partially sectioned views of Figs. 2 and 3 
illustrate how the mechanism fulfills 
these requirements. 


Those parts of the brake needed to meet 
the first three requirements are shown in 
Fig. 2. With no power applied to the . 
brake, a main spring applies force to a 
“whiffletree’’ linkage which in turn applies 
equal forces to the two brake shoes which 
are in contact with the brake wheel. 
Application of power to auxiliary coil PB 
causes the auxiliary armature to close and 
compress the main spring allowing the 
retractor springs to release the shoes from 
the wheel and permitting the drive to 
run. Whenever current to the auxiliary 
magnet is interrupted, either by action 
of the operator or upon power failure, the 
armature will be released and the main 
spring will set the brake. It should be 
noted that the auxiliary armature is not 
coupled directly to the shoe linkage 
mechanism but is arranged to compress 
the spring and to decouple the spring from 
the mechanism. If it were not for the 
retractor springs, the closing of the 
auxiliary armature would result in 


physical separation of the armature from 
bearing surface P on the linkage. 
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MAGNET HOUSING AUXILIARY COIL—PB 


,-RETRACTOR SPRING 


a eee, <Qaz7TITT 
Sl teers | 


p RETRACTOR SPRI NG 


PARKING TORQUE 
ADJUSTMENT 
MAIN SPRING 


Fig. 2. Spring-set magnetically-released portion of brake with frame 


and magnet sectioned 


In Fig. 3, the parts of the brake 
used in the application of controllable 
torque are shown. With the main coil 
de-energized, retractor springs hold the 
shoes away from the wheel. When the 
coil is energized, the main armature is 
attracted, the retractor springs are over- 
powered, and the shoes are set against 
the wheel. Increasing the current in the 
main coil increases the force on the arma- 
ture and the torque developed by the 
brake. 

In the usual operation of the complete 
brake, the auxiliary coil is first energized 
to decouple the main spring, and then 
the brake is essentially an electrically 
applied brake as illustrated in Fig. 3. 


Control 


A d-c constant potential controller and 
a multipoint foot switch used with the 
brake are shown in Fig. 4(A), and the 
elementary diagram for this controller is 
shown in Fig. 4(B). Operator’s controls 
also include a push-button station for 
control of parking torque application. 

In use, the operator first presses a Re- 
lease push button to close relay BRI to 
apply power to auxiliary coil PB to release 
the parking mechanism of the brake. 
After the current in the auxiliary magnet 
has increased enough to release the brake, 


current relay BR2 opens to insert a pro- 


tective resistor in series with the auxiliary 
coil. With the parking mechanism re- 
leased, the drive may be operated nor- 
mally and braking power may be applied, 
as required, by operating the master switch 
to close relays BR3, BR4, and BR5 to 
apply successively increasing amounts of 
excitation to the main coil and therefore 
greater amounts of braking torque. Field 
trials have indicated that three torque 
points are sufficient for most applications, 
but additional points may be added 
when necessary. When the operator 
wishes to apply parking torque, he pushes 
the Set button to drop out relay BR1 
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BEARING SURFACE- 


AUXILIARY ARMATURE 


SERVICE BRAKING 
COIL- SB 


MAIN ARMATURE 


which de-energizes both the auxiliary 
coil and the main coil. The same result 
is obtained when power fails. 

In many cases a normally open electrical 
interlock of relay BR1 is used in the drive- 
controller circuit to prevent energizing the 
drive motor until the parking torque is 
released. Similary, an additional contact 
is sometimes provided in the brake master 
switch to cause de-energization of the 
drive motor before controlled braking can 
be applied. Both of these provisions 
are made to avoid driving through the 
brake. 

In cases where repeated and prolonged 


teu (A) 


Fig 4. A—Brake controller and foot switch. B—Controller elementary diagram 
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Fig. 3. Controllable torque portion of brake with frame and magi 
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“—— RETRACTOR SPRINGS 


sectioned 


braking is necessary, the control has bé 
arranged to establish service dynan 
braking by the motor as the first brakh 
point and to have this operate in parai 
with the controlled braking. 

While the controller is shown with af 
switch, the brake has also been apph 
to a floor-operated crane with cont] 
from a pendant master switch. 
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Design Considerations 


The brake shown in Fig. 1 was | 
signed to have certain dimensions, } 
same as those standardized for spring-# 
brakes with 13-inch-diameter whe | 
listed in the standards of Association 
Iron and Steel Engineers (AISE) 
National Electrical Manufacturers Ad 
ciation (NEMA).*® The standardil} 
dimensions assure that the brake ~ 
mount the same way as a spring-set bri 
and that a wheel suitable for a spring? 
brake can also be applied with the c¢ 
trollable torque brake. ; 
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, 5. A—Service braking torque versus coil 
itation. B—Service braking torque versus 
time 


faximum emergency braking torque of 
brake is 550 pound-feet, the same as 
t of a 13-inch spring-set brake meeting 
SE-NEMA Standards. This torque 
be adjusted down to as low as 150 
nd-feet by turning the calibrated park- 
torque adjustment shown in Fig. 2. 
[he new brake has been designed for a 
ximum controllable torque of 1,100 
ind-feet, twice the torque rating of a 
inch spring-set brake. It is antici- 
red that this maximum torque will be 
tired infrequently and then only for 
t periods of time. Consequently, the 
to provide 1,100 pound-feet maximum 
que has been designed for a 10-minute 
ng. For most applications, a maxi- 
m torque of 550 pound-feet will be 
}quate, and a coil for this torque rating 
) a 30-minute rating. 

he main magnet has been designed 
h arather large air gap. This has two 
antages. First, it permits consider- 
y more than usual shoe wear between 
istments. As the shoes wear, the 


} 
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magnet closes more and more; if the 
magnet were finally to close completely, 
the brake would no longer apply torque 
and would require adjustment. Second, 
the large air gap reduces the effect of 
manufacturing tolerances and that of shoe 
wear on the torque developed by the 
brake This reduction in torque variation 
is highly beneficial when it is desired to 
have two brakes share load equally, or 
when brakes are mounted on each end of 
a moving bridge which is to be stopped 
without skewing. 

Equal clearances on the two shoes, when 
the brake is released, are obtained by turn- 
ing the Shoe Centering Adjustment near 
the top of the auxiliary magnet armature. 
This screw serves as a stop for the top of 
the inner shoe lever as indicated in 
Fig. 3. 

To adjust for shoe wear, a shoe clear- 
ance adjustment is provided as indicated 
in Fig. 83. The arm over the brake wheel 
is hexagonal in cross section and is 
provided with right-hand threads at 
one end and left-hand threads on the 
other so that the shoes can be drawn 
together or moved apart by turning this 
arm. 

The Manual Release shown in Fig. 2 
provides a means to compress the main 
spring without power permitting the 
retractor springs to move the shoes away 
from the wheel. With the brake man- 
ually released, maintenance work such as 
changing brake shoes and removing a 
motor armature and brake assembly may 
be accomplished easily. 


Service Braking 


Service braking torque is shown as a 
function of coil ampere-turns in Fig. 5(A). 
The lower curve shows the torque which 
is developed when the brake is adjusted, 
as recommended, for a shoe clearance of 
1/32 inch. The upper curve shows the 
torque which is developed when the shoes 
have worn to the point where shoe clear- 
ance is 5/64 inch, the point at which 
the clearance should be readjusted to 1/32 
inch. Some force is required to over- 
come the force of retractor springs when 
braking is applied, and since the force 
of the magnet increases as the armature 
moves toward the closed position in 
applying the shoes, the minimum torque 
which can be applied during service 
braking is limited. This minimum torque 
is indicated in the figure. 

Because the brake coil is highly induc- 
tive, some time is required for the current 
to reach its final value, and the applica- 
tion of torque is correspondingly gradual. 
Fig. 5(B) shows torque as a function of 
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Fig. 6. Time to release parking torque 
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Fig. 7. Time to apply torque after power 
failure 


time when full voltage is applied directly 
to the service braking coil and also when 
full voltage is applied to series arrange- 
ments of the service braking coil and 
appropriate resistances to give approxi- 
mately two-thirds and one-third rated 
torque. Maximum service braking torque 
shown on the diagram is higher than 
rated because the tests were made with a 
relatively cold coil. However, the times 
to reach all except the highest values of 
torques are nearly the same when the 


Controllable torque brake mounted 


Fig. 8. 
on bridge drive of overhead travelling crane 
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brake coil is hot as when it is cold because, 
although the ultimate torque is less, the 
inductance to resistance ratio is smaller 
when the coilishot. It is evident that the 
time to develop torque is acceptably small 
for ordinary applications, and field tests 
have confirmed this. 

To obtain the torque-versus-time char- 
acteristics, strain gages were mounted on a 
connecting lever between the service 
braking armature and the inner shoe lever, 
and oscillographic records of strain versus 
time were made. Strain in the connecting 
lever is directly proportional to shoe force 
and, assuming constant shoe friction, 
to brake torque. 

At the same time that a strain pro- 
portional to torque was measured with 
strain gages, service braking coil current 
was also recorded. A  torque-versus- 
time curve was then derived from this 
curve as well as the statically determined 
relationship between torque and current. 
Agreement between the directly obtained 
torque-versus-time curve and the curve 
based on the current-versus-time curve 
was good, indicating that there is no 
appreciable time lag between current and 
the torque corresponding to this current 
under steady-state conditions. Also, it is 
clear that derivation of the torque-versus- 
time curve from the easily determined 
current-versus-time curve is acceptable 
for practical purposes. 


Parking and Power-Failure 
Braking 


The time to release the parking mech- 
anism of the brake is shown as a function 


of parking torque setting in Fig. 6. The 
middle curve shows the characteristic 
when rated voltage is applied to a hot coil 
and when the brake is adjusted for the 
recommended shoe clearance of 1/32 
inch. The other two curves show the 
extremes of the characteristic for a volt- 
age range of 80 to 110% of rated voltage 
and for a shoe clearance range of 1/32 
inch to 5/64 inch. 


Fig. 7 shows emergency set time as a 


function of parking torque for the 
rated voltage range of 80 to 110%. 


Emergency set time is practically im- 
dependent of shoe clearance since most of 


the time is required for initial release 
and relatively little time for movement 
of the linkage. 


Both release and set times of the 


parking section of the brake are compar- 
able to those of spring-set magnetic-re- 
lease brakes which are widely used, and 
thus 


these characteristics should be 


generally acceptable. 


Applications 


Fig. 8 shows one of two controllable 
torque brakes applied to the bridge drive 
of a manned trolley crane located in a steel 


mill shipping department. Here, accurate 


spotting of loads on trucks is accomplished 
by controlled braking. Since the opera- 
tor moves with respect to the bridge, the 
brakes are connected to the control 
through collector bars. 

The controllable torque brake has been 
applied to a coke pusher where accurate 
stopping of the pushing machine is neces- 
sary to line up a ram with an oven door. 


In this case, the auxiliary magnet is ¢ 
energized each time the pusher mach’ ; 
is lined up with an oven to hold t 
machine against movement while Ht 
pusher ram is passing through the oven 
Another typical field application i 
floor-operated crane controlled fron t 
pendant master switch. 


Conclusions 


Development of an electromagnel 
brake with controllable torque provic 
industry with a means to obtain smoq 
and accurate stopping of large drives w 
control from a small electric pilot devir 
Release of parking torque is accomplish 
quickly and with very little effort by 1 
operator. Substitution of electric Wii 
for hydraulic or air lines increases fl 
bility and reduces maintenance. FI 
formance predicted by tests at the ei: 
have been confirmed by apphcoaa 
the field. 
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Random Linear Systems: A Special Case 


ARTHUR R. BERGEN 


ASSOCIATE MEMBER AIEE 


N CONTROL SYSTEMS, elements 
are frequently used whose ‘‘gain’”’ is a 
function of variables outside of the 
ptimary control loop. For example, the 
gain of an armature controlled servo- 
motor depends on the motor field current. 
The gain of a hydraulic valve (with 
respect to flow rate as output and valve 
displacement as input) depends on the 
static head on the valve. The effective- 
ness of a control surface in turning a body 
moving in a turbulent medium varies with 
time: The gain from control surface de- 
flection to turning torque may be assumed 
to vary randomly with time. 
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Bergen—Random Linear Systems: 


It is interesting to know how a random 
gain affects a control system’s perform- 
ance. Even in those cases in which the 
external variables ffecting system gain 
are regulated, the effects of imperfect 
regulation are of interest. Thus, since 
regulators are ordinarily more effective in 
attenuating disturbances which are slowly 
varying rather than rapidly varying, the 
related gain variations in the primary con- 
trol loop may be assumed to be broad 
band in some cases of interest. This 
assumption, basic in the analysis which 
follows, suggests a tractable model for a 
random parameter system. 


A Special Case 


Analysis 


| 
| 
| 


The system considered is linear | 
time invarient except for a single 1 
domly-varying gain. Fig. 1s hows” 
system with the input and output of 
gain element exposed. The system 1 
not be a feedback system, but use iS 3 
case of interest. | 

The random gain k(t) is asceren t 
a stationary random process with meat 
and varying component a(t). 
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My 


JULY 


k(t) = k+a(t) 


ACTUAL SYSTEM 


9 u(t) = a(t) v(t) 


EQUIVALENT SYSTEM 


. 1. Equivalent representation of random 
gain 


) has zero mean. In Fig. 1 an equiv- 
nt system is shown; the variable y/(#) 
the equivalent system is the same as 
t of the original system. A new 
riable u(t) is defined which can be 
ated as a system input. Except for 
s dependent input, the system now 
pears to be a linear time-invariant 
tem. 

uuppose now the system error is to 
computed. Treating w(é) as a system 
ut, by superposition 


— ‘ae her(t—7) r(7) dr+ 
Ape woltey(t—7r) u(r) dr (1) 


‘0, this time, writing u(t) explicity as 


) a0) 


= Se bor(t—z) (7) det 
Jc heali=r) alr) u(r) dr (2) 


equations 1 and 2 /;;(¢) is the response 
the closed loop system at the point 
elled j to a unit impulse applied at the 

t labelled 7, with all other inputs zero. 
Phe solution of equation 2 is intractable 
use of the presence of the a(t) term. 
wever the mean-square error may be 
md readily if the following assump- 
is are made: 


The input vr(t) and the varying com- 
ent a(t) of the random gain are statisti- 
y independent stationary random proc- 


Gt) and a(t) are statistically independ- 
tandom variables for all ¢ not equal to 
(t) is a process with a constant spectral 


y. 

ho(0)=0, i.e., the application of an 
ulse at u does not cause a response at v 
antaneously. 


these assumptions the second is the 
trestrictive. Nevertheless, the prob- 
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lem is of interest since its solution (which 
is easily obtained) should suggest the types 
of behavior possible for systems in which 
the gain variation is not white but only 
broad band. 

With these assumptions it may now 
be established that E[a(é)o(é)r(4)]=0 
for all 4, E denoting the expectation. 
The reason follows. 

From equation 2, v(t) can depend only 
on past and present values of a/(é). 
But, by Assumption 3 only past values of 
a(t) can affect present values of (2). 
Then the present vaiue of v(f) is causally 
independent of the present value of a(d). 
Finally, since by Assumption 2 past and 
present values of a(t) are statistically in- 
dependent random variables, present val- 
ues of a(t) andy (¢) also are statistically 
independent. In this case 


Ela(t)o(t)r(t:) |= Ela(t)|Z[v(t)r(:)]=0 (3) 


since a(f) has zero mean. 

By the same token (assume f>#,), 
a(tz) is independent of all the other vari- 
ables and 


Elu(é a(t )u(te )a(te) =0ia 
E [v(t a(t )o(te )a(te) | = E [v2(t,) JE [a?(t)] 
b=h (4) 


Assuming that the spectral density of 
the a(t) process has a constant values So, 
it may be presumed on the basis of equa- 
tion 4 that 


Elo(t a(t )o(te)a(te) |= E[v*(t1) )Sod(t2—t1) (5) 


and this representation may be estab- 
lished. 

The expected value of v?(#) may now 
be found as a preliminary to finding the 
mean-square system error. Squaring 
equation 2 and using equations 3 and 5 


E[v?(t)]=02 +50 S — oheu2(t—7)X 
E[v*(r)|dr (6) 


the square and average of the first term 
of equation 2 results in the constant term 


denoted by w’?. It is a constant in- 
dependent of time since r(/) is assumed 


a stationary process. vw” is the mean- 
square value of v(t) if a(¢) is zero. The 
cross terms of the square of equation 2 
average to zero by virtue of equation 3. 
The square and average of the second 
term of equation 2 yields a double integral 
which reduces to the single integral in 
equation 6 by virtue of equation 5. 

Equation 6 may be solved directly for 
E[v?(t)] if it is known to approach a 
constant value in the steady state. This, 
however, may not be assumed a priort; 
the stationariness of 7(é) and a(t) do not 
insure this. 

It is thus necessary to investigate the 
limiting behavior of E[v?(¢)]. For this 
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purpose assume that although the input 
r(t) has been turned on sufficiently far 
in the past to obtain steady-state condi- 
tions by the time ¢=0, a(é) is only turned 
on at t=0. Then for positive time 
E[v*(é)] starts from the value m? and 
obeys the equation 


E[p%(t)]=00?-+ Sy ho2(t—7)Ev%(r) dr (7) 


Equation 7 differs from equation 6 in 
the lower limit of the integral and is 
appropriate for considering the transient 
behavior of E[v?(#)]. 

The solution of equation 7 is found by 
using the Laplace transform. Let V*(s) 
be the Laplace transform of E[v*(#)] while 
H™*(s) is the Laplace transform of h?,(é). 
Transforming equation 7 and solving for 
V*(s) 

Uo? 


ne) ~ s(1— SoH*(s)) 


(8) 
The behavior of E[v?(¢)] depends on the 
roots of 1— S,H*(s). If all the roots arein 
the left-half plane (excluding the imag- 
inary axis) there are transient terms which 
decay exponentially; E[v?(t)] tends to a 
steady-state value. If any of the roots 
are in the right-half plane E[v?(t)] 
diverges; the system is unstable. With 
roots on the imaginary axis, except at the 
origin which leads to divergence because 
of the double pole, an oscillatory behavior 
is possible. 

In Appendix J it is shown that a suffi- 
cient condition for all the zeros of 1— 
SoH*(s) to be in the left-half plane is 


SvH*(0)<1 (9) 


while in Appendix II this condition is also 
shown to be necessary. If 


SoH*(0)>1 (10) 


there is at least one root at the origin or 
on the positive real axis and the system is 
unstable in the mean square. Clearly a 
sustained oscillatory mode is not possible 
in this system. 

Assuming equation 9 is satisfied, a 
constant quantity v2 may be defined and 
evaluated using the final value theorem 
of the Laplace transform 


v= lim E[v*(t)] =———_— 
vic tet BOS Tee rea) 
U0? 


=—__-———_ (1 
1—iSi Ui hele ae 


This result can also be obtained (now that 
the steady-state constancy has been 
established) from equation 6 by assuming 


the limiting value v?. 


Equation 11 is valid as long as v? is 
positive. Nonrealizable negative values 
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wu(t)= a(t) v(t) 


for v? on the other hand imply instability ; 
see equation 10. 

Nowitis possible to compute the steady- 
state mean-square error. Taking the 
square of equation 1, noting equations 3 
and 5, and substituting the value of v 
found in equation 11 


22Sobo hheu2(r)dt 
(nS Phat de 


e-=ey (12) 
Here ¢e” is the mean-square error if a(#) is 
zero. This result is valid provided the 
denominator is positive. It may be noted 
that e? increases monotonically as a func- 
tion of Sy. A random gain always 
deteriorates the system’s performance. 

In the special case that v= Kye equa- 
tion 12 may be simplified. 


=~ — 1 
e? = 62 ——__——__ 
1 Sof.” Fou%(7)d7 


(13) 


Example 


A second-order servomechanism with a 
random gain is shown in Fig. 2. The 
equivalent representation for the random 
gain has been introduced. Assume the 
input is a rectangular wave with values 
+8, the input makes independent random 
crossings from one value to the other. 
The autocorrelation function in this 
case is! 


Oo ad a 


where v is the average switching rate and 
8 is the rms value of the input. If a(é) is 
zero (the nonrandom gain case) the sys- 
tem mean-square error is 


(K,+20)T+1 


PAN 3 
ih ae ae eer, 


(14) 
K,=KkKeis the velocity constant. _ This 
result is derived by Newton, Gould, 
and Kaiser.? Curves of the mean-square 
error are shown in Fig. 3 labeled Sp=0. 

Assume now the gain has a varying 
component a(t) which satisfies assump- 
tions 1 and 2 of this paper; Assumption 3 
is satisfied by the servomechanism under 
discussion. Equation 13 may therefore 
be used to evaluate the mean-square 
error. 

The computation of the integral in 
equation 13 is facilitated by using tables 
of integrals, for example, those in refer- 
ence 1. Assuming a positive velocity 
constant K,, the poles of Hyy(s) are in the 
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Fig. 2 (left). Servo- 1.4 
mechanism of the 
example 


rms /B 


.O 


Oo 
oO 


RMS ERROR/ RMS INPUT, e 


Fig. 3 (right). Effect 
of random gain on 00! 
system rms error 


left-half plane and there are also more 
poles than zeros; the integral in equation 
13 therefore converges and Parceval’s 
theorem may be used to evaluate it, i.e., 


dha hou®(r)dr 


jo 
1 
sees Hu(s)Houl —s)ds 
27J —jo 
Substituting 
—K, Ke 
Hyu(s)= Tete, 


and using the tables 


3 (KiK:)? Ke 
he r)dr = — = = 
i) tu (7) T 2K, QR? 


By use of equation 13 
aa lets: 

a= 15 
i SK, (15) 


Qk? 


é2 is given by equation 14. Curves of 
equation 15 are shown in Fig. 3. Note 
that there is now a minimum which occurs 
for 


2 
Kaas 
sya 
The larger the ratio of So/ k? the smaller 
the realizable K,. S)/k? is a measure of 


the relative variability of gain and is an 
appropriate normalization-of So. 


Conclusions 


The mean-square error of a linear con- 
trol system of arbitrary order subjected 
to a random input and in the presence 
of a single random gain within the loop 
may be evaluated by conventional tech- 
niques under the following circumstances: 


1. The input 7(#) and the varying com- 
ponent a(#) of the random gain are statistic- 
ally independent stationary random proc- 
esses, 


2. The random gain a(t) and a(t.) are 
statistically independent random variables 
for all f2 not equal to t. 
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0. 


clusions may be drawn. 

The mean-square error may converge 
a constant value or diverge but sus 
oscillations of the mean-square error 
not possible. The effect of the vail 
component of gain is to increase the 
tem mean-square error. Any app 
decrease of mean-square error is unrealil 
ble and indicates that instability in~ 
mean-square has occurred. 


Appendix | 


ue 


SoH*(0)=SoSo” h2(r)dr<1 


then 1—S,H*(s) has allits zeros in the] 

half plane. H*(s) is the Laplace transfd 

of h(t). 7 
Starting with the definition 


H*(s)=Jo° W2(r)e~ "dr 
= So? (ree 
Using the triangle inequality and : 
that on the imaginary axis or in the rigj 
half plane where a is positive or zero 
| H*(s)| <n r)e— “dr < 
So 12(1)dr =H 
Finally, introducing the hypothesis, 
So| H*(s)| <SoH*(0)<1 
on the imaginary axis or the right-half p! 
As a consequence 1—S)H*(s) cann 
zero in the right-half plane or ima: 


axis; all the zeros must be in the lef i 
plane. 


Appendix Il 
If 


SoH*(0)= Sy So h2(7)dr >1 


then 1—S,H*(s) has at least one rea f 
the right-half plane or at the origin. 


Joy 1 


ssume a real-axis zero of 1—S)H*(s) oc- 
;fora,areal value ofs. Ifs=aisaroot 


*(a)=1 
hypothesis 
0) >1 

1 


0) > H*(a) 


Noting that 


d @o 
~H"(a)= — f rh*(r)e-*%dr 


is negative for all a, and that H*(«) tends 
to zero as w tends to infinity, there is an 
a>0 for which S)H*(a)=1. This proves 
that 1—SoH*(s) has a real root at the origin 
or in the right-half plane if S)H*(0)>1. 


The sufficient condition for 1—.S,;H*(s) to 
have only left-half plane zeros is therefore 
also necessary. 
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transients. As a result of these diffi- 
culties, the classical method is to work 
directly with the open-loop function.1 A 
compensator is chosen to reshape the Bode 
or Nyquist diagram according to certain 
tenuous relationships between the fre- 
quency and time domain. For contin- 
uous control systems, the easy applica- 
tion of this method compensates for the 
uncertain final performance of the closed- 
loop system. 

However, in the case of sampled-data 
systems, the introduction of a compensa- 
tor, not employing a sampler, requires the 
calculation of a new z-transfer function for 
the compensator and the parts of the 
plant into which it is cascaded. Not only 
does this involve a great deal of additional 
labor, but the choice of the compensator 
and its effect on the open-loop transfer 
function is more difficult to predict than 
for continuous systems. 

These difficulties indicate that design 
of sampled-data systems should be based 
on a broader and more systematic pro- 
cedure than that described, and it should 
include elements of synthesis rather than 
trial and error, since the latter involves 
much numerical calculation for sampled- 
data systems. 


Systematic Design Methods 


A design procedure is shown in the 
form of a flow chart in Fig. 2. The var- 
ious blocks represent steps, and the 
letters beside each arrow identify the 
operation required to advance to the next 
step. The chart is headed by blocks rep- 
resenting available data at the start of 
the design: the time specification, and 
the plant. 
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Thus the procedure may be 


initiated from either block depending 
on the emphasis. If final time per- 
formance must be secured, the chart 
should be entered with the time-speci- 
fication block. If the plant involved 
must be preserved in its entirety, initia- 
tion should be from the plant-function 
block. In most cases, however, a com- 
promise procedure can be used. This is 
discussed later. 

If the chart is entered with the time 
specification, the specifications are inter- 
preted in terms of a closed-loop pole-zero 
configuration on the z-plane by means of 
correlation theorems, discussed in Appen- 
dix I. From the pole-zero configuration 
the open-loop root locus yields the corre- 
sponding open-loop pole-zero configura- 
tion. The term open-loop root locus has 
been introduced to distinguish this root 
locus method, which finds the open-loop 
poles from the closed-loop function, from 
the conventional root locus technique 
first introduced by Evans.?? This latter 
method of finding the closed-loop poles 
from the open-loop function is referred 
to in this method as the closed-loop root 
locus. The rules for preparing an open- 
loop root locus are prescribed in Appendix 
II. 

With the open-loop function which will 
satisfy the time requirements known, the 
designer can either: 


1. Proceed directly to determination of 
the compensator, or 

2. Attempt to conciliate the differences 
between the poles and zeros of the open-loop 
function and those of the plant so that, at 
least, the major poles and zeros of the plant 
are included in the final open-loop function. 


By the first method, the procedure is 
straightforward. The desired open-loop 
z-transform is converted into its equiv- 


A compensated sampled-data control 
system 


Fig. 1. 
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PLANT 
FUNCTION 


a 


CLOSED -LOOP 
POLE-ZERO 
CONFIGURATION 


OPEN-LOOP 
POLE-ZERO 
CONFIGURATION 


COMPENSATION 
FUNCTION 


e. CONVERT AND DIVIDE f 


f, NETWORK SYNTHESIS 
COMPENSATION 
NETWORK 


Fig. 2. Synthesis flow chart 


a. CORRELATION THEOREMS 

b. OPEN-LOOP ROOT LOCUS 

c. INSPECT AND SHIFT 

d. CLOSED-LOOP ROOT 
LOCUS 


alent Laplace transfer function. The 
transfer function of the compensator 
may now be found by simple division 


Ga(s) 
G,(s) 


where G(s) is the open-loop or “‘desired”’ 
transfer function, corresponding to the 
selected closed-loop function. The actual 
compensating network may be determined 
by standard network synthesis procedure. 

A study of equation 1 reveals that the 
compensator will, in general, cancel out 
each pole and zero of the plant, as well as 
add the poles and zeros of the desired 
function. In practice, complete cancella- 
tion of the poles and zeros of physical 
components such as electric and hydraulic 
motors is impossible, since their precise 
characteristics are never known over the 
full range of operation. Fortunately, 
system performance is not significantly 
changed as long as the pole and its 
cancelling zero lie near each other. This 


G(s) = (1) 


situation gives rise to a pole and zero . 


placed close together in the closed-loop 
function, and the resulting residue is 
negligible. 

A more serious objection to method 1 is 
the impossibility of cancelling out equip- 
ment time constants and expecting the 
equipment to function in a linear mode 
over a large range of input values. Can- 
celling out a motor time constant and 
replacing it with a much smaller one can- 
not be expected to result in an increase in 
the actual maximum obtainable accelera- 
tion, and the motor and gear combina- 
tion will not be able to follow a ramp 
which exceeds its maximum speed. For 
these reasons, the designer is strongly 
urged to proceed according to the second 
method. 

In method 2, the designer compares the 
poles and zeros of the desired open-loop 
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with that of the plant. By shifting loca- 
tion of the closed-loop poles-zeros, the 
major poles and zeros of the plant can 
usually be brought into the desired open- 
loop transfer function. The open- and 
closed-loop root locus furnish a graphical 
way of bringing about this correspondence. 
The impossibility of incorporating the 
major components of the plant into the 
desired open-loop function does not imply 
failure of this design method. Actually, 
this is a signal to the designer that the 
fixed plant is incapable of providing the 
performance demanded by the specifica- 
tions. This means that the plant must 
be improved or, if this is not practical, 
the specifications must be changed to 
concede the realities of the situation. 
Thus, method 2 increases the designer’s 
circumspection. 

If the compensator can be cascaded to 
the plant through a sampler, the pre- 
ceding steps may still be taken, and the 
compensator can be found directly from 


C(t) 


RISE 
TIME 
SETTLING TIME 


Fig. 3. Typical transient response to unit 
step input 


the z-transform of the system, that is, 


Ga(2) 
Gp(2) 


An alternate method is to start the 
design with the plant, first forming the 
open-loop pole-zero configuration and 
then the closed-loop root locus. The 
optimum performance possible without 
compensation can be determined by use 
of the correlation theorems. In general, 
additional open-loop poles and zeros 
will be required to satisfy the time re- 
quirement, their effect on the time re- 
sponse being ascertained readily by use 
of the root locus and correlation theorem 
technique. 


G,(z)= 


(2) 


Time Response and Basic Pole-Zero 
Configuration 


The critical step in this procedure in- 
volves choosing the proper pole-zero con- 
figuration. The typical time require- 
ments for a transient response to a step 


is shown in Fig. 3. The steady-state | 
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error to a ramp and parabolic inpu + 
shown in Fig. 4. Numerous pole-zex 
arrangements can be found to satisfy 
particular set of time specification 
However, the correlations theorems fi 
transient response listed in Appendix j 
are based on a dominant pair of comple 
poles. A pole-zero spaced close toge 
and near the point 1+ 70 in the z-pl 
provides control over the steady-stag 
error. With this, the basic pole-zex 
configuration in the z-plane is that show 
in Fig. 5. 

While all specifications cannot ~ 
satisfied by third-order systems, the sy¥ 
thesis can be handled initially in th 
way in many cases. As a last step in tlt 
design, more poles or zeros may have 3 
be added if the compensator is to 1 
physically realized. This is discusse 
in greater detail in Appendix III. If tt 
poles and zeros have large real parts - 
the s-plane, they will give rise to 2-plan 
poles and zeros near the origin and hay 
little effect on the previously determina 
results. Fig. 6 shows a typical ‘| 


\i 


zero configuration. 

Control of time response extends ont 
to the sampling instants. This general 
is satisfactory. However, a comple d 
description of the output between the ck 
instants may be secured by employity 


the modified z-transform technique 7 


the final system. i 

Plants with poles and zeros in the righl 
hand half of the s-plane require speci} 
consideration. The poles will ha 
images outside the unit circle in the | 
plane, as will the zeros in most case 
Rather than attempt to cancel these po! 
and zeros, they should be included 
part of the final open-loop transfer fv 
tion. 

The general methods of this present 
tion have been extended to sample 
data systems which do not possess — 
sampler in the error channel.> Thi 
class of system does not enjoy a close 


| 
| 


| 


Cc Cc i 
Ess | 

Es ae | 

TIME TIME 

(A) (8) | 

‘ 

Fig. 4. Steady-state errors to ramp — 


parabolic input 


A—Ramp input 
B—Parabolic input 


tee se se 


Jury 1 : 


p transfer function independent of the 
ut function. 


istrated Example 


As an example of using the proposed 
ign procedure, consider the design of a 
tem to certain specifications. 


The plant has a sampler in the error 
nel operating with period of 0.1 second 
1 a transfer function given by 


4 K(s+19) 
$(s+3.7)(s+8)(s+25) 


The time domain specifications are 


$0.35 


(3) 


wo items which must be decided be- 
e selection of the closed-loop pole-zero 
figuration are the difference in the 
mber of z-plane poles and zeros of the 
nt and how many of each are to be 
ained in the desired transfer function. 
the first item, the plant has one sam- 
t, no delays, and a numerator 3 de- 
es less than the denominator; hence, 
z-transfer function of the plant will 
e one more pole than zero. A dis- 
sion of the required minimum number 
z-plane poles and zeros and the re- 
red difference between the number of 
is given in Appendix III. 

or the second item, the retention of 
mt poles and zeros is frequently de- 
dent on the physical equipment. In 
§ case, an infinite value for K,* re- 
res the pole at the origin to be retained 
, arbitrarily, the pole at 3.7 will also 
‘etained. With these constraints, the 
plest pole-zero configuration can be a 
of complex poles and one real zero. 
p first choice of pole-zero selection is 


Z-PLANE 


5. Basic closed-loop pole-zero con- 
figuration 
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made with poles at 0.5—j0.45 and a 
zero at 0.3. Using the correlation 
theorems 


M=0.21 
Np =3.0 
Ns =8 

K,*=13 


The transient response requirements 
are met with sufficient margin to permit 
realization of a higher velocity error- 
constant by introduction of a dipole. 
First, however, the open-loop poles must 
be ascertained by the open-loop locus 
shown in Fig. 7. The gain is fixed by 
including one pole at point 1.0, and the 
remaining pole is found at 0.64. 

The s-plane pole of 3.7, which is to be 
included in the final open-loop function, 
corresponds to a z-plane pole of 0.691. 
Since the first try has produced a pole, 
not too far removed from this location, a 
closed-loop locus is now formed in Fig. 8. 


Z-PLANE 


Fig. 6. Typical closed-loop pole-zero con- 
figuration 


Choosing on the locus poles at 0.55+ 70.46 
gives the time specifications 


M=0.28 
Np =3 
Ns =9 
K,*=13 


With the gain determined, the open- 
loop function is 


0.591(z—0.3) 


(z—1)(2—0.691) (4) 


G(z)= 
A dipole is now employed to raise the 
steady-state error constant. But the re- 
quired close spacing of the pole and the 
zero will render graphical methods im- 
practical. An alternate procedure, avoid- 
ing this difficulty, works directly with the 
open-loop function. An open-loop pole 
at 0.99 and a zero at 0.95 increases the 
steady-state error constant to 65. 
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Z-PLANE 


Fig. 7. First selection for closed-loop 
configuration and open-loop root locus 


Realization of the open-loop function as 
a continuous one, without time delays, 
requires a zero at the origin. To retain 
the difference of one between the number 
of poles and zeros, an additional pole is 
added at 0.001. This pole was purposely 
placed near the origin to cancel the effect 
that the zero may have on the transient 
response. The desired open-loop system 
is thus described by the transfer function 


KA(z) 
‘B(2) 
0.5912(z—0.95)(z—0.3) 
~ (g—1)(z—0.99)(z—0.691)(z—0.001) 


(5) 


Inclusion of the additional poles and 
zeros will result in a shift of the complex 
poles and a change in the transient re- 
sponse. The new pole positions could be 
determined by analytical methods and the 
correlation theorems employed once again. 
But a simpler method would be to form 
the output z-transform function C(z) toa 
unit step directly from G(z) as 

a KA(z) 


sal KAG)EBG) (6) 


G(z)= 


C(z)= 


By simple division, the output at the 
sampling instant is, in this case, 


C(z) =0.591m(t —0.1)-+1.087u0(¢-2)-+ 
1.33u(f 3) +1.34%0(t —0.4)+ 
1.21210( + —0.5)+1.065m0(-+ —0.6)-+ 
0.988x0(t —0.7)-+0.935u0(t—0.8)-+ 
0.95510(t —0.9)-+0.995u0(t — 1.0) 


The transient specifications are now 
M=0.34 
Np=4 
Ns=9 


Converting the z-transfer function to 
its equivalent Laplace transform, with 
poles in the primary strip, is the next 
step. This involves taking a partial- 
fraction expansion of G(z)/z, multiply- 
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Fig. 8. Adjusted pole-zero configuration 
and closed-loop locus 


ing through by z, converting term by term 
to Laplace transform with the aid of 
tables, and, lastly, collecting the terms 


over a common denominator. In this 

case the final result is 
20.95(s-+0.512)(s+15.94) 

QS SS (7) 


5(s+0.1)(s+3.7)(s+69) 


The compensation is now determined 
by dividing the derived transfer function, 
equation 7, by the plant equation 3, 
yielding 


G,"(s) 
_ 20.95(s+0.512)(s +8)(s+15.94)(s+25) 
= (s+0.1)(s+19)(s+69) 
(8) 


The network is not realizable in its 
present form, but addition of a fourth 
pole far out on the negative axis has a 
negligible effect on the transient response, 
and yields a realizable compensation 
function 


G(s) 
_ 620(s+0.512)(s-+8)(s+15.94)(s+25) 
~  (s+0,1)(s+19)(s-+30)(s+69) 
(9) 


This function is realizable in terms of 
resistance and capacitors by standard net- 
work synthesis procedures. 


Conclusion 


The foregoing operation furnishes an 
integrated and systematic method of 
designing sampled-data feedback control 
systems to required time response through 
use of root locus techniques and correla- 
tion theorems. Its flexibility permits em- 
phasis on either the fixed components of 
the plant or on the time-performance of 
the closed-loop system. The method 
increases the designer’s circumspection 
into the effect of closed-loop and open-loop 
poles and zeros on the time response. 
This is particularly important in the 
design of sampled-data systems because of 
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the additional numerical work involved 
in taking the z-transform and the com- 
plication of determining a new trans- 
form function whenever a compensator is 
added without a sampler. 


Appendix |. Correlation 
Theorems 


If the transfer function is denoted 
KA(z) 
Biz) 


The closed-loop transfer function W(z) 
can be denoted 


KA(z) 


G(z)= 


(10) 


KQ(z) 


MO)" 5)+KAG) PC) fe 
or 
Wejo eae --(2—) (12) 


(2—pi)(—p2). . .(2—Pn) 


Then the main characteristics of a 
transient sequence of samples, as shown in 
Fig. 8, can be determined from the closed- 
loop poles and zeros with a predominant 
pair of poles, zo, as*® 


ny=2| Zang Q(z0)+ Pee) | (13) 


where 


imag (20) 
=tan-} ——~— 
pe real (29) 


M =(product of distance from all poles to 
the point 1+70, excluding distances from 
the two predominate poles to the point 1+ 
j0)/(product of distances from all poles to 
predominate poles excluding the distance be- 
tween predominate poles) X (product of 
distances from zeros to the predominate 
pole)/(product of distance from all zeros 
to the point 1+ 0) |z9| "2. 


log D 


af a log |20| 


(14) 


where D is the difference between unity and 
the accepted settled value, for a unit step 
input. (D is 0.05 in the illustration.) 

The steady-state error constants K*p, 
Ky*, and Kq* for a unit step, unit ramp, and 
unit parabolic input5 are in terms of the 
closed-loop poles and zeros. 


K-1 


K,*= ™ aa (15) 
H (1—-p1)-K Tg) 
i =e 
: (16) 
eS ae. i —9k 
t=1 
K,*= ee (17) 


SER 
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In terms of open-loop poles and zeros KH 
error constants are® 4 


K,*= lim G(z) ft 
2—1 q | 

ley | 

K,* == lim (2—1)G(z) (14 
T44 j 

iy 5 

Kq* == lim (z—1)?G(z) (20 
T* 2-1 | 


Open- Loop Ro ° 
ocus 7 
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For a system such as Fig. 1, the opes 
loop function C(z) and closed-loop functie 
W(z) are related by 7 

( 

G(z) 
1+G(2) 

Solving for G(z) gives 

W(z) .| 
¢e) = (2 
©) Twe) 1 


Hence, the poles on root of G(z) ange 


solution of the equation | 
W(z)=1 (21 


Alternately, equation 23 may be wri! 
as 


and W(z)=0+2hr 


W(z)= 


where # may be any integer or zero and. 
|W(z)|=1 


The approximate root locus genetag 
can be found quickly by applying 1 
similar to those for the conventional 
locus and by employing a spirule. 

Some helpful facts regarding approximal 
construction of the open-loop root loa 
are> 


1. The loci are continuous curves, star | 
at a pole for zero gain and terminating i 
zero or at infinity for infinite gain. 
2. The loci exist on a part of the real a d 
where an even number of poles plus ze | 
are found to the right. | 
3. If the number of poles w exceeds S 1} 
number of zeros m, there are w-m branch} 
which gotoinfinity. Further, the direct ih 
of the infinite asymptotes with respe ct | : 
the real axes are 


Qhr 
w—-m 


hao en a 


4. All asymptotes intersect at a poin 
on the real axis. , 


m w | 
» as, Ik | 
es c= 


w—m ii 

"a 

5. The angle at which a locus emanad 

from a complex pole or terminates ups Y 
complex zero is found by summing the at 

to all other poles and zeros. When emana 


4 


from a pole, the angles from the poles 8 
taken negative. When terminating on} 
zero, the angles from the zeros are tabi] 
negative. 
| 
| 
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Fig. 9. Various types of plants 


—Plant is continuous network 

—Plant includes hold network 

—Plant contains pure delay 

Parts of plant are separated by samplers 


Appendix Ill. Required 
Number of Poles and Zeros 


The general location of the poles and zeros 
achieve various requirements have been 
cussed. However, at the start of the 
thesis procedure, knowing the minimum 
ber of poles and zeros needed by a given 
plication is required. Although math- 
atically permissible, a practical compensa- 
cannot possess poles at infinity. Conse- 
ently, the excess of poles over zeros for the 
pensated or desired open-loop function 
st be at least equal to the excess for the 


plant Gy(s). Since the closed-loop function 
zeros are the same as those of the open-loop 
and the number of poles are the same for 
either open- or closed-loop function, it can 
be said that 


(number of poles—number of Zeros ) yz) 


(28) 
=(number of poles—number of Zeros) Gp) 


It would be advantageous to avoid the 
the lengthy process of taking the z-trans- 
form from the Laplace transform in order to 
determine the number of zeros of the plant. 
This number usually can be found by in- 
spection, as will be illustrated by consider- 
ing the open-loop systems shown in Fig. 9. 
When the continuous system, as shown in 
Fig. 9(A), has two poles in excess of the 
zeros, output will be available at the second 
sampling; ie, n=1. 

Since the value of the output can be deter- 
mined at the samplinginstant by dividing the 
denominator of G(z) into the numerator, an 
output at this sampling instant requires that 
the number of z-plane zeros be only one less 
than the number of z-plane poles. The 
number of zeros equals the number of poles 
when the Laplace transform poles exceed the 
zeros by only one. Also, in these cases, a 
zero is found at the origin. If a hold net- 
work is employed, as shown in Fig. 9(B), the 
same process of reasoning indicates that 
the number of poles exceeds that of zeros by 
one; however, a zero need not necessarily be 
at the origin. A plant including a delay, 
such as Fig. 9(C) of less than one sampling 
period, will give rise to the same z-plane 
pole-zero difference as discussed for Fig. 9(A) 
except the zero generally will not be found 
at the origin. 

For every complete sampling period that 
the delay involves, poles are placed at the 
origin without any additional zeros. The 
over-all transfer function for Fig. 9(D) is 
the product Gpi(z) Gpo(z), and the difference 
between the poles and the zeros of each func- 
tion is additive in the final transfer function. 


If the poles and zeros differ by one, then the 
over-all transfer function will have an ex- 
cess of two poles more than the zeros. This 
is in accordance with physical expectation, 
since the second network does not receive 
an input until the second sampling instant. 
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Discussion 


T. J. Higgins (University of Wisconsin, 
Madison, Wis.): The discusser has read 
this interesting and clearly written paper 
with profit. A rational means of designing 
compensating networks for sampled-data 
systems to meet prescribed response, rather 
than dependence on earlier used trial-and- 
error procedures, is most desirable. The au- 
thors are to be congratulated on their suc- 
cessful effort. 


The Synthesis of Optimum Systems 


from Nonideal Components 


J. K. WOLF 


NONMEMBER AIEE 


HE USEFULNESS of statistical tech- 
niques in the synthesis of optimum 
trol systems has been amply justi- 
d~7 and a voluminous literature 
ists on the subject. A wide class of 
‘h synthesis problems are those where 
basic circuit is fixed and the designer 
choose values for a finite number 
parameters. It is usually assumed that 
se parameters relate to components 
possess exactly their nominal 

es. In such cases, once a design 
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criterion has been established, the optimi- 
zation procedure is straightforward. 

It is the purpose of this paper to con- 
sider situations where the actual values 
of the components differ from the 
nominal values in some random fashion, 
i.e., the components have either associated 
tolerances or associated drift characteris- 
tics. In the following discussion, such 
components will be called ‘“‘nonideal.”’ 

If the components to be used are 
nonideal, then the optimization param- 
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eters can be treated as random variables 
whose probability densities have one or 
more variable moments. For example, 
consider the case where one of the 
parameters is represented by the value 
of a single resistor which has a +10% 
tolerance. If there is an equal proba- 
bility of the resistance R falling any- 
where within +10% of the nominal 
value R’ (and zero probability of falling 
elsewhere), the following density function 
would be assigned to the random variable 
representing the resistance value: 
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Manuscript submitted October 14, 1960; made 
available for printing December 16, 1960. 
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COMPENSATION 
NETWORK 


Fig. 1. Block diagram of servomechanism 


p(R)=5/R’, 0.9R’<R<1.1R’ (1) 
=0, elsewhere 


For this example, R’ is the mean value 
of the random variable R and is the 
variable moment specified in designing 
the system. In other cases, other mo- 
ments or characteristics of the density 
function might be used. The particular 
moment or characteristic which corre- 
sponds to the nominal value of the 
parameter to be specified will be called 
the zero tolerance variable (ztv). In 
general, all the moments of the density 
function could be functions of the ztv, 
but it will be assumed that the density 
function can be written as an explicit 
function of this parameter. For purposes 
of notation, the ztv corresponding to the 
parameter a; will be denoted by a’. 

If more than one parameter is involved 
in the optimization, a joint density 
function must be specified. However, 
in many cases, the tolerances associated 
with the different parameters will be 
statistically independent so that the joint 
density function can be written as the 
product of the individual density func- 
tions. 


Analysis 


In any statistical synthesis procedure, 
a criterion for optimization must be 
established. In accordance with the 
great majority of literature in this area, 
the minimum mean-squared error (mse) 
criterion will be used. 

It is necessary first to examine briefly 
the problem of designing a system with 
ideal components. The system will be 
assumed fixed except for the values of g 
parameters, ai, a2, ..., a. The mse 
can then be calculated in terms of these 
g parameters provided that the correlation 
functions of the input signals and dis- 
turbances are known. Formally, there 
is no difficulty in calculating the mse; 
in fact, the values of certain integrals 
which occur frequently in these calcula- 
tions have been tabulated. The opti- 
mum values, ayo, ag, ..., ago, Of the 
parameters ai, a2, ..., @g are then 
determined by standard minimization 
techniques. 

For systems with nonideal components, 
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it is assumed, as before, that there are 
q parameters to be specified and, in 
addition, that these parameters corre- 
spond to components with statistically 
independent tolerances. The density 
function for the 7th parameter a; will 
be written p(ai; ai’) where a;’ is the 
ztv of this density function. The average 
mse* is then 


ag )=S oe fe 


(Qi nen 


(au, oe 


q 

+» &q) II (ai; az’ )day.. .dag 
a 

(2) 


where ¢?(a1,..., a) is the mse that would 
result if the components were ideal. 
Thus, in equation 2, the average mse is 
used as the new optimization criterion 
and the statistics of the random param- 
eters must be taken into account as 
well as those of the signals and dis- 
turbances. 

The remaining procedure is, in prin- 
ciple, straightforward. The quantity 
€q” is minimized with respect to the ztv 
a1’, ae’, ..., ag’ by the usual methods 
of the calculus. The values of the ztv 
which produce this minimum will be 
denoted by aio’, 20’, hon Gao’ 


Equivalence of Solutions 


It is interesting to determine the 
conditions for which the minimum 
average mse choices of the ztv are 
numerically equal to the minimum mse 
choices of the parameter values for ideal 
components, i.e., the conditions for 


O40 = ain’ t= i 2, «++ (3) 
It is shown in the Appendix that, for 
components with statistically independent 
tolerances, the following conditions are 
sufficient for equation 3 to hold: 


1. The probability density functions 


are symmetric about their respective ztv. ° 


2. The mse (aq, . 
components is symmetric about its 
minimum (i.e., about the point-defined 
by Qi=Qio, =. 2; A on} q). 


.., @q) for ideal 


These conditions are sufficient but 
not necessary for the equivalence of 
solutions. 


Example 


The following example is taken from 
the control literature? where it first 


* Although the term “average mse’ may appear 
redundant, it is used to indicate that the squared- 
error is averaged over both the statistics of the 
random inputs and the statistics of the random 
parameters. 


G.(s)= 


appeared, with ideal components 
sumed. A simple compensation ne 
is to be designed for a servomecha 
The input 0;(¢) is a series of step dis 
ments of random amplitudes occu 
at random time intervals. The ra 
process defined by the time derivatiy 
this input is a white noise” with 1 
constant power spectrum Sp. | 

The block diagram of the sy. 
given in Fig. 1. In this figure, 4 
turbance 6,,(¢) is shown which is assui 
to be independent of the input signa 
to have the flat power spectrum . 
The motor is represented by the di 
tial equation 


J6(t) + Ko6)(t) aaa K;v(t) 


i 


so that its transfer function is 


Ks 
s(Js+K2) 


The compensation network has a Bi 
time constant 7; and an adjustable gi 
Ky; thus 


Ky 
ese 


Gi(s) — 


The problem is to specify the 
of K, which minimizes the stati 
expectation of the square of the | 
e(t). 

It is convenient to deaue the follow Ww 
parameters: 


aA(T,+J/K2)KiK3/Ke 


Ti(J/K2) 
Clie) 


LA No. 
~ So Ti +J/Ke)? 


K.A(M+J/K2)So/4 


If all the components are ‘ideal 
mse is shown® to be 


1+a(1—A)+La? 


Ke) Ee ea ae 


If @ is outside the interval (0, 
the system is unstable and the m 
infinite. All parameters are assur 
fixed except for K, which occurs linea 
in a. Thus equation 11 can be - m4 
mized with respect to a. The result 


ao= [4-+(A+L)/2]4 


which specifies the optimum value of 
for the ideal case. The mse? is 


@min= Ko[A +2(A +L)'7241] 
Let us now consider the non 


case and assume that a is a ra 
variable described by the density fur 


bla, «')=1/28, « '—BSa<a'+p 
=0, elsewhere 


: 
t 


Juuyg 


. 2. Optimum nominal gain as function 
of half-width of density function 


0% 


. 3. Percentage change in nominal gain 
function of half-width of density function 


Here the ztv 
The 


ere 6 is a constant. 
4s the mean of the distribution. 
erage mse is now 


ig SE +o(1—A ead 


.- ee pl meee f= Mies) 


(15) 


ich becomes after integration 


el itt) 
B 2e| \a2t 4)” 


Da Ale = 8)] 
[l—A(a’+8)] 


+loge a 
a — 8 


log, 


A 


ation 16 is a minimum with respect 
a’ when 
1 A-((A+L)+6%42?—A-L)*]'7 


= 17 
A?*—-A-L eke 


ao’ satisfies the inequality 


= 1 
ag’ Se (18) 
If inequality 18 is not satisfied, a 
ution may exist but the fact that the 
e is infinite outside the limits of 
bility must be considered. If 6 is 
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too large, there is no solution and an 
infinite mse results independent of a. 

For small 6, i.e., small tolerances, an 
expansion of equation 17 in a power 
series in 6 yields 


sone it _(4?-A-1) 

TAPE: ory? 
(A2-A-L) 
8(A+L)*/2 Br... (19) 


For 6=0, it is seen that ao’= ay. 

The results of the calculations in 
connection with this problem are 
shown in graphical form in Figs. 2 
through 4 for the case where the two 
time constants 7; and J/K, are equal 
(as assumed in the original problem") 
From equation 8 it is seen that, for this 
case, the parameter A has the numerical 
value 0.25. To consider reasonable 
values of L, the ratio Sp to No was taken 
to be 4.0 seconds~*. (The choice of the 
ratio So to No has been previously made 
by Newton, et al.1”) For this ratio and 
for typical time constants in the range 
0.05-1.0 second, the parameter L lies 
in the range 0.0625 to 25.0. 

In Fig. 2, the optimum nominal value 
ao’ of the parameter a is plotted versus 
the half-width 6 of the density function 
of equation 14, for L=0, 0.2, 1.0, 5.0, 
and 25.0. The optimum value ay’ 
increases with an increase in pf. It 
should be noted that the y-axis intercepts 
are the optimum values ap for the ideal 
case. Fig. 3 gives the same information 
as Fig. 2 except that the ordinate is the 
percentage change in the optimum value 
ao’ from its value for ideal components. 
The largest percentage change in ap’ for 
a fixed 8 occurs for the largest value of 
the parameter L. 

Fig. 4 shows, for the same five values 
of L, the percentage increase in average 
mse as a function of the tolerance param- 
eter 8. It can be seen that large 
tolerances (large values of 8) are re- 
quired to give appreciable error since 
the feedback system is relatively in- 
sensitive to gain variations. It should 
be kept in mind that the average mse 
has been calculated on the assumption 
that the optimum value of a, ie., ao’ 
given by equation 17, is used in the 
system. For any other value of a, 
e.g., ao given by equation 12, the average 
mse would necessarily be greater. 


- Conclusions 


A method has been presented for the 
minimum mean-squared error synthesis 
of systems built from components with 
associated tolerances or drift charac- 
teristics. The departures of component 
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Fig. 4. Percentage increase in mean-squared 
error as function of half-width of density 
function 


values from their nominal values have 
been treated as random variables. Syn- 
thesis is shown to involve not only the 
signal and noise power spectral densities, 
as in the usual case where the component 
values are assumed to be exact, but 
also the density functions characterizing 
component value variations. 


Appendix 


From equation 2, the average mse eg? 


can be written in terms of the mse e? for 
ideal components as 


+) aq’) = W ko A Ba ve 

we q 

CME a hes ay) plaj, x's | (20) 
j=l 


€q*(an', a 


Let the quantity g(a’, 
aq’) be defined by 


TOME =. pag yh flea 


qd 
Aa, . 5 o| pts w/e (21) 


jAt 


niet ME ee 


i Qa e 


Then, equation 20 can be written 


MHP ILE 


Qiy..., ay’ )p(ai, a; )day 


eq2(an', oe 
(22) 


The average mse will be a minimum with 


respect to the variables ai’, ..., ag, when 
de? 

= = 0) ME 1 ee (23) 
Oax" 


If, in equation 22, the partial derivative 


of €,2 is taken with respect to a,’, the 
result is 


foe) 

OE)? . ; 

ia CHS O ener oa unten) 
i —o 


ae [p(ai, ai’) da; (24) 
a4 


Now, if 2 is symmetrical (the symmetry 
must be even since e? is an extremum at 
axo) about the point ax =azo for k=1, ...,.q,. 


then g(a1’, ..., @i, ...; a’) has even 
g 


LoL 


symmetry about the point a;=aj) and can 
be expanded in the Taylor series 


(ai— ain) 
O(a een ee nag) = (2m)! 
m=0 
2m 
oe | (25) 
0a; =ci9 


From the binomial theorem, the term 
(ai— ayo)?” can be written 


2m 
(aj— ag)?” = > (ai— arg’) "( aes! — 
n=0 


The substitution of equations 25 and 26 
into equation 24 results in 


| |x 
Aj = Aig 
2m 


» et aa 
(2m—n)\n! 


n=0 


(oo) 


den? - og 
day! oe 0a;2” 


m=0 


eco Po) 
[tama 3a,” plas, a4!) dag 
(27) 


If plai, as’) is symmetrical about the 
point a;=a;', then the integral 


a r) 
[ (ai—ai’)” sat 


[p(ai, ai’) |dax 


is equal to zero for even values of m. Thus 
equation 27 can be expanded: 
Be? . 
=a(a;’ — ayo) +a3( a4’ — aio) + 
Oa,’ 
as( ai! — ajo)’ +. aie (28) 


It is apparent from this expression that a 
root of 


=0 (29) 


is aj=aj. In other words, the optimum 
value of a,’ We is written ajo’) is equal 
to aio. 

Since _ index 7 is arbitrary, it is true 
that aio’=az for i1=1, 2, ..., 9g, which is 
the desired result. 
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On the Design of A-C Servo 
Lead Networks 


GERALD WEISS 


MEMBER AIEE 


ARALLEL-T and bridged-T R-C (re- 

sistance-capacitance) networks are 
often used for the lead compensation of 
carrier-type servomechanisms. Presently 
available design methods require the use of 
cumbersome charts and do not readily 
permit evaluation of the relative gain. 
This paper carries out an analysis based 
on considerations of realizability. For 
each type of structure a realizability curve 
is determined. The analysis leads to 
equations relating the zero-frequency gain 
of the R-C networks to that of the R-L-C 
(where L is inductance) resonant damper. 
Finally, a very simple design method is 
derived, based on the conventional low- 
pass prototype, with no curves required. 


The Bilinear Compensating Function 


A great many feedback compensation 
problems can be handled by the bilinear 
transfer function or doublet 
pra 
p+b 


called a leading doublet for a<b and a 


G(p) =k$— (1) 
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lagging doublet for b<a. The function is 
easily realized by the simplest electric 
network, a single R-C L-section. For the 
lead network the maximum value of the 
constant multiplier & is unity; for the 
lag network it is b/a. Viewed in the 
complex plane, the lead network is used 
as a cancellation network; i.e., the zero 
at (—a) is used to cancel the most ob- 
jectionable real axis pole of the loop trans- 
mission. The lag network is used as a 
dipole close to the origin in type 1 control 
systems; it does not affect the root loci 
except in a minor way, but it increases 
the zero-frequency gain. 

For purposes of design, the frequency 
response point of view is probably more 
practical. At real frequencies p=jw, the 
network function becomes 


atjw 


oe es (2) 


The phase angle @ of G(jw) is related to 
frequency by 


w(b—a) 


tan hate RE A (3) 
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| 


The maximum phase Jee occurs at 
=V/ab 


and ee 


tan 0m = ioe == (hf) | 


The zero-frequency gain of the lead nq 
work is a/b, and the infinite-freque 
gain of the lag network function i is H j 
Hence the behavior of these network 
completely determined by the rea 
span a/b and the (geometric) cente 
quency +/ab. 
In the frequency domain, the s i 
lag network is basically a low-pass : 
: 


its ultimate purpose being to reduce 
gain at phase crossover without affec 
the latter. The phase lag of the nety 
interferes with this stabilizing ac 
Hence the attenuation must be ia 
duced at a very low frequency. | 
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Fig. 1. Comparison of lead networks on basis of gain 


I—“Optimum” low-pass prototype 


II—Modified low-pass prototype 


III—Approximation 


xample, consider a simple instrument 
srvomechanism with a mechanical time 
onstant of 1/30 second. The lag net- 
ork values of a and 6 should be about 3 
adians per second and 0.3 radian per sec- 
nd respectively. 

The lead network is used to reduce the 
Ihase lag at gain crossover; hence the 
alues a and 0b straddle the crossover 
oints. The compensation requirements 
re then usually couched in the following 
brim: We require @, degree of phase lead 
t @, radians per second. In other words, 
st one point on the frequency charac- 
pristics of the network is specified. 


Equation 3 then becomes a relationship 
etween the unknown network parameters 
and 6 and the specified 6; and «. 
is is one equation in two unknowns, 
ence the required 6, and w are satisfied 
y an infinite number of networks, e.g., 
es I and II of Fig. 1. All these net- 
orks have phase lead 6; at frequency , 
ut they differ in their maximum achiev- 
ble phase lead and in their zero-frequency 
ain. Out of the infinity of networks 
here is one (Fig. 1, curve I) which accom- 
lishes the desired phase lead in the most 
sonomical way, with the maximum zero- 
equency gain; this network has its 
aximum phase lead at the desired fre- 
uency. 6; and w are then given by 
quations 4 and 5 respectively, and these 
e the design formulas for a d-c lead net- 
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GiGjio) = 


jo+a 
jo+b 


jontai 
ed A Se 
Gj ») ae ie 


jotai 


G(ja) = 


work. Solving for 6; and w yields 


a 1—sin 6, 


b 1+sin 6, (6) 
r 1—sin 0, 7A 
== (0) 1) 
; cos 6; MES, 
1+sin 6 
ee (7B) 
cos 0; 


A-C Compensating Networks 


Compensating networks for a-c servos 
must have an envelope transfer function 


b-a 


(A) 


of the form of equation 1. It is readily 
shown! that a linear network cannot have 
such a transfer function. The familiar 
low-pass band-pass transformation 


CR PINGS 

yields a realizable network with transfer 

function 

Gas (9) 
S*-+-2bs-+w,? 

The corresponding envelope transfer 

function is a biquartic, closely approxi- 

mating the desired bilinear envelope func- 

tion. In the case of the lead network, 

a<b, the biquartic function has one 

dominant zero near (—a).? 

The function represented in equation 9 
is always physically realizable; an R-L-C 
realization can be immediately obtained 
from a low-pass prototype.* The prac- 
ticability of the realization depends on the 
effective Q, which is w,/26 or w,/2a for lag 
and lead networks respectively. With the 
use of a 400-cps (cycle-per-second) carrier, 
the previous example would require Q’s of 
the order of 4,000 for a lag network and 
40 for a lead network. Hence a-c lag 
networks are completely ruled out, and 
lead networks are the only passive a-c 
compensating networks used in current 
a-c servo practice. Note that the wider 
the servo bandwidth, the lower the re- 
quired Q, hence the easier the realiza- 
tion. These networks therefore find their 
main use in the compensation of wide- 
band servos; they ought to be called 
“wide-band” networks but, paradoxically, 
they are often referred to as ‘‘narrow- 
band” networks because the low-pass 
band-pass transformation is a “narrow- 
band transformation.” 

The usual R-L-C realization of the a-c 
lead network is shown in Fig. 2. An un- 


b-a 


NM 


(B) 


Fig. 2. R-L-C L-network realization 


A—Low-pass prototype 
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B—A-c network 


s?+-9as+-a,? 
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Fig. 3. 
A—Low-pass prototype 


p+a 
p+b 


G(p)=k 


grounded network with lower coil Q and 
smaller multiplier k, first described by 
McDonald,‘ is shown in Fig. 3. 

In quantitative terms, a useful measure 
of the required lead network bandwidth 
would be the ratio of servo bandwidth to 
the carrier frequency. The frequency at 
which the desired phase lead 6, is specified 
is of the same order as the servo half- 
power frequency. Hence it is plausible to 
use as a measure of network ‘“‘relative 
bandwidth” the ratio w;/w,. For the 
optimum lead network (maximum gain) 
w, is given by Vab. The following nor- 
malized quantities will be useful in the 
subsequent analysis: 


a 
ee zero-frequency gain (10) 


so SAU neg eee 


We 


(11) 


R-C Lead Networks 


Function 9 is R-C realizable when the 
poles are real, ie., for b>w,, which is 
equivalent to 


B>Va 


The carrier frequency w, must be larger 
than the lower corner frequency a. 
Hence the relationship between all the 
coefficients is 


(12) 


a<ae<b (13) 


This makes the zeros of G(s) complex, 
so that the function is not realizable as a 
simple ladder. If the function is mini- 
mum phase (a>0), the only case con- 
sidered here, the network design can be 
carried out by Dasher’s procedure® and 
the resulting structure is a 7-element 
bridged-parallel-T network. 
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Ungrounded R-L-C lead network? 


B—A-c network 


2521 Dastwe 

BO a connk a 
Ath (eg fay 
coil hey aye ar 


For the simpler structures in common 
use, the 4-element bridged-T and the 6- 
element parallel-T, condition 12, is no 
longer sufficient; larger values of b are 
required. The determination of the 
realizability conditions of these various 
networks forms the subject of this section. 
In the (a, 8) plane, the curve 


B=Va 


forms the boundary between the R-L-C 
and R-C regions. The boundaries of the 
bridged-T and parallel-T regions are de- 
fined by other equations in a and £. 
Finally, each specific structure, e.g., the 
symmetrical bridged-T, or the parallel-T 
with three equal capacitors, is character- 
ized by a realizability equation, an equa- 
tion relating a and 8 or, alternatively, a, 
b, and we. 

The literature dealing with these simple 
bridged-T and parallel-T networks is 
considerable,6 but apparently only 
Sobezyk’ and Benner® have considered 
realizability. The presentation in this 
paper is based on Sobezyk. 


(14) 


mR 


(A) 


Fig. 4. Bridged-T networks 
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BRIDGED-T REALIZABILITY 


The bridged-T networks under con) 
sideration are shown in Fig. 4. Sobez 
has derived the parameters m and 
comparing the transfer function of thes: 
networks with equation 9, term } 
term: 


4a(b—a) —w,? zd 
1S ee a (to 

We 

_. Meo | 
TE hah Over 5 


Pare 
~ 2(b—a) 


=RC 


positive; hence 


we?<4a(b—a) (1 
or 

1 > 

8>—_ == a 

24/1—0 Bs 


Equation 19 is plotted in Fig. 5, td 
gether with equation 14. Any constraay 
placed on the bridged-T network, 1.e., any 
specific relationship between the networ 
elements, defines a particular curve in th] 
realizability region. For example, tt] 
symmetrical bridged-T is defined by n= 


| 
3 t 
leading to : 


| 

wo? =2a(b—a) ¢ 
or 
f 

jo b | 
A/a) . || 


; 
iT 


This equation is also plotted in Fi g. 
The symmetrical bridged-T with eleme 
values is shown in Fig. 6. t 


PARALLEL-T REALIZABILITY 


The transfer function of the parallel 
network of Fig. 7 is a bicubic, rather thal 
a biquadratic. By comparing it termt 
term to the function 


nR R 
mC 


(B) 


JULY rk Lt) 


Fig. 5. Bridged-T realizability 
regions for 
0.9 
21 9 asta,2 
2] © oe j= 
me ones 
Gee 
ee = oundary 
“3 0.7 B=Va 
5 BRIDGED-T REGION a I|—Bridged-T boundary 
Q 0.6 
ae = “Nine 
a 2 q 
Be es 2V1-—« 
= RC REGION Ms “T I[I—Symmetrical bridged - T 
F 0.4| (NOT BRIDGED- T) curve 
2 p= 
P- LC REGION V2(1—a) 


_(s +245 +ee?)(s +N) 
(s?-++2bs +a,?)(s +N) 


obezyk has derived the value of the 
ircuit elements. 


(s) 


=RC=— (23) 
We 

_2(b—a)(2a-+N) 

* N2+2aN +0," 24) 
_e2[  2aN N +2aN +e, 

“5 [1+ @2 2(b—a)N | a) 
_2(b—2)(2aN +e062)N _ a 

1 (N2+2aN +0,2)eo,? ee 
N — N2424N +02 

ae cote 2(b—a) | os 


For realizability all parameters must be 
Ositive. This establishes two inequali- 
eS: 


(N) =N?—2(b—2a)N +[w,?—4a(b—a)] <0 
(28) 


we R 
az 


2 


0.3 


foN) = [we?—4a(b—a)] N2— 
2(b—2a)w?N +u*<0 (29) 


To meet these inequalities, the parameter 
N must meet the following conditions 


(b—2a) —V/b?—w,2< N<(b—2a)+ 
. Vb —w. (30) 

We? we? 
(b—2a)+-Vb2?—we? = (b— 2.) —-V/9—w,? 
(31) 


The biquadratic transfer function is 
realizable as a parallel-T if a parameter V 
can be found which satisfies simultane- 


ously the two inequalities 30 and 31. It 
can be shown? that for 
w,?—4a(b—a)<0 (32) 


a parallel-T realization is always possible. 
This is not an important case, however, 
since the simpler bridged-T structure is 
also realizable; see equation 18. 


(B) 


Fig. 6. Symmetrical bridged-T realization 
T=RC=a/w,? 
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Fig. 7. Parallel-T network 


For the case of 
@.*—4a(b—a)>0 (33) 


the parallel-T structure is realizable for? 


2V/a(b—2)<a,< (-2) ee 
(34) 
or 


a 1 


et ———— 
(0.5—a)+/0.25+ta—a? 2r/i—a@ 
(35) 


Fig. 8 shows the realizability regions. 

Realizability equations have been de- 
rived? for seven specific parallel-T net- 
works shown in Fig. 9; the first six of 
these networks are taken from Sobezyk,’ © 
and the seventh from White.? The var- 
ious equations are shown in Table I and 
plotted in Fig. 10, where the symmetrical 
bridged-T is also shown. 


Bandwidth Limitations 


Two facts emerge from the realizability 
curves of Fig. 10 and TableI: First, each 
network must obey a constraint between 
a, b, and w,. This is incompatible with 
the specifications, where the optimum a 
and b are determined directly by the 
specifications 0, and w; (equations 6 and 7); 
the carrier frequency is independent of 
them. 

Second, these specific constraints call 
for very large values of b and B. Con- 
sider the case a=0.1, which corresponds 
to a maximum lead angle of 55 degrees; 
this is a very common specification. 
R-C realizability alone requires a 6B of 
0.32. Realizability by a symmetrical 
parallel-T requires a 6 of 0.61, and a sym- 
metrical bridged-T a 8 of 0.71. The 
latter figure would mean a 400-cps servo 
with crossovers in the vicinity of 300 
cps. Such wide-band servos are never 
encountered in current practice. 

It is concluded that R-C lead networks 
are not practically realizable if one insists 
on minimum zero-frequency attenuation. 
They can be realized, however, if one 
starts from a new low-pass prototype: 


jotta, 


Gi(jo) OED: (30) 
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RC REGION 
PARALLEL-T 


Fig. 8. _ Parallel-T 
realizability regions 
for 
+2%astw,? 
2+ 9bs+u¢" 


I—R-L-C/R-C 


G(s) = 


boundary B= Va 
||—Parallel-T 
boundary 


RELATIVE BANDWIDTH 


as shown in Fig. 1. The value of b; can 
be made sufficiently large so that realiza- 
bility by the desired network is achieved. 
The price paid is increased zero-fre- 
quency attenuation. 

The transfer function of the a-c net- 
work is now 


s?+2a15 +e? 
s?+2b15 +a? 


New, normalized parameters a, and f; 


G(s)= (37) 


(B) 


(c) 


A—Symmetrical parallel-T T=1/w, 


B—Degenerate case, three equal C’s, two equal R’s 
C—Degenerate case, three equal R’s, two equal C’s 


T= V2/we 
T=1/V/ 20 


D—Three equal capacitors 
E—Three equal resistors 
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Va 

(0.5—a)+ 
/0.25+a—a* 

|1|—Bridged-T 


boundary 
4 


Soi 


B= 


can be defined by reading equations 10 
through 35, the relationships in Table I, 
the element values in Figs. 6 and 9, and 
the co-ordinate in Figs. 5, 8, and 10, in 
terms of a1, b;, a, and fi. 


Determination of Parameters 
The problem now is to derive the values 


of a, and 6b, from the original specifica- 
tions, the phase lead 6, at frequency a. 


Fig. 9. Various parallel-T networks 


F—Full-rejection symmetrical parallel-T, by- passed with input div of 


T=1/2a 
T=2a/w,? 
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We start with the fwo low-pass proto. 
types of Fig. 1, the optimum prototype — 


ih jotta s 
ee St 
and the modified prototype Gi; of eq Ta 


tion 36. For G, the point (4, a) | i 
merely one point on the curve; thus 


i. 
| 


For G the point (#1, «) is the point 
maximum phase lead; thus 


tan 6,= = =A 1-2) by 


and j 
w= Vad (4 


Substituting equations 40 and 41 inti 


equation 39 yields % 
Vab(bi= a1) _ b-—a ( | 
a1b,—ab 2/ab | 


Equation 42 must be solved simultang 
ously with the realizability equation i 
the particular structure, and this giv 
relationship between the pairs (a, bd) io 
(a, 61). The equation so obtained x 
cubic. To sidestep the difficulty of sat 


R R @ 


iN 


ae - ig a En EE a ant 


/] 


T=1/e; RaXZin 


G—full-rejection symmetrical parallel-T, by-passed with brid i 


resistor? T=1/w, 


Y 


PSS EAE ANS 


Baier! 


RELATIVE BANDWIDTH 


A 

: 
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. 

é 

‘ 
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Fig. 10. Realizability curves for various R-C structures 


l—Symmetrical bridged-T 

II—Degenerate case, three equal C’s and two equal R’s, or three 
equal R’s and two equal C’s 

I]|—Parallel-T, three equal C’s or R’s 

IV—Symmetrical parallel-T 

V—Full-rejection circuit with input divider 

V|l—Full-rejection circuit with bridging resistor? 

VII—R-L-C/R-C boundary 

Vill—Parallel-T boundary 

IX—Bridged-T boundary 


Table |. Parallel-T Realizability Equations 


Condition on 
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Fig. 11. Gain-phase curves for 55-degree phase lead at 15 cps, with 


400-cps carrier 
A—Asymptotic gain curves 
B—Phase curves 


I—R-L-C network 
!|—Symmetrical parallel-T with unmodified lower corner frequency 
I1I—Symmetrical parallel-T, final design 


ing the cubic, consider’ first a simpler 
equation obtained by assuming bi>>5), 


so that 

Gee (43) 
by 

tan 6; eet (44) 
ay 


Reieapine The corresponding polar plot is the dashed 


Parameters, _Circuit, Curve, line of Fig. 1. Substituting equations 
Network Fig. 7 Fig. 9 Realizability Equation Fig. 10 40 and 41 into equation 44 yiel ds 
| = b-a ‘ab 
=m=1 A = ave IV Pela (45) 
Ss psal'e a) aual{e(s ejisiia/aei\s) e) = THY TA ec eevee fhe ves escs (1—3a)+/1+10a—7a? rare /ab ay 
1 or 
qual C’s and 2 equal R’s..... m=ne=1........ EB ore oievopers Baa i=on SSCS OPO RICE Cer Arete II 
Pane 2ab 2a 2a 
a= = = (46) 
1 b-—a a 1l-—e 
qual R’s and 2 equal C’s..... m=me=1....... CE comin B POUT ae Piow re aware, cor etaae 3 II ae 
m1 =ne2 
2 3/2 = for 
Shobedcien Guadss Coe Ons m=ng=1........ Datel 2 a are RETEST sarciorete : 
— 1 =2. 
E B By 28 III >: Cs ae 
Da OHATEIGTESE CLANS Oe m=mg=1.......E........ SS sn Hi Oi 
1—3a)+V1+18a—15a? 3 é 
is = Equation 47 yields the remarkable re- 
mejection, input divider........6002.seeeeeee Uiewonopae p=2V/ a. sie petinccs a chy os Vv Slt chatcthe lower Gaines frequency of 
2Va S : 
rejection, bridged®..............000eeeeeeee: Glacaee B EVE the optimum low-pass prototype ought to 


be doubled. The physical significance 
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is this: With the optimum prototype 
a set of values a and 6 can be found for a 
given maximum phase lead 6, at the 
modulating frequency w;. If one now 
increases the upper corner frequency b 
merely to obtain R-C realizability, the new 
phase lead at the frequency w will be 
somewhat larger. This excess of phase 
lead is not harmful, but it is obtained at 
the cost of additional attenuation, the 
zero-frequency gain being a/b. The sur- 
plus phase lead is traded in for gain by 
increasing the lower corner frequency a. 

As a numerical example, consider a 400- 
cps carrier-frequency servo to be compen- 
sated by 55 degrees of phase lead at 15- 
cps modulating frequency. The corner 
frequencies of the optimum low-pass 
prototype are at a=4.74 cps and at b= 
47.4 cps, for a zero-frequency gain @ of 
0.1. If the value of the upper corner fre- 
quency 0 is now changed to 800 cps for 
symmetrical parallel-T realizability, the 
new phase lead at 15 cps will be 71.1 de- 
grees, much larger than required. By in- 
creasing the lower corner frequency a to 
9.48 cps the phase lead at 15 cps becomes 
55.7 degrees, close to the desired value. 
The corresponding gain and phase-versus- 
frequency curves are shown in Fig. 11. 

The results obtained have been based 
on a simplified analysis using}:>>b. A 
more exact analysis requires simultaneous 
solution of equation 42 with the appro- 
priate realizability equation. The result- 
ing sets of cubics have been plotted and 
analyzed. For servo parameters en- 
countered in current practice, say a<0.25 
and 8<0.25, one again obtains the ap- 
proximate relationship a,= 2a. 

It is also to be noted that the phase lead 
calculations are based on the low-pass 
prototype, not on the exact envelope 
transfer function. This approximation 
yields slightly optimistic results. The 
true envelope lead is somewhat smaller 
than the number indicated by the low- 
pass prototype, but the error is less than 
0.5 degree for the modulating frequency 
ranges used in current servo practice.? 


Relative Gain 


The R-C networks are obtained only 
at the expense of zero-frequency gain. 
It appears highly desirable to be able to 
determine the relative gain of various 
types of network structures in advance of 
carrying out the detailed network design. 
This information is obtained from the 
ratio ai/a, the ratio of relative gain for a 
particular R-C network to the maximum 
relative gain obtainable with the R-L-C 
network. The ratio a;/a can again be 
obtained in a rigorous way by analysis of 
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the cubics resulting from simultaneous 
solution equation 42 with the realizability 
equations. A simpler and equally valid 
procedure is to solve the realizability 
equations simultaneously with the use 
of equation 47. 


BrRIDGED-T 


From equation 18 the realizability 
equations for the family of 4-element 
bridged-T networks is 


Wo” — ha;(by —a) (48) 


where h<4 for bridged-T realizability; 
for the symmetrical bridged-T h is 2. 
Rearranging this equation yields 
ha? 
=——— 49 
ha, Hee! Sed 
Substituting equation 47 into equation 
49 yields 


4ha? Aha , 
a es 50 

apraeacnie prema 2 a) 

Hence 

& =4hp? (51) 

For the synmetrical bridged-T 

= =88? (52) 

(34 


PARALLEL-T 


In Table I and Fig. 10 all parallel-T’s 
except curve II (which behaves similarly 
to the bridged-T and is therefore of no 
practical interest) are characterized by the 
approximate realizability relationship 


hB =o (53) 


where h<1 for parallel-T realizability; 
for the symmetrical structures h is 1/2 
and for the three equal R’s and C’s, h is 
4/2/3. Equation 53 can also be written 


we =hbr (54) 

Combining equations 47 and 54 yields 
2h = 

Aa 9/56 (35) 

by We 

Hence 

a, 2hg 

a Va 2 

and for the symmetrical parallel-T’s 

CH fh 

meee | e) 

CoMPARISON 


The relative zero-frequency gain of 
R-L-C versus symmetrical parallel-T 
versus symmetrical bridged-T is obtained 
from equations 57 and 52: 
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1 to E to 86? (58) | 
As a numerical example, take the case 
of a 400-cps servo, with a 55-degree phase 
lead required at 40 cps. This is equiy- 
alent to a=6=0.1. The relative gain : 
of the three networks is 12.5 to 4 to 1 
The lower the bandwidth, the more n- 
favorable the ratio for the R-C networks 
but, of course, the larger the Q of thes 
R-L-C network. The Q for the com} 
ventional resonant damper which is given 
gain, is 
1 7 | 
=o = 59) ) 
ee 5: 2V/ap (59) ) 


For the above numerical case, the Q is& 


16 | 
Network Design - 
§ 


The network design requires the follow- 4 
ing steps: : 


é 
* || 
1. Calculate a and 6B from the speci: | 


fications. $ 

2. Determine the relative gain of the d 
various types of structures using equation 4 
58. Determine the Q of the R-L-C ee | 
using equation 59. 


decision to be based on results of steph 28 
plus other factors, such as availability ands 
price of components, structural factors,} 
environmental factors, etc. | 


For the R-L-C network the remaining: 
steps are: ay 


4. Determine a and 6 from the specifica-} 
tions, using equations 6 and 7. 


5. See network element values in Fig. 27} 
within a constant multiplier. 


6. Determine desired impedance level! 
and multiply all element impedance by theé 
appropriate factor. : 


For the R-C networks, the remaining ste pss 
are: | 


4. Determine a, from the specification, 
using equations 7(A) and 47. 
5. See appropriate network element values 
in Fig. 6 and Figs. 9(A), (D), (E), (F), andi 
(G). Element values given within a con 
stant mere 


constant. Equations for impedances 
R-C networks are given by Sobczyk.? 


Concluding Remarks 


ratio”’ ae to @ in this paper) ai ane i 
“notch width.” These methods 


ther cumbersome and do not permit 
lvance determination of the relative gain 
hievable by the various structures. 
sing the analytic methods here pre- 
nted, an equation for relative zero- 
equency gain is derived, for R-L-C 
rsus parallel-T versus bridged-T. A 
cision on which network is to be used 
i then be based on all the pertinent 
chnical and economic factors, without 
le need for extensive tentative designs. 
nce this decision has been made, the 
‘tual network design is almost trivially 
mple: from the specifications one de- 
gns a low-pass prototype as one would 
ra d-c servo, calculating the lower corner 
quency only. This value is doubled. 
he upper corner frequency is not needed; 
is automatically determined by the 
Oice of structure. The accuracy of the 
sign procedure is well within that re- 
aired by the usual phase lead specifica- 
n. 

The symmetrical networks considered 
this paper have realizability curves 
nich lie well within the interior of the 
spective realizability regions; see Figs. 
and 10. It might be argued that asym- 
trical networks located closer to the 
undaries are superior. Equations 51 


Discussion 


Iph J. Kochenburger (University of 
pnnecticut, Storrs, Conn.): This paper 
tforms an excellent service, both of 
mmarizing methods of carrier compensa- 
m and of proposing new and practical 
teria that might be used in design. The 
ngrams shown in the author’s Figs. 5, 8, 
d 10 are of special interest to the designer 
ause they point out the regimes of opera- 
m where certain specific types of carrier 
mpensation networks will prove to be 
2 most practical. 

his discusser wishes to call attention to 
me additional considerations that should 
kept in mind by engineers planning to 
p these techniques: 


The author refers to the common 
p of phase lead networks as cancellation 
works. This viewpoint is a common 
2 and has the advantage of providing 
ople mathematics once an objectionable 
e has been cancelled and replaced by 
bther pole of lower time constant. How- 
er, the optimum selection of the com- 
sating zero (the 0} in the author’s nota- 
in); i.e. of the low-frequency break point, 
not always correspond to exact can- 
lation. Very often, it is found that a 
ver b (i.e. higher low-frequency break 
nt) than that corresponding to pure 
icellation will be desirable. 


In the author’s brief discussion of 
ase lag compensation he states that such 
npensation has a minor effect on the root 
. Actually, the root loci are modified 
such networks in a very significant way. 


LY 1961 


and 56 indicate a possible doubling of rela- 
tive gain for asymmetrical networks. 
Unfortunately, however, the input and/or 
output impedances of the networks become 
zero or infinite on the realizability bound- 
ary.” Hence the actual transmission of 
these asymmetrical networks is likely to 
be smaller than that of the symmetrical 
networks. This does not show up in 
the present analysis where both gen- 
erator and load impedances have been 
neglected. 

The assumption of zero source imped- 
ance and infinite load impedance is not an 
unreasonable one when the actual load 
impedance is at least one or two orders of 
magnitude larger than the source imped- 
ance. This condition is not always met 
in practice; for example, many transistor 
amplifier and magnetic-amplifier circuits 
have low input impedance. An exact 
analysis taking into account both source 
and load impedance can be expected to 
show that the zero-frequency transmission 
is maximized by using an asymmetrical 
network. However, the optimum is not 
likely to be very sharp, and the symmet- 
rical structures considered in this paper 
have at least the merit of simplicity and 
economy. 


4 


The low-gain portion of the loci are swung 
toward the imaginary axis; a conditional 
degree of stability is caused; and if the 
designer fails to take advantage of the 
gain increase afforded him by the use of 
such networks, a slow response with very 
poor damping will result. 


8. Again, with reference to phase lag 
compensation, the difficulties stated by the 
author are true but are not the only reasons 
why carrier compensation is impractical 
for this purpose. A very important con- 
sideration is that of carrier-frequency shift. 
Phase lag networks are especially sensitive 
to such shifts because of the very small 
notch widths involved. Notthoff! has 
shown how carrier shifts can cause poor 
operation when they become an appre- 
ciable fraction of the notch width. For 
example, a typical notch width associated 
with phase lag compensation rarely would 
exceed 0.1 cps; 400-cps carrier sources 
with a frequency regulation better than 
0.1 cps are “few and far between.” 


4. Probably the most significant point 
that should be added to the author’s state- 
ments is that operation at high relative 
bandwidths is undesirable even though, 
paradoxically, such operation would present 
the least difficulty from the standpoint of 
network synthesis—because of the data- 
sampling effect of carrier data transmission. 
This effect is somewhat similar to discrete 
data sampling although it is more difficult 
to analyze. As in the instance of discrete 
data sampling, useful control information 
cannot be transmitted when the signal 
frequency becomes too high a fraction 
of the carrier frequency. In practice, a 
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B of 0.3 is a good practical upper limit. 
For example, if a servo bandwidth of 300 
cps were desired, 400-cps carrier would 
no longer be adequate and at least 1 ke 
would be recommended. Because of this, 
the practical regimes of carrier compensa- 
tion will fall in the region of B<0.3. From 
the author’s conclusions, this would corre- 
spond to the recommended use of R-L-C 
and parallel-T forms of compensating net- 
works. 


5. Even phase-lead carrier compensa- 
tion networks are subject to carrier fre- 
quency shift, although not as severely as 
the phase lag networks. Superposed on 
this difficulty are the other practical design 
difficulties mentioned by the author. The 
necessity of finding very high Q inductive 
elements when R-L-C schemes are em- 
ployed and of accepting excessive attenua- 
tion at zero frequency when parallel-T 
and other R-C networks are employed 
are important disadvantages. For this 
reason, I must admit that I have developed 
a prejudice against carrier compensation 
in general. Fortunately, in most problems 
encountered, it has been possible to avoid 
the issue. Usually some method has been 
available whereby, without the introduction 
of excessive new hardware, noncarrier 
compensation of both phase lag and phase 
lead types can be introduced in a-c servos. 
Most often this has been accomplished by 
feedback methods of compensation. 
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Gerald Weiss: Dr. Kochenburger’s compli- 
ments are gratefully appreciated, and his 
additions and corrections are very much to 
the point. 

Inasmuch as the effect of carrier frequency 
shift on these networks depends on their 
notch width, a large notch width becomes 
desirable. In that respect, then, the 
bridged-T network is the best, the parallel- 
T network follows, and the R-L-C damper 
is the poorest, that is, the most sensitive 
to variations in carrier frequency.1 This 
surprising result can be rigorously verified 
by envelope transfer function analysis.? 

An additional interesting and practical 
problem is the design of an R-L-C damper 
with a low Q coil. Consider the example 
in the paper, where a=6=0.1 with a 400- 
cps carrier; the Q required for the optimum 
network was found to be 16. Now assume 
that available inductors have a Q of only 
12. It ought to be possible to design a 
resonant damper with this coil, and accept 
some extra attenuation. The alternative 
would have to be a parallel-T with relative 
attenuation of 3.2. 

Realizability of the lower coil Q requires 
an increase in the lower corner frequency 
a [the shunt resistor in Fig. 2(A) of the 
paper]. This reduces the available phase 
lead, and the upper corner frequency 6 has 
to be increased by a larger ratio in order to 
obtain the desired phase lead at the specified 
frequency. Note that this is the opposite 
procedure from the R-C network design 
where the upper corner frequency 0 is 
increased for R-C realizability, and then 
the lower frequency a is increased to trade 
in surplus phase lead for gain. 
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Fig. 12 (left). | Realizability 
region for reduced coil Q 


Fig. 13 (right). Relative gain 
versus relative coil Q for reso- 
nant damper 


This problem is solved by starting out 
with the relationship among the coefficients 
of the optimum and modified prototypes, 
equation 42 of the paper, 


bi —a 
ab— aby 2ab 


(60) 


One can then define two values of coil Q, 
the value Qo for the optimum network 
(maximum gain), and the value Q,; of the 
available coil: 


On 5 ak (61) 
oN arm (62) 
Qi_@ 

Os an (63) 


Substituting a: and 6 from equations 61 
and 62 into equation 60, one can determine 
the relative gain a,/a and also the relative 
value of the series resistor [see Fig. 2(A) 
of this paper], (b:—a)/(b—a). 
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These quantities must always be posi ti 
hence one obtains the realizability conditio: 


Qo 
which is plotted in Fig. 12, where 1 
relative coil Q for half-gain is also showr 
Equation 64 is also shown in Fig. 13 fo 
a=0.1. The diagrams indicate that ; 
reduction in coil Q up to 50% can 
achieved. 

I concur emphatically with Dr. Kochen 
burger’s strictures against carrier le 
compensation. In my opinion it is the le 
compensation per se, rather than the ca 
compensation, which is the main villai 
Lead networks increase the system bane 
width, require high amplifier gain, 2 
render the system more vulnerable to noise 
The resulting bandwidth is often far i 
excess of that required by the applicatior 
This weakness of the lead compensatio 
scheme is accentuated in carrier-frequene 
systems because of the presence of add: 
tional high-frequency noise in the form | 
carrier frequency and harmonics. In add 
tion to this, the R-C carrier lead network 
require a further increase in amplifier gai 
with concurrent reduction in signal+t 
noise ratio. d 
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